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ANALYTICAL AND EXPERIMENTAL INVESTIGATION 
OF HIGH ENTRAINMENT JET PUMPS 

by Kenneth E. Hickman, Gerald B. Gilbert, and John H. Carey 


SUMMARY 


The use of jet pumps is of increasing interest for boundary layer control 
or control force augmentation in V/STOL aircraft. In typical applications, a small 
mass flow of primary air at pressures up to 400 psia can be used to entrain a much 
larger mass flow of secondary air at ambient conditions. The primary nozzle flow 
is supersonic while the secondary flow is subsonic. The jet pump system design ob- 
jectives may be maximum entrainment, maximum thrust augmentation, or some 
combination of the two. Little information is available in the literature to guide the 
designer of jet pumps for such applications. 

In this investigation, a simple analytical model was developed to predict 
the performance of high-entrainment compressible flow jet pumps with constant area 
mixing tubes. While the model is suitable for hand calculation, a computer program 
was prepared to facilitate calculation of jet pump performance curves and allow com- 
parison of different jet pump designs. Analytical techniques were developed for 
matching the jet pump design to its associated duct system in order to achieve maxi- 
mum entrainment or thrust augmentation. 

The validity of the analytical model was confirmed by an extensive test 
program using a multiple-nozzle jet pump with two different mixing tube lengths. 

The primary-to- secondary flow area ratios were varied from 0. 0013 to 0. 0067. 

The primary flow pressure ranged from 55 psia to 350 psia and the primary flow 
temperature ranged from 200° F to 1200° F. The observed entrainment ratios varied 
from 15 to 37. The performance of each jet pump geometry was measured over a 
very broad range of operating conditions in order to develop performance maps for 
comparison with the analytical predictions. 


Section 1 


INTRODUCTION 


1. 1 Background 

Jet pumps have been used for many years in industrial applications where 
a high-pressure gas such as steam is used to pump a lower-pressure gas. The jet 
pump is a simple low-cost device with no moving parts and is particularly convenient 
for use with troublesome fluids such as two-phase flows, high-temperature gases, 
or corrosive gases. Jet pumps are usually employed as low-pressure-rise devices 
and their thermodynamic efficiency is low, i.e. , under 20%. Because they are low- 
cost devices of limited performance potential, there has not been a strong incentive 
for research and development work on industrial jet pumps. 

In recent years, applications of jet pumps to boundary layer control sys- 
tems have become of increasing interest for STOL aircraft. Systems have been pro- 
posed which use jet pumps to entrain a large flow of secondary air which is then 
directed over a deflected flap for lift augmentation. In a configuration patented by 
F. G. Wagner (references 1 and 2), a jet pump is used to entrain air from one sec- 
tion of the trailing edge of a wing (boundary layer suction upstream of a deflected 
flap) and then to discharge it over a deflected flap. In this way, the inherent ineffi- 
ciency of the jet pump is partially balanced by the double employment of the entrained 
air for boundary layer control. Jet pumps may also have application in VTOL air- 
craft for direct lift or control force augmentation. The primary, high-pressure flow 
for the jet pumps can be provided by a bleed from the main engine compressors or by 
an auxiliary power unit. 

The use of jet pumps as primary components of V/STOL aircraft sys- 
tems places new emphasis upon development of design techniques for these devices. 

It is essential to be able to minimize the size of jet pumps for particular primary 
and secondary flow conditions, and to be able to predict the performance of jet pumps 
over a broad range of operating conditions. However, systematic design and analy- 
sis procedures are not available for high-entrainment-ratio compressible-flow jet 
pumps. 
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1.2 Previous Work 


A number of investigators have carried out analytical and experimental 
studies of air-to-air jet pumps, primarily for applications requiring high pressure 
rise or thrust augmentation. The entrainment ratios developed by these jet pumps 
are low, generally less than 10. Thus, this work is not directly applicable to the 
high entrainment requirements of V/STOL aircraft systems. Nevertheless, this 
work provides useful guidance for the development of performance prediction tech- 
niques and design rules for high-entrainment jet pumps. A brief review of some of 
the principal air-to-air jet pump papers follows. 

The performance of constant-area jet pumps was analyzed by McClintock 
and Hood for a range of design and operating conditions in reference 2. The analy- 
sis was prepared by assuming incompressible flow but the influence of compressi- 
bility was discussed in qualitative terms. Empirical coefficients derived by testing 
were included in the theory. The influence of mixing length and the use of various 
multiple-nozzle primary flow geometries were studied experimentally. The jet 
pumps treated had entrainment ratios of 10 or less; the design goal for the study was 
achievement of maximum thrust augmentation. 

A one-dimensional method of analysis of jet pumps was developed by 
Keenan, Neumann , and Lustwerk in reference 4. The analysis was applied to both 
constant-area and constant-pressure mixing processes. Test results were obtained 
for jet pumps with secondary-to-primary area ratios up to 100, primary-to-secondary 
pressure ratios up to 200, and a primary-to-secondary temperature ratio of 1. 0. 

The various regimes of operation of jet pumps with supersonic primary flows and 
both supersonic and subsonic primary flows were described. The analytical re- 
sults given were for jet pumps developing substantial stagnation pressure rises 
(e.g. , pressure ratios from 2 to 10) at low entrainment ratios (under 10). 

Fabri carried out a number of experiments on jet pumps which supple- 
ment the results reported by Keenan, et al. Fabri's results, with supporting analy- 
sis, are given in references 5 and 6. These tests also were confined to low entrain- 
ment ratio jet pumps. Excellent agreement was obtained between analytical predic- 
tions and measured jet pump performance. 
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An extensive analytical and experimental program was conducted at the 
University of Minnesota Rosemount Aeronautical Laboratories on jet pumps with 
secondary-to-primary area ratios up to 36, primary-to- secondary pressure ratios 
up to 32, and primary-to-secondary temperature ratios up to 3.3. The results are 
presented in references 7, 8, and 9. Single-nozzle primary flows with constant- 
area mixing tubes were used in these studies. High entrainment ratios were not an 
objective of the jet pump design; typical entrainment ratios reported were less than 
3. A number of analytical results showing the influence of duct matching upon jet 
pump and system performance were included in the reports. 

An analytical procedure for constant-area jet pumps with subsonic pri- 
mary flow was developed in reference 10. A computer program was prepared for 
use in optimizing jet pump design for particular application requirements. An analy- 
sis applicable to supersonic primary flows is given in reference 11. The analysis 
was compared with test results, but was not otherwise applied for jet pump design or 
optimization. 

The performance of a high-entrainment jet pump was measured in the 
Wagner "Jet Induced Lift" boundary layer control system (figure 1 ). These tests 
were performed by the present investigators under NASA contract No. NAS 2-25 18. 
The results are reported in reference 12. The jet pump component of the system 
employed a variable-area mixing tube (designed in an effort to obtain constant pres- 
sure mixing) and a 9-nozzle cluster for the primary flow. The jet pump was tested 
in the system at secondary-to-primary area ratios ranging from 150 to 800, primary- 
to-secondary pressure ratios up to 26, and primary-to-secondary temperature ratios 
up to 5.5. The desired constant pressure distribution in the mixing tube was not 
achieved. The entrainment ratios predicted for the complete s;ys tem were notattained. 
The results of the NAS 2-25 18 program showed that the methods used to design the 
high-entrainment jet pump and to match it to the duct flow characteristics were in- 
adequate. These results provided the impetus for the present study. 
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1. 3 Purpose 


The objectives of this investigation were as follows: 

• to develop analytical procedures for predicting the 
performance of high-entrainment-ratio jet pumps 

• to demonstrate the application of these procedures 
to match a jet pump design to its connecting duct 
system 

• to verify the analytical procedures by testing a jet 
pump over a broad range of operating conditions and 
jet pump geometries. 

The analysis and experimental workwere confined to constant-area mixing 
tube geometries because both analysis and construction are simplified by this choice. 
The only other mixing process which can be analyzed without complication is the con- 
stant pressure case. However, no reliable methods are available for designing a 
mixing tube which will actually achieve constant pressure mixing. Furthermore, 
this condition can be achieved at only one operating point for a jet pump of fixed 
geometry. 


Some of the test results obtained under contract NAS 2-2518 indicated 
that the design of the nozzle cluster and its position in the mixing tube may have 
created either high pressure losses in the secondary flow or poor mixing conditions 
at the mixing tube inlet. The mixing process did not seem to be completed within 
the length of the mixing tube used. These effects were thought to be partially respon- 
sible for the difference between predicted and measured performance for the com- 
plete system. Therefore, an additional objective of the present investigation was to 
test two alternative "low-drag" nozzle cluster designs and an extended mixing tube to 
determine whether these design changes would lead to significant performance im- 
provements. 
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Section 2 
SYMBOLS 


2 . 1 

A 

A 

C 

C 

w 

d 

f 

So 

H 

J 

k 

K 

L 

m 

m 

max 

M 

P 

P 

p* 

q 

q 

r 

r 

o 

R 


Symbols Used in the Analysis and Test Results 

,.2 

area, ft 

area ratio = A,, C /A , dimensionless 
th w m 

gas specific heat at constant pressure, Btu/lbm-°F 
nozzle flow coefficient, dimensionless 
duct diameters , inches 
friction coefficient, dimensionless 

2 

dimensional constant, 32.2 lbm-ft/lbf-sec 
stagnation enthalpy, Btu/lbm 
conversion factor, 778.2 ft-lbf/Btu 
gas specific heat ratio, dimensionless 
duct loss coefficient , dimensionless 
duct length, inches 

entrainment ratio = W /W , dimensionless 

s p’ 

maximum entrainment ratio, dimensionless 

Mach number, dimensionless 

pressure, psfa or psia 

pressure ratio = P /P , dimensionless 
po so 

pressure ratio = P /P , dimensionless 
r m so 

dynamic head = pV^ /2 g Q , psfa 
dimensionless dynamic head = pV /2g Q P go , 
radius, inches 

outer radius of mixing tube cross section, inches 
gas constant for air = 53.35 ft-lbf/lbm-°R 
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T 

T 

ir 

v 

v 

v* 

w 


temperature, °R 


temperature ratio = T /T , dimensionless 

po so 

system thrust augmentation = w m ^/WpVp 2 , dimensionless 
velocity, ft/sec 


velocity ratio = v s 2 / v p 2 ’ 
velocity ratio = V m /Vp 2 , 
mass flow rate, Ibm/min 


dimensionless 

dimensionless 


O' 

p 1 

a 

6 

AP 


AP 


ext 


AP 

s 

AP * 
s 


parameter defined by equation (34) 
parameter defined by equation (35) 
parameter defined by equation (45) 
parameter defined by equation (46) 
pressure change, psf 

ambient pressure rise imposed upon jet pump system = 

P, - P , psf 
b ao’ K 

static pressure rise = P - P , psf 

m so’ 

dimensionless static pressure rise = (P - P )/P , 

^ m so 3 * * * 7 so’ 


AP. 


AP t * 


stagnation pressure rise = p mQ “ p so > psf 

dimensionless stagnation pressure rise = (P -P )/P 

b mo so so 


3 

p density, lbm/ft 

r thrust augmentation = W m V m AVpVp 2 , dimensionless 

r w wall shearing stress, psi 

(f> impulse function 
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Subscripts 

a 

b 

d 

e 

1 
j i 

m 

MT 

0 

P 

s 

SD 

th 

1 

2 


atmospheric condition 
blowing duct exit section 
original conical diffuser exit section 
test rig diffuser exit section 
suction duct inlet section 
overall duct loss coefficient 
section at end of mixing region 
mixing tube 
stagnation value 

primary stream variable at nozzle exit 

secondary stream variable at nozzle exit 

section duct 

primary nozzle throat 

section at primary nozzle exit 

section at end of accommodation region 


Superscripts 

( ) ' value of parameter at end of frictional mixing tube 

extension (Section 3.1.5) 

( ) "mass-momentum" averaged value from test 

results 
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2.2 Symbols Used in the Computer Program 


Formulation 

Name 

Computer 

Name 

Definition and units 

Constants 



So 

GO 

o 

dimensional constant, 32.2ft-lbm lbf-sec 

J 

CONV 

conversion factor, 778 ft-lbf/Btu 

R 

R 

gas constant, ft-lbf/lbm-°R 

Variables 



A 

ABAR1 

A,, C /A , dimensionless 
th w m’ 

A th 

ATH 

2 

nozzle throat area, ft 

A 

P 

AP 

2 

nozzle exit area, ft 

A 

m 

AM 

mixing tube area, ft 2 

Not Used 

ABAR2 

A /A , dimensionless 
p m 

C 

w 

CW 

nozzle flow coefficient, dimensionless 

C p2 

CP2 

primary specific heat (C ), location 2, 
Btu/lbm-°R p 

C s2 

CS2 

secondary specific heat (C ), location 2, 
Btu/lbm-°R P 

C m 

CM 

specific heat (C ) at location m, Btu/lbm-°R 
P 

H p2 

HP2 

primary stagnation enthalpy, station 2, 
Btu/lbm 

H s2 

HS2 

secondary stagnation enthalpy. station 2, 
Btu/lbm 

K , 
sd 

FDUCT 

suction duct friction coefficient 

K , 
mt 

FTUBE 

mixing tube friction coefficient 

M „ 
p2 

PMOK 

primary Mach Number, location 2, 
dimensionless 
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Formulation 

Computer 

Definition and units 

Name 

Name 


M s2 

SMOK 

secondary Mach Number, location 2, 


dimensionless 

M 

m 

EMOK 

Mach Number at location m, dimensionless 

m 

ENTR 

entrainment ratio, dimensionless 

P 

PBAR 

P /P , dimensionless 
po so’ 

P 

po 

PPO 

primary stagnation process, psi 

P . 

SOI 

PSOI 

secondary stagnation pressure at duct 
inlet, psi 

P 

so 

PSO 

secondary stagnation pressure at nozzle 
exit, psi 

p 

m 

PM 

static pressure at location m, psi 

P 

mo 

PM TOT 

stagnation pressure at location m, psi 

T 

TEAR 

T /T , dimensionless 
po so 

T 

po 

TPO 

primary stagnation temperature, °R 

T 

so 

TSO 

secondary stagnation temperature, °R 

T 

p2 

TP2 

primary temperature, location 2, °R 

T 

s2 

TS2 

secondary temperature, location 2, °R 

T 

m 

TM 

temperature at location m, °R 

V p2 

VP2 

primary velocity, location 2, ft/sec 

V s2 

VS2 

secondary velocity, location 2, ft/sec 

V m 

VM 

velocity at location m, ft/sec 

W P 

WP 

primary mass flow rate, lbm/min 

W 

s 

ws 

secondary mass flow rate, lbm/min 

W 

m 

WM 

total mass flow rate, W + W , lbm/min 

p s 

P - P 

mo so 

DELP 

stagnation pressure rise, psi 
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Formulation 

Name 

Computer 

Name 

Definition and Units 

P - P 

mo so 

P 

so 

DDELP 

dimensionless stagnation pressure rise 

P - P 
m so 

DSTAT 

static pressure change, psi 

P - P 
m so 

Pso 

DDSTAT 

dimensionless pressure change 

Pm V 

2g o 

ENER 

2 

kinetic energy at location m, lbf/ft 

p V 2 
m m 

2 g P 
°o so 

DENER 

dimensionless kinetic energy 

T 

AUG 

W V /W V n , momentum ratio, 
m m p p2’ 

dimensionless 


Arrays 



CHARI 

store solutions for dimensional table 


CHAR2 

store solutions for dimensionless table 


COUNT 

store number of solutions for plotting 


ENT 

store entrainment ratio for plotting 


PRISE 

store DDELP values for plotting 


PTITLE 

store title for plotting 


SUM 

store number of solutions for table 


THROS 

store AUG values for plotting 


TTITLE 

store title for plotting 


Section 3 


ANALYSIS 


3. 1 Formulation of the Mathematical Model 

In this section an analytical model is developed to predict the flow be- 
havior in a compressible flow jet pump with a constant-area mixing tube. The anal- 
ysis is intended to provide a complete description of the important flow parameters 
at specific locations within the jet pump and to describe the overall operation of the 
jet pump as an entrainment and thrust augmentation device. The analysis was pre- 
pared for air-to-air jet pumps. The parameters used in the analysis are listed in 
section 2. The geometrical parameters are shown in figure 2. 

The fundamental purpose of the analytical model is to develop the per- 
formance characteristics of the jet pump directly. These performance character- 
istics can be represented by plots of jet pump pressure rise and momentum ratio as 
functions of entrainment ratio for a number of values of primary jet pressures and 
temperatures and for various area ratios. The performance characteristics in this 
form are analogous to head vs. capacity curves or performance maps which are 
commonly used for pumps and compressors. 

The equations describing jet pump performance and flow behavior include 
the entrainment ratio as an independent parameter. The assumption of a particular 
value for the entrainment ratio (together with the inlet flow pressures and tempera- 
tures and the primary-to- secondary flow area ratio) allows calculation of all of the 
performance and flow parameters for that operating point. Then another value is 
assumed for the entrainment ratio and the calculation procedure is repeated. Suc- 
cessive points on the jet pump performance curves are determined in this way until 
the complete curve is traced out. In the present calculations, the entrainment ratios 
were limited arbitrarily to the range from 10 to 40. In some cases, the Mach num- 
ber of the flow in the accommodation region reached 1. 0 for an entrainment ratio 
less than 40. Higher entrainment ratios cannot be achieved in such cases because 
the constant-area mixing tube chokes when the secondary flow Mach number reaches 1. 0. 
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3. 1. 1 Assumptions Used in the Analysis 


The following assumptions are made to simplify the analysis without 
seriously compromising its accuracy: 

1. The values of specific heat at constant pressure (C) 
and the specific heat ratio (k) are expressed as func- 
tions of temperature; otherwise the gas is considered 
to be a perfect gas. 

2. Wall shear forces are assumed to be negligible when 
compared to the pressure forces and the momentum of 
the primary and secondary streams. (This assumption 
is reviewed in section 3. 1.4.) 

3. No heat is transferred across the wall of the jet pump. 

4. The mixing tube is assumed to have a constant cross- 
sectional area along its entire length. 

5. When the primary nozzle is operated at an off-design pres- 
sure ratio, the primary jet is assumed to expand or contract 
isentropically until the primary and secondary streams have 
equal static pressures. This adjustment process is assumed 
to take place in the accommodation region between sections 1 
and 2 (see figure 2 ) and is assumed to be completed before 
any mixing takes place between the two streams. 

6. The stagnation temperature of the primary flow is assumed 
to be sufficiently high that moisture condensation shocks do 
not occur as the flow expands. 

3.1.2 Analysis of the Accommodation Region 

The geometrical parameters and flow conditions in the jet pump are de- 
fined as shown in figure 2. The primary stream enters the accommodation region as a= 
very high velocity jet; its Mach number maybe as high as 3. 5. The large momentum of the 
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primary jet induces a secondary flow. In the region defined as the accommodation 
region, it is assumed that the primary and secondary jets do not mix, but the prim- 
ary jet expands or contracts until its static pressure matches that of the secondary 
stream. At the point where the static pressures are equal, denoted as section 2, 
the accommodation process is assumed to be complete and the flows are parallel. 
This accommodation process is generally accompanied by a series of oblique expan- 
sion and contraction shock waves as the primary flow area adjusts to match the local 
static pressure outside the jet. However, if the jet pump is operated close to its de- 
sign conditions, the degree of accommodation is small and the losses caused by the 
shock waves will be small. A simplified oblique shock analysis indicates that a noz- 
zle designed for 350 psia supply pressure can be operated down to 200 psia with 
a total pressure loss due to shock waves of only 3%. For the values of supply pres- 
sure to be considered here, the error introduced by treating both streams as isen- 
tropic flows is negligible. In fact, one of the aims of this research was to show that 
an assumption of isentropic flow during the accommodation process will produce 
good results even when the system is operated at conditions quite far from the de- 
sign point. 


As long as the ratio of (Pp o /P so ) is sufficiently high to guarantee a super- 
sonic primary flow, as in the cases being considered here, the primary mass flow 
rate may be calculated directly using equation (1). 


k + 1 


W 


kg c 

R 


k+ 1 


k-i 144 xP (60xA th C %v ) 


E2_ 

V^i 


(1) 


po 


By specifying an entrainment ratio, m = W /W , the secondary stream 

s p 

mass flow rate can also be found directly. 


W = m W 
s p 


(2) 


In a perfect gas, the local values of total and static pressure are related 
to the local Mach number by the following equation. 
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P o /p = (1+ ~ M 2 ) 


k- 1 


(3) 


At the end of the accommodation region, the static pressures of the primary and 
secondary streams are equal. Therefore, at the end of the accommodation region 
the following relation must be satisfied. 


k-i 9 k-l 2 lcl 

KJ <i + m „ 9 2 ) = p so / (i + ^m s2 ) 


po 


p2 


(4) 


The mass flow rate per unit area for an isentropic flow is given by the 
relation below: 


W/A = 



M 




1 + 


k-1 ^2 1 
— M ) 


k + 1 
2 (k- 1) 


(5) 


m 


The geometry of the constant area mixing tube requires that A p2 + A g2 = A 
Using equation (5 ) to represent A p2 and A g2> and inserting the appropriate unit 
conversion factors, the geometry condition becomes as follows: 


k + 1 


W 




k-1 .. z. 2 * 1 - 1 ’ 


144x60 


M P 2~> 


nr 


V) 


k+ 1 
2 (k - 1) ' 


P M „ 
po p2 


P M „ 
so p2 


■ = A 


( 6 ) 


R 


When W , m, P , P , T , T and A_ are specified, equations (4) 
p’ ’ po so po’ so m ^ ^ n 

and ( 6 ) can be solved simultaneously to obtain M p2 and M g2 . Equation (5 ) can be 

used to find A n and A 

p2 s2 

Since the flow of both the primary and secondary streams is assumed to 
be isentropic in the accommodation region, values of static pressure, temperature, 
and velocity can be obtained for each stream at location (2) by employing the follow- 
ing equations. 
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p2 


P po /(1 + V M p2 2) 


k- 1 


( 7 ) 


T = T / ( 1 + M 2 ) 

p2 po ' 2 p2 ' 


( 8 ) 


T s2 = V (» +i r- m s2 2 > 


( 9 ) 


V 


p2 


V 2k g o 

1 k- 1 


R (T - T ) 
' po p2 ; 


( 10 ) 


V 


s2 


-V?Lh. 

” k- 1 


R (T - T _) 
' so s2’ 


( 11 ) 


3. 1. 3 Analysis of the Mixing Region 

The primary and secondary streams enter the mixing region with equal 
static pressures and parallel velocities. In this region, complete mixing takes place 
and a uniform flow with constant properties across the channel is obtained at section 
m. Treating the mixing region as a control volume with completely specified enter- 
ing flows, the following equations can be applied. 


Continuity Equation: 

Mass Flow Rate: 

Equation of State ; 

Momentum Equation: 

I44x60xg A (P _ - P ) = 
& o m ' s2 m/ 


w 

m 

W (1 + m) 
P v ' 

(12) 

W 

m 

60 x p V A 
m m m 

(13) 

P 

m 

p R T /144 
m m 

(14) 

W V - 
m m 

W V _ - W V _ 
p p2 s s2 

(15) 
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Energy Equation: 


2 


V 


m 


W H 0 + W H 0 = W C T + 0 
p p2 s s2 m m \ m 2g JC 

' fa o m 


(16) 


where H is stagnation enthalpy; 


V 


H 


p2 


= C 0 T 0 + 


p2 


P 2 p2 2g Q J 


H 


s2 


V 


- C_ n T~ + 


s2 


's2 s2 2 g J 
”0 


The equation for the specific heat at constant pressure for air as 
given in the Gas Turbine Engineering Handbook , Gas Turbine Publications, Inc. , 
1966, page 4, is presented below. 

C. = . 24916 - ,482xl(T 4 Tj + .681xl(T 7 T. 2 (17) 

Equations (13) and (14) are combined to give 


T = 60 x 144 x (P V A /R W ) 

JVJ ' m m rv» m ' 


m m m 


m 


( 18 ) 


Equation (15) is written in the form 


V 


m 


W 

= E. 

W 

m 


V „ + m V 0 
p2 s2 


g A 

+ 60 x 144 x - 9 -— m (P _-P ) 
W ' s2 m 


(19) 


and equations (16) and (18) are combined to yield 


RW 


m 60 x 144 x V A C 

m m m 


[ H P 2 +mH s2] 


RW V 
m m 


2x60xl44xg JC A 
& o mm 


( 20 ) 


Since the values of T „ and T „ are known at the inlet to the mixing 

p2 s2 

region, C „ and C 9 can be evaluated directly using equation (17). Equations (12), 
pz s z 
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(17), (18), (19), and (20) represent five equations for the five unknown exit parame- 
ters, W , T , C. P_ and V . Therefore, these five equations represent a com 
m m m m m 

plete set which describe the properties of the flow leaving the mixing section. In 
most instances, T g an< ^ T g ^ are sufficiently close in value to permit the use of a 
constant specific heat, C g = C g 2 = C^. This simplifies the solution by removing 
equation (17) from the equation set. 

3.1.4 Jet Pump Performance Parameters 

The jet pump performance parameters which are of particular interest 
are the stagnation pressure rise and the momentum augmentation. These parame- 
ters can be evaluated by using the following expressions. 


and 


A P = P - P 
t mo so 


k 

k- 1 

AP. = P (1 + ~rr- M 2 ) - P 

t m 2 m ' so 


( 21 ) 


r = outlet momentum /primary momentum 


T = 


W V 
m m 

W V . 
P P2 


(m + 1) V. 


m 


V 


( 22 ) 


p2 


During the jet pump test program, the measurement of the stagnation 
pressure rise produced by the jet pump proved to be difficult to accomplish. It was 
much easier to measure the stagnation-to- static pressure rise in the jet pump. 

This pressure rise can be used to define another jet pump performance parameter, 
AP g , which serves as an alternative to the stagnation pressure rise parameter, A 


AP = P - P 
s m so 


(23) 


All of the values needed to compute AP,, A P ,and r are provided by the 

t S 

analyses of sections 3. 1. 2 and 3.1.3. 
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3.1.5 Inclusion of Mixing Tube Wall Friction 



The analysis of the mixing region presented in section 3.1.3 neglected 
the effect of wall friction; this simplifies the equations describing the mixing process. 
In this section, a procedure is developed to include the effects of wall friction in the 
jet pump analysis. 

The mixing region was assumed to extend from the point where the pri- 
mary and secondary stream pressures are equal to the point where they have merged 
into a uniform flow with constant properties across the channel. In reality, the wall 
friction effects occur in conjunction with the mixing process. Unfortunately, it is 
difficult with the current state of knowledge to predict wall friction losses accurately 
in the mixing region. Therefore, rather than adjust the mixing region analysis to 
include the wall friction effects, we considered that it would be preferable at this 
time to treat the mixing process and wall friction as independent effects by imagining 
the mixing tube to extend as shown in figure 3 beyond the point where the mixing 
process is complete. The flow phenomena occurring in the mixing portion, segment 
I, is a mixing process without wall friction as analyzed in section 3.1.3. The fric- 
tional portion, segment II, represents the effect of wall shear forces upon a uniform 
adiabatic flow. The hypothetical extension of the mixing tube is meant to represent 
the friction occurring within the actual mixing tube. 

The effect of wall shear forces in ducts is commonly represented by a 
coefficient of friction defined as follows: 

f = 7 ^— (24) 

PVV 2g o 


where t w is the shearing stress exerted upon the stream by the wall. The corres- 
ponding stagnation pressure loss is given by equation (25). 


AP 

o 




(25) 
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where 


L = duct length 

D = duct diameter 


In many flow analyses, the value of the coefficient of friction is taken 
from test results for pipe flow. The friction coefficient is a function of wall rough- 
ness and Reynolds number. However, in the case of the mixing tube, the wall fric- 
tion occurs in a very non-uniform flow which has a high level of turbulence. The 
value of the coefficient of friction for the mixing tube cannot be accurately determined 
from pipe flow data. Therefore, it is convenient to represent the mixing tube wall 
friction loss in terms of a head loss factor, Kjyjrj., which must be determined experi- 
mentally. 


AP 

o 


mixing tube 


K 


MT 


p V 
m m 

2 S n 


(26) 


The equations for an adiabatic flow in a constant- area tube with a stagna- 
tion pressure loss are given below. Primed variables denote parameters at the end 
of the hypothetical extension of the mixing tube. 


Momentum Equation: 
k 


J ( 1 + ~ 
m\ 2 


144 P 1 + — rr M ~J 
m ' 


V ^ 


'( 1 + nr M ™' 2 ) 


k 

k- 1 


144 P '1+ _ 

mV 2 


= K 


P V 
m m 


MT 2g r 


(27) 


Energy Equation: 


V 


C T + 


m 


V 


,2 


mm 2 g 


C T ’ + -jp- 
m m 2 g 


(2 8 ) 


Continuity Equation: 


p V = p ’ V ' (29) 

m m m m 
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State Equation: 


P 

m 


P RT 
m m 


P ' 
m 


p 'RT ' 
m m 


Definition of Mach Number: 



V 


2 


S r 


m 

k R T 


m 



V 


,2 


Sr 


m 

k R T ' 
m 


(30) 

(31) 

(32) 

(33) 


These equations can be combined to give the following results: 

k 


p Yi + m ' 2, \ 

ml 2 m ) 


k- 1 


= a 


(34) 


where 




(35) 


(36) 


P 





k- 1 
2 



(37) 


Eliminating P ' between equation (34 ) and equation (35) yields equation 
m 

(38) which can be used to determine M 
' ' m 

k + 1 

vh ( 1 + iL r- M m ,2 ) 2<k 1>= “/•#■ (38) 

m ' 
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After equation (38) is solved for M ', equation (34) can be used to determine P 

m m 

Combining equations (29), (30), (31), (32), and (33), an equation for T ' is obtained. 

T ' = T (P '/P ) 2 (M '/M ) 2 (39) 

m m ' m nr ' m nr ' ' 

Equation (33) is then used to determine V 1 . 

m 

The equations developed in this section can be used to compute the values 
of P m \ T \ M ', and V ' for a jet pump when the "ideal" analysis of sections 
3. 1. 2 and 3.1,3 is completed and the value of is known or assumed. Alterna- 

tively, these equations may be used to deduce the value of K..,_ when values of P 

MT m 

T ', and V ' are known from test results and values of P , T , and V are com- 
m m m m m 

puted by using the "ideal" analysis. 

3.1.6 Dimensionless Formulation 


In this section the equations describing the jet pump operation are for - 
mulated in terms of dimensionless variables. The non-dimensional formulation is 
valuable for two reasons: 

• According to the principles of dimensional analysis, a 
solution in terms of independent non-dimensional groups 
is a general one. The same solution may be applied for 
jet pumps having great differences in individual design or 
operating parameters so long as the independent non- 
dimensional groups are identical. For example, one 
such group is the primary-to-secondary flow area ratio, 

AfhC ^/A ; if all other non-dimensional groups are the 

same, a large-scale and small-scale jet pump having 
identical area ratios will have identical non-dimensional 
performance characteristics. 

• The non-dimensional formulation permits identification of 
the minimum number of independent non-dimensional groups 
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which are required to completely specify a jet pump design 
and its operating characteristics. 


In the derivation which follows, unit conversion factors are not included, 

and all values of the specific heats (C , C , C 1 are assumed to be constant and 

' p s nr 

equal. 

The equations which apply in the mixing region are given below. 
Continuity Equation: 

Mass Flow Rate: 

Equation of State : 

Momentum Equation: 


w 

m 

W 

P 

(1 +m) 

(12) 

W 

p 

V A 

(40) 

m 

m 

m m 

P 

m 

P m 

R T 

m 

(41) 


gA (P „ - P ) = W V -WV 0 -mWV 0 
o m v s2 nr mm p p2 p s2 


(42) 


Energy Equation: 


T +mT = (m + 1) C 
po so v ' 


V 


2 


m 


T + 

. m 2 e 0 c -> 


(43) 


The dimensionless variables to be used are defined as follows. 


V* = V /V P* = P /P , P = P /P , 
m p2 m so po so 


T = T /T , A = A., C /A 
po so th w m 


(44) 


Using these variables, equations (12), (40), (41), (42), and (43) can be combined 
as follows. 

p * - < p s2 /p S o> + r [i + “ < V S /V ‘ < m + i> v *] 


(45) 
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v* 2 + _JL_ 2 P* V* _ 6 m 

(k - 1) y (m + l) (m + l) y 


0 


(46) 


where 


y = W V 0 /g A P 

p p2/ °o m so 


(47) 


and 


2 g C T /V ‘ 
& o po p2 


(48) 


Equations (45) and (46) show that P* and V* are related through y , 6, 


P s2 /P so- V s2 /V p2- m ' k ’ ” d T ' 


The equations which govern the flow in the accommodation region are 
developed next. For an isentropic primary stream, 


k- 1 


k- 1 


T JT = (P J P ) 
pZ po p2 po ; 


± P s2 

- P P SO. 


(49) 



(A.. C ) 
' th w' 


(50) 


and 


p2 


'2 g D k R 
k- 1 


T (1 - T J T ) 
po v p2 po' 


(51) 


Equations (49), (50), and (51) can be combined to give 


k +1 
k- 1 


y = ? 


Usl 

k - 1 k + V 
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( !k, 

p 

so -J 


(52) 


24 



Similarly, equations (44) and (51) can be used to obtain a value for 5: 



Thus, y and 6 are shown to be functions of P, A, k and P /P 

s2 so 


The secondary stream velocity is given by equation (54): 


V 


s2 


V 2g ° kR 

k- 1 



Equations (49), (51), and (54) can be combined to yield 



From the definition of total or stagnation pressure, 



1 / 




k 

k- 1 


(54) 


(55) 


(56) 


Considering equations (51), (53), (55), and (56) together, it can be seen 
that the parameters y , 5, V g9 /Vp 2 > ■ P s2 / ^ P so are ft™ 0 ** 0118 °f M g2 , P ’ T ’ A ’ 
m and k only. Equation ( 5) of section 3. 1. 2 can be written for the secondary 
stream as follows. 


M s2 

k + 1 

f k-1 2 \ ^ ^ - 1) 

, 1 + V M S 2 


w V t"” 

s Y so 
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(57) 
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Equation (57) can be combined with equation (50) and written in terms of 
dimensionless variables as given below. 



(1 + 


k- 1 
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s2 


k+l 
2 (k - 1 ) 


m (- 


•k+l 


k+ 1 

2 (k - 1 ) 

) P A 


yr ^-yv 


(58) 


For high entrainment ratio jet pumps, the term A p /A m is of order 0. 01. 
Thus, the term (1-A /A^) can be approximated as unity. Equation (58) can now be 
written as follows. k+l 

M s2 , 2 * PA 


k+l 

k -1 2 W- 1 ) 

< i + V m s2 > 


= m (! 


k + 1 


) 


yr 


(59) 


Equation (59) shows that M g2 can be determined from m, P, A, T, and k. This in- 
dicates that 7 , 6 , V „/V „, and P J P are functions of P, T, A, k and m. Re- 
SZ pZ sz so 

turning to equations (45) and (46), it can be seen that the performance of the jet pump 
depends upon the parameters P, T, A, k and m. For given values of P, T, A and k, 
a complete dimensionless solution can be obtained for each specified value of entrain- 
ment ratio. 


The jet pump performance parameters can be experssed in the form of 
dimensionless groups using the fundamental dimensionless parameters. The momen- 
tum augmentation, r , is already a dimensionless group: 


T = 


W V 
m m 

W V 0 
P P2 


= (m + 1) V* 


(60) 


The dimensionless stagnation pressure rise parameter is defined as follows: 

k 

AP ^ 1 

* k -1 2 

AP t = -p 1 = P * < 1+ V M m ) - 1 

so 


(61) 
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To evaluate AP, *, it is necessary to know M . Using the dimensionless variables, 
2 ^ 

M m is given by equation (62). 



m 


kg RT 
to o m 


V* 
p * 


, m + 1 
( k 


) y 


(62) 


Thus, AP^.* is a function of the fundamental dimensionless variables which deter- 
mine V*, P*, and 7 . The static pressure rise parameter, AP , can be expressed 

s 

in dimensionless form as follows: 

P -P 

AP * = = P* - 1 (63) 

s P 

so 


Equations (60) through (63) show that the dimensionless jet pump per- 
formance parameters are functions of the fundamental independent dimensionless 
variables P, T, A, k, and m. Five such independent variables and only five have 
to be specified in order to determine the jet pump performance characteristics in 
dimensionless terms. (This conclusion is restricted to jet pumps which satisfy the 
assumptions listed in section 3. 1. 1 and the additional assumption that the specific 

heats of all of the streams are equal, i. e. ,C=C n =C=C .) 

p 2 s 2 m 


It is possible to use a different set of five independent dimensionless 
variables. For example, a velocity ratio V = V g2 /V 2 can use( ^ * n Pl ace T to 
complete an alternative set of five independent variables, P, V, A, k and m. An- 
other possible set is P, T, A, k, and V. The velocity ratio, V, was one of the basic 
design parameters used to select the jet pump geometry for the boundary layer con- 
trol system tested under contract No. NAS 2-2518. 


These remarks can be summarized by the expressions below: 


Jet Pump Performance 
Characteristics 
(dependent variables) 


Design and Operating 

Conditions 

(independent variables) 


AP, *, AP * and r are functions of 
1 / s 


P, T, A, k, m 
or 

P, V, A, k, m 
_ or _ 
P, T, A, k, V 
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3.1.7 Maximum Entrainment Ratio 


Equation (59) can be used to determine the maximum possible entrain- 
ment ratio. This maximum occurs when the secondary stream Mach number reaches 
a value of unity. Setting M^ = 1, equation (59) can be written as follows. 


Vt 7 


P A 


(64) 


Thus, the maximum entrainment ratio corresponding to the choking of the secondary 
stream can be determined directly from the dimensionless initial conditions. 


3. 2 The Computer Program 


A computer program was prepared to predict the performance charac- 
teristics of constant area jet pumps using the analytical concepts formulated in the 
preceding sections. The program was written to develop both dimensional and di- 
mensionless solutions. Values of P, T, A, P , T , and A are read in as initial 

so so m 

conditions. Values of k, R, and m are included within the program. For each value 

of m, values of t , AP., AP , AP * , and AP * are calculated. The values of r , AP * 

L S t S t 

and AP g * depend only on the dimensionless data, while AP^. and AP g depend also upon 

P , T , and A . 
so so m 


The program is written in Fortran IV language. The machine used was an 
IBM System 360/65 with an SC4020 plotter. Automatic plotting of the performance char- 
acteristics was obtained by using the subroutine EZPLOT developed by the Missile Sys- 
tems Division of Avco Corporation, Burlington, Massachusetts. 

A block diagram of the computer program is shown in figure 4. For each 
set of initial conditions, solutions are obtained for values of entrainment ratio between 
10. 0 and 40. 0 in steps of 3. 0. The results are printed as each solution corresponding 
to a particular value of entrainment ratio is determined. The results are also stored 
in arrays for plotting and for presentation in tabular form. A printout of the entire pro- 
gram is presented in Appendix A. A discussion of the program by blocks is given in 
Appendix B. Appendix C provides typical computer solutions which indicate the form of 
the output data. 
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Certain of the blocks shown in figure 4 and described in Appendix B are 
denoted as optional, indicating that they can be removed from the program without 
interfering with the operation of the remaining blocks. Instructions for removing 
these blocks are given in Appendix B, section B-2. 

When frictional effects in the suction duct and mixing tube are to betaken 
into account in the performance predictions, values of the loss coefficients K g ^ (de- 
fined in Appendix B) and must be provided as input data for the computer solu- 

tion. These loss coefficients are functions of the flow Reynolds numbers. If values 
of K gc j and have been established for ducts of one size and the computer per- 

formance predictions are to be used for ducts of much larger or smaller sizes, it 
may be advisable to adjust the values of the loss coefficients used by the computer 
to account for the Reynolds number change. 

3. 3 Solutions for a Ran g e of Jet Pump Designs 


The computer program was applied to develop jet pump performance 
plots for a broad range of geometries and operating conditions. The range of solu- 
tions was selected to encompass all of the test conditions used in this investigation (sec- 
tion 4) and also the range of conditions of interest to NASA for boundary layer con- 
trol systems and momentum augmentation. The performance plots were developed 
for use in preliminary design of jet pump systems for matching the jet pump to a 
duct system and for predicting the resulting system performance characteristics. 
Techniques for applying the solutions to system design are described in section 3. 4. 

The range of conditions used to obtain the performance plots were ini- 
tially defined in dimensional form as follows: 

Tp o primary flow stagnation temperature 450° F to 3500° F 

Pp Q primary flow stagnation pressure 100 psia to 400 psia 

T go secondary flow stagnation temperature 20° F to 120° F 


P secondary flow stagnation temperature 1500 psfa to 2116 psfa 
s o 


The range of values selected for the corresponding dimensionless parameters were 
as follows: 
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T = 1. 5 to 8. 0 

P = 5 to 40 

The nozzle and mixing tube geometries available for the test program had area 
ratio values (A) ranging from 0.00125 to 0.0067. The range of values selected for 
the performance plots is given below: 

A = 0. 001 to 0. 007 

The ranges of values of P, T, and A given above were used to prepare 
9 sets of performance plots showing AP^* vs. m and r vs. m for various values 
of P with T and A fixed; these plots are indexed in table 1 and given in figures 6 
through 23. Table 1 includes the values of maximum entrainment ratio attainable 
for each combination of P, T, and A values. This maximum entrainment ratio is 
set by choking of the secondary stream as given by equation (64) of section 3. 1. 7. 

Typical computer output sheets for one of the solutions are reproduced 
in Appendix C. The printed output includes values of jet pump parameters not 
shown in the plots but required for the jet pump- duct matching techniques 
described in section 3. 4. These parameters are given in dimensional form based 
upon standard secondary stream inlet conditions, P = 2102 psfa, T = 70 D F. 

q SO so 

and A = 0. 08726 ft' . 
m 


In the last six cases, the higher values of P cannot be attained because of 
choking of the flow in the mixing tube. For the cases with T = 1.5 and 3.5, choking 
occurs in the secondary flow (M g £ = 1) as discussed in section 3. 1. 7. When T was 
set at 8.0, choking was predicted to occur first at the mixing tube exit, i. e. , M m = 1. 

Cross-plots showing AP^* vs. m and r vs. m for various values of A 
with T and P fixed are presented in figures 24 and 25. Lines of constant mixing tube 
exit Mach number are also shown. These cross-plots provide additional insight on 
the effect of the area ratio upon jet pump performance. 
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3.4 Jet Pump-Duct Matching Considerations 

The previous sections have developed analytical techniques and a com- 
puter program which allow prediction of the performance characteristics of high- 
entrainment- ratio jet pumps. The performance characteristics take the form of 
plots of jet pump pressure rise and momentum ratio as functions of entrainment 
ratio for a number of values of primary jet pressures and temperatures and for var- 
ious area ratios. These performance characteristics are analogous to head vs. ca- 
pacity curves or performance maps which are commonly used for pumps and com- 
pressors. The actual point (i.e. , entrainment ratio) at which a jet pump will oper- 
ate when connected to a particular system of inlet and discharge ducts is dictated by 
the geometry of the duct system. 

The resistance curve of a duct system is roughly parabolic as shown in 
figure 26. A typical jet pump characteristic is also shown on the figure . The actu- 
al operating point of the jet pump-duct combination is defined by the intersection of 
the two curves. The duct characteristic curve is set by the duct geometry and is 
essentially independent of the jet pump operating conditions. Therefore, if the duct 
geometry is not changed, the operating point of the system for any jet pump primary 
flow condition must be located on the parabola. Figure 27 shows how the operating 
points for a system can be determined if the jet pump performance at various pres- 
sure levels is known. 

This section establishes a procedure for use to determine the operating 
points of a jet pump in a duct system when the loss characteristics of the duct sys- 
tem are known. This procedure can be employed as shown by example to match the 
jet pump design to the duct system so as to achieve maximum entrainment ratio or 
maximum thrust augmentation for given primary flow conditions. 

3.4.1 Representation of Duct Loss Characteristics 

When the analytical model developed in sections 3. 1 and 3. 2 is supplied with 
values of P, T and A, the performance parameters such as stagnation pressure rise 
and momentum ratio (thrust augmentation) can be calculated as a function of entrain- 
ment ratio. In order to determine the specific value of entrainment ratio which will 
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be obtained during operation of a given jet pump of fixed geometry, the associated 
duct system flow characteristics must be taken into account. Using the notation 
shown in figure 2 , the stagnation pressure rise may be written as follows: 


A P, = 


mo 


- P = (P 


SO 


ao 


-P ) + (P - 
so ' mo 


P bo> + 


bo 


ao 


( 65 ) 


The two bracketed terms together represent the stagnation pressure loss of the duct 
system including the inlet duct (first term) and blowing duct (second term). 

At low Mach numbers, the stagnation pressure loss due to friction in a 
duct is proportional to the kinetic energy of the flow. For high entrainment ratio 
jet pumps, the mass flow rates in the suction and blowing ducts are nearly equal. 
Therefore, the total pressure loss of the entire ducting system may be related as a 
first approximation to the kinetic energy of the blowing duct inlet flow by equation (66) 


(P - P. ) + (P - P ) = K. P m V m 
' mo bo 7 ao so' £ 


2 S, 


( 66 ) 


This type of expression has been shown to be accurate for representing frictional 
losses in duct systems of various shapes. 

3.4.2 Evaluation of the Loss Coefficient, K 

The loss coefficient depends on the geometry of the particular ducts 
being used and the Mach number level (ref. 13). At sufficiently low Mach numbers, 
(i.e. , under 0.3), compressibility effects can be neglected. For flows at higher Mach 
numbers, the value of can be corrected for compressibility effects. 

Loss coefficients have been presented for a number of duct configura- 
tions in references 13 through 26, The configurations reported include ducts of rec- 
tangular and circular cross section with varying amounts of diffusion or acceleration. 
Bends and elbows having a number of different angles of turn are included in these 
references. The loss coefficients reported were measured for subsonic flow cover- 
ing a range of Mach numbers up to 1. 0. 
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To provide an example of typical loss data, the influence of Mach num- 
ber level upon the loss coefficient for straight, conical diffusers is shown in figure 
28. The variation of with inlet Mach number is significant for a diffuser of speci- 
fied geometry. 

Because the value of is closely related to the duct configuration and 
the Mach number level, and values of are readily available for only a few simple 
duct shapes, the designer of a jet pump and duct system generally will not be able to 
look up an accurate value for for a new duct design. If optimum matching of the 
jet pump and ducting is required, the loss coefficient of a new duct geometry will 
have to be determined experimentally. Testing can be done by using either a full- 
scale or reduced scale model of the duct. The tests must cover the Mach number 
range which will be encountered by the actual duct when operating with the jet pump. 
Flow tests of ducts sometimes have additional value; regions of flow separation or 
undesirable velocity profiles may be revealed. When the duct geometry is modified 
to eliminate these problems, the loss coefficient is usually reduced. 

3.4.3 Development of System Performance Equations 

At the outlet of the blowing duct, the static pressure in the flow must be 
equal to the local "atmospheric" pressure. The use of a blowing duct having the 
same cross-sectional area as the constant-area mixing tube will limit the entrain- 
ment ratio which can be achieved in the jet pump system. Higher entrainment ratios 
can be obtained with the same jet pump if a diffuser is added to theblowing duct. The 
diffuser allows higher velocities and flow rates in the mixing tube. The mixing tube 
pressures can be sub-atmospheric; the diffuser decelerates the flow to increase its 
static pressure up to the atmospheric pressure level at the blowing duct exit. 

A calculation method can be developed for use to determine the actual 
operating point (i. e. , entrainment ratio) for a jet pump system as a function of the 
area ratio selected for the blowing duct diffuser. The calculation method makes use 
of the generalized jet pump performance characteristics developed by the computer 
program described in section 3. 2 and Appendix B. 
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The loss coefficient defined in equation (66)is inserted in equation ('65) 
with the following result: 


AP t = 


< P bo- 


P ) 
a o' 


+ K„ q 


“£ m 


(67) 


where 


q m = 


P V 
m m 
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Equation (67)may be rewritten as follows: 


A P, 
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m ' bo b ; 


( 68 ) 


The term (P^-P^) represents the external ambient pressure difference imposed 
upon the jet pump system. This term will be called A P 


AP 


ext 




P ) 

ao' 


(69) 


The value of AP was zero for the experimental jet pump since both 
the discharge static pressure (P^) and inlet stagnation pressure (P ) for the jet 
pump system were equal to atmospheric pressure. This term is not necessarily 
zero for jet pump systems which operate in the presence of an external velocity 
field. For example, a jet pump used for boundary layer control at the trailing edge 
of a wing will have its inlet pressure (P ) established by the flow behavior in the 
suction slot entry passages and by the local pressure acting on the wing. The dis- 
charge pressure (P^) will be set by the local pressure field on the wing and by the 
flow behavior from the slot to the deflected flap. 


The term (P^ q - P^) in equation (68) represents the dynamic head of the 
flow at the blowing duct exit. This term is related to the blowing duct exit Mach 
number as shown in equation (70) : ^ 
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(70) 
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In order to evaluate this term, the value of must be calculated using the selected 
blowing slot area A^, the exit pressure P^, and the flow conditions at the discharge 
of the jet pump mixing tube: 


M b 


V 


1 + 


k- 1 
2 




V T mo ' r m 

P b * kg o 


(71) 


This equation is based upon the assumption that the flow in the blowing duct is adia- 
batic. The values of T mo and W are output values from the jet pump performance 

calculations described in the previous section. The evaluation of M, can be made con- 

b 

veniently by using figure 29. 


Equations (70) and(71)can be replaced by equation(72) when the blowing 
duct exit Mach number, M^, is less than 0. 3. 


P - P 
bo b 


^b V b 

2g„ 


(72) 


The error in using equation (72) in place of equations (70) and f71)is less than 2% of the 
true (P^ - P^) difference when is less than 0. 3. 

A particularly simple jet pump-duct system matching equation can be 
derived when equation (72) is used. Equation (73) is the continuity relation for the 
blowing duct: 


p V A = p, V, A, 
m m m b b b 


(73) 


The influence of the blowing duct area ratio is introduced when equation(73) is com- 
bined with equation (72): 


bo 





(74) 


The jet pump-duct matching equation, equation (75) , is derived by combining equa- 
tions (68), (69) and (74) ; 
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AP fc = 


AP . + q 
ext m 


K„ + 


m 


A, 


m 


(75) 


for 




0. 3 


For preliminary design purposes, the value of P m /Pj. } can be taken as 
1.0. A more accurate value can be determined as follows: 

For a perfect gas, 


m 


m 


(76) 


m 


With less than 0. 3, equation (77) holds with an error of less than 2%: 


T k = T, 
b bo 


(77) 


Since the flow in the blowing duct is adiabatic, its stagnation temperature remains 
constant; 


T. = T 
bo mo 


The last relation required is equation (79): 


(78) 


T 

mo 


1 + 


k- 1 
2 



(79) 


When equations (76) through (79) are combined, an equation for calculating P m /P b 
is derived: 





k- 1 
2 



( 80 ) 
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Equations (74), (75), and (80) can be used with small error only if 
is less than 0. 3. If is greater than 0. 3, the jet pump-duct matching equation, 
equation (81), is derived by combining equations (08), (69), and (70): 


AP 


t 


AP ext + 


K. q + 
£ m 



(81) 


for all values of M, 
b 


where is computed by using equation(71) and figure 29. 

There are two figures of merit which are of interest in the evaluation of 
jet pump systems for boundary layer control or thrust augmentation applications. 
These figures of merit are the entrainment ratio and the thrust augmentation. The 
equations developed above can be used to determine the entrainment ratio at which 
a jet pump system will operate. Several additional equations are required in order 
to calculate the thrust augmentation obtained from a jet pump system. 

The thrust augmentation obtained with the complete system is defined in 
equation (82): 


IT 

TT 


system thrust augmentation 


W V, 
m b 

W V . 
P P2 


(82) 


The thrust augmentation produced by the jet pump alone was designated as t in section 
3.1.4 (equation 22). The computerized jet pump performance analysis of section 3.2 
provides as output data values of r as a function of entrainment ratio. Thus, once 
the entrainment ratio is known for a jet pump system, the value of r is known and 
the system thrust augmentation can be calculated as follows: 


so 


T = 


W V 
m m 

W V „ 
P P2 


V, 


TT 


V_ 


m 


(83) 
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The value of V^/V^ can be related to the blowing duct area ratio by using equation 
(73): 



(84) 


The value of p^/p^ is given by equation (80) when ^ 0. 3. When exceeds 
0. 3, the value of is given by equation (85): 



(1 + 


k- 1 

2 


(1 + 


k- 1 
2 


M 2 ) 
m ' 


V) 


(85) 


where the blowing duct flow has been assumed to be adiabatic and the value of M^is 
determined by using equation(71) and figure 29. 

The equations given above can be used to compute the thrust augmen- 
tation parameter once the operating point of the jet pump is known. The next sec- 
tion establishes a procedure for determining the operating point. 

3.4.4 Calculation Procedure for Determining the Operating Point of a Jet 
Pump in a Duct System 

The operating point of a given jet pump and duct combination can be de- 
termined as follows: 


Required Initial Data : 


Given: Jet pump design and performance characteristics: 


basic jet pump 
parameters 


so 


po 

T 

so 

T 

po 

A,, C 
th w 


m 
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jet pump performance 
curves; output data from 
computer program 


AP t 

vs 

m 

q m 

vs 

m 

M 

m 

vs 

m 

P 

m 

vs 

m 


values characterizing the 
particular duct system < 



Solution technique if ^ 0.3: 

The specified values of , A P ext , P^, and are inserted into equa- 
tions (80) and (75) . The jet pump performance curves are used to find associated 
values of A P^., q , M^, and P m which satisfy equation (75) . This is a trial-and- 
error process which is begun by assuming a value for entrainment ratio, m. The 
corresponding values of and P m are determined from the jet pump performance 
curves and entered into equation (80). The resulting value of * s entered, 

together with the value of q m from the jet pump curves, into the right-hand side of 
equation (75). If the resulting value of A P does not agree with the curve value, a 
new value of m is assumed and the process is repeated. The iteration process is 
simplified by graphical solution techniques which are described in the section en- 
titled "Sample Calculation". This calculation process finds the value of m at which 
the jet pump system will operate with the selected value of A^. 


Solution technique if M^> 0. 3: 


AP gx j, P^, and A^ are inserted into equa- 


The specified values of K 
tions (81) and (71), The jet pump discharge flow rate, W , and the stagnation tem- 
perature of the discharge flow, T mo > are plotted as functions of entrainment ratio, 
m. Values of can be determined as a function of m using equation (71)and figure 
29. The solution technique is a trial-and-error process which is begun by assuming 
a value for m. The jet pump performance curves are used to find values of q m and 
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A Pj. for each value of m assumed. Corresponding values of are determined as 
above and entered with the q m values into the right-hand-side of equation (81). If 
the resulting value of AP^ does not agree with the curve value, a new value of m is 
assumed and the process is repeated. Graphical solution techniques, described in 
the section below, can reduce the number of iterations required. 

Evaluation of thrust augmentation: 

The solution techniques described above yield the value of entrainment 
ratio at which a jet pump will operate in a selected duct system. The performance 
data provided by the jet pump computer program allows determination of the values 
of the following jet pump performance parameters at the operating point: r , P , 
Mm’ ^m’ anc * ^mo' These values allow calculation of the thrust augmentation, IT, 
by using the equations presented at the end of section S. 4. 3. 


Additional Comments: 

The values of K . and P are not constant for all values of entrainment 
£ so 

ratio. At high entrainment ratios, the Mach number levels within the ducts may 
become sufficiently high that the influence of compressibility upon must be taken 
into account. Similarly, the value of P , which is a non-dimensionalizing param- 
eter in the jet pump performance analysis, varies slightly as shown in equation (86) 
when the entrainment ratio (and secondary stream flow rate) changes. 

P = P -K. q where q varies with m (86) 

so ao i mn mi v ' 

The variations of K„ and P with entrainment ratio are generally second-order in 
£ so 

magnitude. These variations can be neglected in preliminary design calculations, 
then included for final design if and K. are known as functions of the Mach num- 
ber M = M.. 

m l 

3.4.5 Sample Calculation 

The use of the procedure described above to determine the operating 
point of a jet pump-duct system is illustrated by the sample calculation which follows 
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The jet pump design data is: 


p 

300 psia J 


po 

c p = 

20. 13 

p 

so 

14. 9 psia V 


T 

1200° F ) 


po 

} T = 

3.074 

T 

80° F \ 


SO 



A 

0.08726 ft 2 f 


m 

2 } A = 

545.4 

A..C = 
th w 

0. 000165 ft J 



(Test value for the Case 4 
nozzle) 


The computer solution for the jet pump performance yielded the values given in 
table 2 . The values of AP^, p M^, an d c l m are plotted against entrainment 
ratio in figure 30. 

The duct system design conditions were assumed to be as follows: 

AP ext = 0 P b= P ao> 

K l = °' 1 

P =14.9 psia 

ao ^ 

The blowing slot discharge Mach number, M^, was assumed to be less than 0. 3. 
Using these values in equation (75) equation (87) was derived. 



The calculations were begun for a blowing duct diffuser area ratio, 
(A, / A m ) , equal to 1.0. Three values of entrainment ratio, m = 13, 15, and 17, 
were selected arbitrarily. The corresponding values of A P^, M m , q m? and P m 
were read off from figure 30. These values were used to compute p /p^ from 
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equation (80) and then to 
hand side will be called 


compute the right-hand side of equation (8?) ; the right- 
A Pt ^ trial' resu ^ s are gi ven in the table below. 


For A, = 1.0: 

b 


m 

AP t 

(psf) 

M 

m 

V 

(psf) 

13 

99. 6 

. 205 

65. 3 

15 

98. 5 

. 235 

87 

17 

97.5 

. 270 

109 


P p m AP,) . . , 
m t J trial 

P 


psia 

P b 

(psf) 

15. 15 

1.02 

74. 9 

15. 00 

1.01 

96. 8 

14. 80 

1. 01 

121 


The A P^ trial va ^ ues can plotted against entrainment ratio as shown in figure 31 
The intersection of the A P^. j trial curve W3 th the A P^ jet pump performance curve 
represents the solution of equation (87) for the selected value of A,/A m - This in- 
tersection is the operating point of the jet pump in the specified duct system. 


3.4.6 Influence of Blowing Duct Area Ratio and Duct Losses Upon Entrainment 


Similar calculations were carried out for values of A^/A m equal to 2, 3, 
and 4. The results are shown in figure 31. The entrainment ratio increases as the 
blowing duct area ratio is increased; the trend is more clearly shown when the results 
are replotted as in figure 32. For the particular jet pump and system design condi- 
tions assumed for this sample calculation, the maximum entrainment ratio is 
achieved when the mixing tube is choked, i.e., when = 1.0. 

The influence of the duct loss coefficient was explored by setting = 

0. 2 instead of 0. 1 as previously assumed. The duct matching calculations were re- 
peated using equation (75); the results are shown in figure 32. Only a small increase 
in entrainment ratio can be obtained by increasing the area ratio from 4 to 5. This 
is a consequence of the fact that, by using a sufficiently large area ratio in the blow- 
ing duct diffuser, the term (P^ - P^) in equation(68) can be reduced to almost zero. 
In that case, equation (75) takes the following form: 
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( 88 ) 


A P, = A P , + K„ q 
t ext £ m 

for small values of A /A, 

m b 

The maximum value of entrainment ratio is the value for which equation (89) holds: 

K ! V - A V iP ext < 89) 

In the present example, with = 0. 2 and AP gxt = 0, the limiting value of entrain- 
ment ratio is found by use of figure 30 to be 34. The corresponding mixing tube exit 
Mach number, M , is 0.615. 

The results show that the maximum entrainment ratio which can be 
achieved in a duct system driven by a particular constant-area jet pump is set by 
one of two conditions: 

• by choking at the mixing tube outlet or the suction duct 

inlet (i. e. , M = 1.0 or M n = 1.0) if the duct losses are 
v m s2 ’ 

or sufficiently low 

• by the duct loss limit which is represented by equation 

(89) if the value of M remains below 1.0. 

m 

The form of equation (75) is such that, along curves representing con- 
stant values of , an increase in blowing duct diffusion always yields an increase 
in entrainment ratio until the limiting value is reached. In practice, is a vari- 
able which depends upon the diffuser area ratio. In jet pump systems with low-loss 
inlets, the value of is determined primarily by the blowing duct loss coefficient 
which increases as the area ratio increases. This trend is shown for conical dif- 
fusers in figure 28. Examples of the effect of the variation of are shown by the 
dashed curves in figure 32; these curves represent the loss characteristics of 15° 
and 20° conical diffusers. The peak entrainment ratio for the 20° diffuser is 
achieved by using an area ratio of 3; higher area ratios lead to reduced entrain- 
ment because of increased losses. 
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3.4.7 


Influence of Blowing Duct Area Ratio and Duct Losses Upon Thrust 
Augmentation 


Figure 33 shows the thrust augmentation parameter, r , for the jet pump 
itself. This curve is taken directly from the computer calculations for the jet pump 
selected in section 3. 3. 5. In order to determine the values of the system thrust 
augmentation parameter IT in relation to the blowing duct diffuser area ratio A^A 
and the loss coefficient , figure 32 was used to determine the entrainment ratio 
corresponding to selected values of A^/A^ and K £ . Then figure 33 was used to find 
the associated values of r. The equations of section 3.4.3 permitted calculation of 
T T . 


m 


The variation of system thrust augmentation with blowing duct diffuser 
area ratio is shown in figure 34. The curve for = 0 yields maximum thrust aug- 
mentation when the mixing tube is choked, i.e, for M m = 1.0. Even with a very low 
loss in the duct system (K^ = 0. 1), the thrust augmentation reaches a maximum 
value at a mixing tube Mach number less than 1. 0. The curves for =0.1 and 
0.2 show that the thrust augmentation does not fall off rapidly if the diffuser area 
ratio is made larger than optimum. This suggests that, when designing a duct sys- 
tem without complete data on duct losses, it is preferable to err on the side of in- 
creased diffusion. 


The relationship of system thrust augmentation to the entrainment ratio 
is shown in figure 35. The thrust augmentation peaks on the curve forK^ =0.1 and 
0. 2 and then falls off with increasing entrainment ratio. This is a consequence of 
the fact that the thrust augmentation is proportional to the product of entrainment 
ratio and blowing duct exit velocity as follows: 


IT = 


( m + 1) 



where (m + 1) » m for high entrainment ratio jet pumps 

and V „ « constant 

p2 


so 


TT « mV, 


(90) 
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In order to achieve entrainment ratios higher than the value at the peak of their 
curve, the diffuser area ratio must be increased. This has the effect of reducing 
the exit velocity faster than entrainment increases. The net effect is a reduction 
of the product m and thus TT. 

The maximum entrainment ratio attainable is set by the choking limit 
for the = 0 and = 0. 1 cases. For the = 0.2 case, the maximum entrain- 
ment value is set by the duct loss limit as represented by equation (84) . At this 
limit, the diffuser area ratio is very large and the duct exit velocity is zero. Con- 
sequently, TTmust be zero as shown by equation (85). This illustrates the general 
rule that the thrust augmentation in a jet pump system is always zero at maximum 
entrainment unless the jet pump mixing tube is choked. 

3.4.8 Conclusions 


The two previous sections have shown the influence of duct losses and 
blowing duct diffuser area ratio upon the entrainment ratio and thrust augmentation 
obtained in a jet pump-duct system. The results shown in figures 32, 34, and 35 are 
quantitatively valid only for the particular jet pump geometry and operating conditions 
which were chosen in the section 3. 4. 5. However, the figures illustrate trends which 
are qualitatively correct for high entrainment compressible flow jet pumps as a gen- 
eral class. 


The results show that the design goals of maximum entrainment and maxi- 
mum thrust augmentation may require different duct geometries; a system designed 
for maximum entrainment may have a low value of the thrust augmentation parameter 
and vice versa. The influence of duct losses is shown to be very strong. Entrain- 
ment ratios and thrust augmentation both can be improved significantly by making 
only minor reductions in the duct loss coefficient. This provides considerable in- 
centive for testing flow models of proposed new duct designs in order to adjust their 
geometry to achieve minimum losses. Accurate estimates of duct loss coefficients 
can be obtained from these tests; such estimates are required in order to predict the 
performance of a new jet pump system and to allow selection of the best diffuser area 
ratio. 
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The design problem for a jet pump system often takes the following form: 


Given: Primary flow pressure, temperature, and flow rate 

Duet system inlet and discharge pressure levels and 
inlet pressure 

Problem: What is the proper mixing tube area and blowing 

duct diffuser area ratio to be used to achieve the 
design goal, e. g. , maximum entrainment or thrust 
augmentation ? 

The information provided in plots like figures 32, 34, and 35, together with duct loss 
estimates, will allow the designer to evaluate the effect of diffuser area ratio upon 
entrainment ratio and thrust augmentation for a selected mixing tube area. By pre- 
paring similar sets of curves for several other values of mixing tube area, the de- 
signer can chose the best combination of mixing tube area and diffuser area ratio to 
meet the design goals. New jet pump performance curves analogous to figures 30 and 
32 will be required for each value of mixing tube area to be considered. Data for 
these performance curves can be obtained by using the computer program described 
in section 3.2 A series of computer solutions covering a broad range of jet pump 
geometries and operating conditions is provided in section 3. 3 for use in preliminary 
design calculations. 
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Section 4 


TEST PROGRAM 


The test program had two major objectives: 

• to provide data for use to evaluate the analytical model 

• to determine whether new, reduced-blockage nozzle 
clusters could be used to improve the performance 
of the original jet pump 

The tests covered a range of operating conditions and nozzle geometries as indicated 
below: 


Primary Flow 

pressure range 
temperature range 
nozzle throat area range 
nozzle cluster 
nozzle geometry 

Secondary Flow 

inlet pressure 
inlet temperature 
mixing tube geometry 
pressure rise 

This section of the report describes the jet 
and the results which were obtained. 


55 psia to 400 psia 

200° F to 1200° F 

1. 1 x 10 -4 ft 2 to 6. 0 x 10 -4 ft 2 

three designs 

four designs 


laboratory ambient 
laboratory ambient 

2 

constant area = . 087 ft , two lengths 
regulated by discharge throttling device 

i test arrangement, the test program. 


4. 1 Test Arrangement 

The jet pump test arrangement with its instrumentation is shown schemat- 
ically in figure 36. The primary flow supply system employed a 2-stage reciprocating 
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compressor capable of supplying 7 lbm/min of air at 400 psia. Electrical heaters 
were used to achieve temperatures up to 1200° F. The primary flow was delivered 
to a multiple-nozzle cluster directed along the axis of a constant-area circular mix- 
ing tube. 


The momentum of the primary flow entrains a secondary air flow from 
the room into the bellmouth inlet and then into the mixing tube. Here, the two streams 
mix together and the stagnation pressure of the secondary stream is increased. The 
flow from the mixing tube passes through a conical diffuser and exhausts to the atmos- 
phere through an adjustable throttling cone. 

The individual components of the experimental jet pump are described 

below: 


1. Calibrated bellmouth inlet section. 

This component consists of a wooden bellmouth, metal 
connecting tube, and fiberglass primary flow inlet section. 

The bellmouth differential pressure was calibrated in terms 
of flow rate by using an orifice and blower available in the 
laboratory. The calibrated bellmouth permitted direct meas- 
urement of secondary mass flow rate for all jet pump tests. 

A ceramic insert was used to protect the fiberglass duct 
from the hot primary flow pipe and flange. 

2. Mixing Tube 

The original variable-area mixing tube from the previous 
investigation (ref. 12) was used for the first tests in order 
to provide baseline performance data. This mixing tube had 
a length of 6. 87" (figure 37). After the initial tests were completed , the 
mixing tube was bored out to a constant inner diameter of 4. 000". 

A mixing tube extension of the same diameter was also fabri- 
cated. The remainder of the test program was completed using 
both the original mixing tube length of 6. 87" and the extended 
mixing tube length of 18. 87". 
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3. 


Conical Diffuser 


The initial section of the conical diffuser had a length of 
10. 98" and an area ratio of 1. 79. This diffuser section was 
previously used during the Wagner BLC system tests. An 
additional section was added to this diffuser to obtain an over- 
all area ratio of 5. 0. The purpose of the exhaust diffuser was 
to maximize the static pressure recovery so that the highest 
possible system entrainment ratio could be achieved. Changes 
in the axial positioning of the throttle cone in the diffuser dis- 
charge produced a variable system resistance. The jet pump 
performance characteristic (pressure rise versus entrainment 
ratio) was generated by varying the system resistance in this 
manner. 

4. Nozzle Cluster Geometry 

The nozzle cluster geometry used in the previous investiga- 
tion (ref. 12) was believed to cause excessive blockage of the 
secondary flow at the mixing tube inlet, thus causing reduced 
performance of the jet pump system. Two "reduced-blockage" 
nozzle cluster configurations were tested in order to determine 
whether improved performance could be achieved. 

The first reduced blockage nozzle cluster was made by placing 
2. 0" long nozzle extensions between the original nozzle cluster 
body and each of the nine nozzles. This change moved the clus- 
ter body back away from the mixing tube inlet in order to reduce 
the velocity level around the cluster by increasing the adjacent 
flow area. This change is shown on figure 38. 

The second reduced blockage nozzle cluster was a completely 
new design consisting of 7 nozzles at the end of small diameter 
tubes (1/4" dia.). This low drag cluster is shown on figure 39. 
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5. Primary Nozzle Geometry 

Five sets of primary flow nozzles were used in the test program. 
The throat size, nozzle type, and design conditions are listed in 
table 3. 


The first four nozzle sets listed in table 3 were used with the original 
nozzle cluster and also the first reduced-blockage configuration (original cluster 
with 2" extension tubes). The fifth set, Case LD#2-4, consisted of seven nozzles 
for the second reduced-blockage cluster. The Case LD#2-4 nozzles were designed 
for the same operating conditions as the Case 4 nozzles of the original cluster. 


The nozzle flow coefficients listed in the table were calculated from test 
results according to the definition below: 


W 


w 


W 


P ideal 


where 



measured nozzle flow rate at design pressure and 
temperature 


P ideal 


isentropic flow rate through nozzle throat at design 
pressure and temperature; based upon one-dimensional 
flow assumption 


4 . 2 Instrumentation and Data Reduction Pro cedures 
4. 2. 1 Instrumentation 

The instrumentation used to determine the performance of the experi- 
mental jet pump is shown on figure 40 and described in table 4. 

The jet pump inlet bellmouth was calibrated for use as a flowmeter. The 
calibration was accomplished by connecting the bellmouth and the suction duct to the 
inlet of a blower. An orifice and a throttling arrangement were included in the blower 
system. The blower permitted calibration of the bellmouth up to a flow rate of 200 
lbm/min. The resulting bellmouth flow equation is given on the following page: 
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229.5 Yp b 


(91) 



W 


where 


Ah^ = P^ differential pressure, inches of water gage 

3 

= inlet density, lbm/ft 

Stagnation pressure traverses were made in the mixing tube and diffuser 
with a Kiel probe. When the short mixing tube configuration was used, traverses 
were taken only in the diffuser at the location shown in figure 40 . When the mixing 
tube extension was used, stagnation pressure traverses were taken in the tube 16.4 
inches downstream from the primary nozzle exit plane. Additional traverses were 
made in the diffuser 26. 4 inches downstream from the primary nozzle exit plane. 

The angular orientation of the traverse locations is shown on figure 40 . 
The same numbering system was used for all traverse locations. 


4.2.2 Data Reduction Procedures 


The measured data were used to calculate the following jet pump per- 
formance parameters. 


m = 


AP 


s 


W 

s 


W 

P 



jet pump entrainment ratio 


jet pump static pressure 
parameter 



P - P 
mo so 


jet pump stagnation pressure 
rise parameter 


V vs. 



po 


velocity profiles 


P vs. distance 
P vs. m and P 


jet pump static pressures 
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The static pressure parameter was calculated using the wall static pres- 
sure measured at the discharge of the short mixing tube or the maximum static pres- 
sure reached within the extended mixing tube. The secondary flow stagnation pres- 
sure at the mixing tube inlet (P ) was determined by subtracting the measured bell- 
mouth inlet and cluster loss (section 4. 2. 3) from the barometric pressure. 


* 

The jet pump stagnation pressure rise parameter, A P^. , is based upon 
the secondary flow total pressure P gQ and upon P . Our first tests showed that it 
was not possible to measure P mo accurately with the "short" mixing tube, i.e. , the 
original 6. 87" mixing tube length. Complete mixing was not achieved by the end of 
this short tube; sharp velocity peaks corresponding to the primary flow jets were ob- 
served at the exit section. These peaks were too sharp to be accurately measured 
by stagnation pressure probes of reasonable size. Therefore, the stagnation probe 

measurements could not be used to determine P with the necessary accuracy. 

mo 

To obtain an approximate value of P mQ for the short mixing tube tests we 
used the following procedure: 


Stagnation pressure probe traverse data was taken at the conical 
diffuser exit section, P^ q in figure 40. This data was used to de- 
termine an average stagnation pressure (P^ ) by procedures des- 
cribed below. Then the value of P mQ was computed by using equ- 
ation (92) which accounts for the diffuser stagnation pressure loss. 


P = P , + K 
mo do 


p V 

m m 


2g. 


The value of K was selected to be 0. 05, a value representative of 
the loss coefficient for a good diffuser. 


(92) 


The procedure used to determine an average stagnation pressure 

from the stagnation pressure probe readings was the "mass-momentum 

integral method". The stagnation pressure was measured along a diameter 

at the diffuser traverse location. A wall static pressure measurement 

was obtained at that cross-section. Using the ratio of the local stagnation 

pressure and the wall static pressure at each point along the diameter, 

2 

the local Mach number was determined. Using a plot of (1+ kM ) vs. 

2 2 

cross-section area (i.e. , r /r ), the impulse function was determined 
by graphical integration: 
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(93) 


V 


impulse function 



+ kM 2 ) dA 


The following equations were used to determine the "mass momentum" 
averaged properties of this non-uniform compressible flow: 


- calculate M from: — — 

1 + kM 2 



<P 


(94) 


- calculate P (static) from: P = 

(1 + kM ) A 


(95) 


- calculate P (stagnation) 
from: 0 


P 

o 


P 



k 

k-1 


(96) 


These "mass-momentum average" values of pressure satisfy the 
measured mass flow rate and integrated momentum of the flow. 
They correspond to values that would be obtained if the actual non- 
uniform flow was mixed to a uniform flow in a frictionless, con- 
stant area duct. 


To calculate P for the extended mixing tube, the "mass momentum 
mo 

integral method" was applied to the traverse data taken in the mixing tube 16.4" 
downstream from the primary nozzle exit plane. No correction for diffuser losses 
was needed in this case. 

The velocity profiles were calculated from the local values of Mach num- 
ber and the measured jet pump exhaust temperature (T ). 

4.2.3 Suction Duct and Nozzle Cluster Losses 

In order to calculate the jet pump pressure rise parameters AP and 
* s 
A P, , the mixing tube inlet stagnation pressure P is required. The value of P 
t SO so 

was determined by subtracting the appropriate suction duct and nozzle cluster 
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stagnation pressure losses from barometric pressure (the suction duct inlet stagna- 
tion pressure). The loss data used for the determination of P go is given below. 

Bellmouth Inlet 


The stagnation pressure losses in the bellmouth, suction duct, and nozzle 
cluster configuration used in the test rig (figure 36) were measured for the following 
four configurations: 

1. Original nozzle cluster 

2. Original nozzle cluster with 2" nozzle extensions (figure 38) 

3. "Low-drag" nozzle cluster (figure 39) 

4. No nozzle cluster 

The loss measurements were made by connecting the jet pump including the initial 
section of the conical diffuser (area ratio 1.79) to the suction line of a blower and 
orifice installation. Air was drawn through the bellmouth and jet pump system by 
the blower. The stagnation pressure at the end of the constant area mixing tube was 
calculated using the measured wall static pressure and the Mach number (computed 
on a one-dimensional basis) at the same measuring station. The stagnation pressure 
loss was set equal to the difference between this value and the atmospheric pressure. 
The loss values for the four nozzle configurations are presented on figure 41. These 
loss values include the suction duct loss, the nozzle cluster loss, and the short mix- 
ing tube wall friction loss (L/d m = 1. 35). 

The results show that the losses are identical for the two reduced block- 
age clusters #2 and 3. A comparison of the losses for the low blockage clusters to 
the losses for the original cluster and the minimum possible losses (the "no cluster" 
curve) shows that the reduced blockage clusters cut the cluster losses by about 42^. 

These results were used for two purposes: 

1. to guide the selection of the "optimum" nozzle cluster 
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2. to permit calculation of P gQ at the jet pump inlet by sub- 
tracting the suction tube and nozzle cluster losses from 
the atmospheric pressure. The P gQ value was used to 
calculate the static pressure parameter AP g * and the 
stagnation pressure rise parameter A Pj*. 

BLC Suction Duct 


The data obtained during tests of the Wagner BLC system under the pre- 
vious contract NAS 2-2518 indicated that the aerodynamic drag of the present nozzle 
cluster was much higher than had been expected. This drag was considered to be 
one of the principal reasons for the difference between the predicted and actual per- 
formance of the jet pump BLC system. To check this point, the NAS 2-2518 suction 
duct (figure 42), the short mixing tube, and the initial section of the conical diffuser 
were connected to the inlet of the blower and orifice system. The original nozzle 
cluster geometry was installed in the suction duct. The tests results are shown on 
figure 41. The BLC suction duct produced higher losses than the bellmouth config- 
urations. 


Still using the NAS 2-2518 suction duct, the nozzle cluster was dismantled 
in 3 steps with a loss test made between each step. The nine nozzles only were re- 
moved first but no decrease in loss was measured. Next, the cluster body was re- 
moved but the elbow was left in place. The measured losses dropped by 13% to the 
"no cluster" curve shown on figure 41. This amounts to a 5 psf loss reduction fora 

jet pump entrainment ratio of 20 (W = 128 lb/min). When the elbow was removed, 

s 

no further reduction of losses was observed. Therefore, the cluster body with its 
extensions produced the nozzle cluster losses. 

This series of loss tests shows that the nozzle cluster losses are much 
smaller than had been deduced from the experimental data obtained on the jet pump 
BLC system. The NAS 2-2518 final report showed suction duct losses to be three 
times the values shown on figure 41. At that time, a large percentage of these losses 
were attributed to the nozzle cluster. The results in figure 41 show that this conclu- 
sion was incorrect for two reasons: 
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1. The calculation of suction duct stagnation pressure loss from 
the BLC system data is sensitive to the choice of the static 
pressure value at the mixing tube inlet to be used for the cal- 
culation. The static pressure data recorded by the first 2 or 
3 pressure taps in the first half-inch of the mixing tube always 
showed a sharp spike of low static pressure at the mixing tube 
inlet. The static pressure variation becomes more gradual 
from about the fourth tap onward in the mixing tube. Testing 
of the mixing tube with the BLC suction duct during the series 
of loss tests of figure 41 showed this same variation of static 
pressure even without the nozzle cluster in place. The spike 
of static pressure must be caused by the local curvature of the 
streamlines as the flow reaches the mixing tube throat. If suc- 
tion duct losses are calculated from the third or fourth static 
pressure tap instead of the first tap measurement, the calcu- 
lated BLC suction duct losses would compare favorably to the 
measured losses shown on figure 41 . 

2. The loss tests on the BLC suction duct showed that the duct it- 
self causes most of the suction duct losses. The nozzle cluster 
accounts for about 13% of the suction duct loss which may amount 
to 4% to 8% of the total system resistance. 

4.2.4 Pressure Loss Due to Wall Friction in the Constant-Area Mixing Tube 

The mixing tube wall friction loss levels must be taken into account when 
comparing the experimental pressure parameters to the analytical predictions which 
do not include these losses. Therefore, the pressure loss in the extended mixing tube 
was measured during the nozzle cluster loss tests which were described in section 
4. 2. 3. The wall static pressure was measured at two points 10. 88 inches apart in 
the constant 4" diameter mixing tube while air was being drawn through the tube by 
the blower. The total pressure change was calculated using the local Mach numbers 
(determined on a one-dimensional flow basis) at the two cross sections. The static 
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pressure change and the stagnation pressure loss as a function of flow rate are given 
in figure 43. The static pressure change is larger than the stagnation pressure loss 
because of compressibility effects. The measured loss levels agree well with pre- 
dictions based upon pipe friction factors for fully-developed turbulent flow. 

When the jet pump is operating, the mixing tube velocity profiles differ 
from the profiles for fully-developed turbulent flow in pipes because of the primary 
flow - secondary flow interactions. Therefore, wall friction losses in a jet pump 
mixing tube can be expected to differ somewhat from the losses predicted for fully- 
developed turbulent flow. No data correlations or analytical procedures are avail- 
able to allow accurate prediction of wall friction losses in mixing tubes. Thus, pipe 
friction factor correlations or equivalent test results as in figure 43 must be used as 
a first approximation in order to estimate mixing tube wall friction losses so that the 
analytical predictions of jet pump performance can be compared to the experimental 
resu Its. 

4. 3 Tabulation of Test Conditions 

The jet pump testing was carried out in three series of tests. The con- 
figurations used and the test objectives are described below. 

Series 1: The jet pump included the calibrated bellmouth inlet, 

the original short NAS 2-2518 mixing tube with varying area, 
and the original nozzle cluster. The purpose of this test series 
was to determine the performance of the jet pump as used in the 
NAS 2-2518 BLC system tests. 

Series 2: The jet pump included the calibrated bellmouth inlet, 

the short constant-area mixing tube, and three nozzle clusters; 
the original NAS 2-2518 cluster, the same cluster with 2" noz- 
zle extensions and moved back from the mixing tube throat (i.e. , 
low-drag cluster #1), and the low-drag nozzle cluster #2 as in 
figure 39. The objectives of this test series were to determine 
the effect of nozzle cluster design upon jet pump performance 
and to determine which cluster would be best for subsequent 
testing. 
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Series 3: The jet pump configurations tested included the calibrated 

bellmouth inlet, both the short and the extended constant-area 
mixing tubes, and the optimum nozzle cluster selected from the 
results of test series #2 (the LD#1 cluster). The purpose of this 
test series was to provide experimental verification of the jet 
pump performance analysis over a broad range of operating con- 
tions. Each of the four nozzle sets described in table 3 were 
used; they were operated not only at their design points but also 
at off-design pressures and temperatures. 

A description of the individual tests, operating conditions 
and configurations is given in table 5. 

4. 3. 1 Presentation of Data 

Tabulated data and graphical results are presented for each of the run 
numbers listed in table 5 . An index to the tables and figures is given in table 6 . 
The test data is presented in tables 8 to 19 and figures 44 to 84. 

4. 4 Discussion of Test Results 

4.4.1 Comparison of Constant- Area and NAS 2-2518 Varying-Area Mixing Tubes 

The original NAS 2-2518 short mixing tube was designed to have a 5. 5% 
contraction in area from inlet to outlet. The tube was intended to have constant- 
pressure mixing at its design point. Since the jet pump has not been able to produce 
a flow rate as high as the mixing tube design flow rate, the constant static pressure 
condition was never achieved. However, a nearly-constant static pressure distribu- 
tion has been approached near the discharge end of the mixing tube at the higher flow 
rates obtained during the test program. 

The short constant-area mixing tube was made by boring out the NAS 2- 
2518 mixing tube to a constant internal diameter of 4.00". 
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Entrainment Ratio Results 


Tables 8 and 9 present the primary and secondary flow rates and the 
entrainment ratio for each test point recorded for the two short mixing tube config- 
urations. For each primary stagnation pressure, the maximum entrainment ratio 
represents the test where the throttle cone was moved out of the exhaust diffuser as 
far as it would go. The two test configurations were then completely identical except 
for the mixing tube. The data shows that use of the NAS 2-2518 mixing tube resulted 
in maximum entrainment ratios about 6% to 8% higher than the constant-area mixing 
tube configuration. A basic difference was observed between the measured velocity 
profiles at the diffuser discharge (see below). 

Jet P ump Sta gnation an d Sta tic Pressure Parameter Results 

The experimental results for both mixing tube configurations are com- 
pared to the analytical predictions in figure 48. Both the stagnation and static pres- 
sure parameters are plotted for four primary pressures. 

The experimental stagnation pressure rise values calculated from a sin- 
gle traverse were found to be 13% to 32% below the analytical predictions. The one 
experimental value calculated from traverses in two perpendicular directions was 
8% below the analytical prediction. The differences between the experimental and 
analytical values are caused by a combination of effects listed below. 

• The experimental value is very sensitive to the accuracy 

and thoroughness of the traversing of the diffuser discharge 

to determine P, . 

do 

• The estimated diffuser losses between the mixing duct exit 

and the traverse station are added to P , to determine P 

do mo 

Since the diffuser inlet flow is highly distorted (the mixing 
tube is too short and primary nozzle jets persist into the dif- 
fuser inlet) the estimated losses may be too low. 

• The suction duct and nozzle cluster losses must be subtracted 
from the barometric pressure in order to determine P gQ . The 
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loss characteristics used were obtained by drawing air past 
the nozzle cluster using a blower at the discharge of the duct 
system (section 4. 2. 3). The flow around the cluster may be 
different when the primary nozzles are in operation. Higher 
inlet losses may be the result. 

• Wall friction losses in the mixing tube (not taken into account 
in the analytical predictions) reduce the jet pump stagnation 
pressure rise (section 4.2.4). 

The experimental static pressure parameter curves are similar in slope 
to the curves for the analytical predictions. The experimental curves are shifted to 
higher negative values of the static pressure parameter. This means that the meas- 
ured static pressure values are too low just as the integrated stagnation pressure 
values were too low. There are two reasons why this occurred. 

• The data shows that the static pressure is still increasing at 
the end of the constant-area mixing tube. This means that con- 
siderable mixing is still taking place. A longer constant-area 
mixing tube improves the agreement of the experimental and 
analytical static pressure parameters as shown in section4.4. 3. 

• The suction pipe and nozzle cluster losses, which are used to 
determine P , may not be accurately represented by our 

SO 

measured loss characteristics. This problem was discussed 
above for the stagnation pressure loss parameter. 

To improve the agreement between the experimental and analytical re- 
sults, the length of the mixing tube was increased using the mixing tube extension 
piece. The stagnation pressure traverses were then repeated at the end of the ex- 
tension tube. 

Velocity Profiles at the Diffuser Discharge 

Velocity profiles were calculated from the stagnation pressure traverse 
data used to calculate the jet pump discharge stagnation pressure. The velocity pro- 
files measured for the two mixing tubes are presented in figures 44, 45, and 49. 
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Figure 44 shows velocity profiles taken at four throttle settings at a pri- 
mary flow pressure of 350 psia using the NAS 2-2518 mixing tube. Three observa- 
tions can be made from this figure. 

• The profiles are generally not symmetrical, 

• a sharp dip in velocity occurs at the center of the flow 
for the higher entrainment ratios, 

• the velocity dip becomes less pronounced as entrainment 
ratio decreases. 

Figure 45 shows two velocity profiles for the NAS 2-2518 mixing tube at 
reduced pressures and one velocity profile measured with an unheated primary flow. 
In all profiles, a sharp dip in velocity occurs in the center of the tube. The unheated 
primary flow velocity profile is nearly symmetrical. The asymmetry of all of the 
heated primary flow velocity profiles is probably caused by a radial shift of the noz- 
zle cluster relative to the mixing tube due to thermal expansion of the elbow which 
feeds hot air to the cluster. This shift has been measured to be about 1/ ^6 inch. 

This situation was corrected in some of the subsequent tests by setting the nozzle 
cluster off-center at room temperature to compensate for thermal expansion at op- 
erating temperature. 

Figure 49 presents velocity profiles measured for the short constant-area 
mixing tube configuration. The velocity profiles are not symmetrical for this mixing 
tube either. For run 24, two perpendicular velocity traverses were obtained. Both 
traverses were asymmetrical and both were shifted in the direction which would be 
expected if caused by elbow thermal expansion (the location of the traverse planes 
relative to the elbow is shown in figure 40). The dip in velocity at the center of the 
tube was much smaller for the constant area mixing tube configuration than for the 
NAS 2-2518 mixing tube. Further discussion of the dip in the velocity profile is pre- 
sented in section 4. 4. 4. The presence of the exhaust cone was shown to have no ef- 
fect on the velocity profile and a wedge probe traverse revealed that a small amount 
of swirl was present in the flow. 





Mixing Tube Static Pressure Distribution 

Static pressure data is presented in figures 46, 47, and 50. Figure 47 
compares the variation of static pressure along the mixing tube and diffuser for three 
configurations: 

1. NAS 2-2518 mixing tube in the BLC duct system. 

2. NAS 2-2518 mixing tube in the Dynatech Test Rig. 

3. Constant area mixing tube in the Dynatech Test Rig. 

The test conditions for each of the three configurations are listed below. 
Slight differences in primary temperature and entrainment ratio existed for the three 
conditions. 


Configuration 

P 

po 

psia 

T 

po 
° F 

W 

P 

lbm/min 

W 

s 

lbm/min 

W 

m 

lbm/min 

m 

1 

350 

1200° 

6.4 

118.5 

124. 9 

18.5 

2 

350 

1085° 

6. 70 

118. 20 

124.9 

17.65 

3 

350 

1130° 

6.65 

118. 25 

124. 9 

17. 8 


The static pressures for configuration 2 are about 2" of water larger 
than for configuration 1. This shift in the static pressure level is a result of re- 
duced losses in the suction duct and bellmouth as compared to the BLC suction duct. 
The 2" of water shift is equivalent to a 10. 5 psf decrease in suction duct losses. Com- 
parison of this number to figure 41 shows that 10. 5 psf is about half of the difference 
between the measured losses for the two inlet geometries. 

The effect of changing from the contracting NAS 2-2518 mixing tube to a 
constant area tube is shown by configurations 2 and 3 in figure 47. The static pres- 
sure rises more rapidly in the constant area tube because of the lower velocity levels 
However, the static pressure recovery obtained in the diffuser is substantially larger 
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for the NAS 2-2518 mixing tube case. This effect may be the result of lower mixing 
tube stagnation pressure rise or reduced diffuser effectiveness for the constant-area 
mixing tube case. 

Figures 46 and 50 show the variation of static pressure at two locations 
in the mixing tube for configurations 2 and 3 as a function of the primary pressure 
and the entrainment ratio. The increase in static pressure for the constant area 
tube above that of the contracting tube was found to exist for all primary pressure 
levels and entrainment ratios tested. 

4.4.2 Reduced Blockage Nozzle Clusters 

Two reduced drag nozzle cluster configurations were tested to determine 
the influence of the cluster drag upon jet pump performance and to allow selection of 
the best cluster for further testing. Section 4. 1 presents the dimensions of each of 
the cluster configurations. Table 5 lists the test conditions and table 6 provides 
an index to the results obtained. 

The jet pump performance with the three nozzle cluster configurations 
(original, original cluster with 2-inch nozzle extension, and the low-drag cluster #2) 
were compared in the following ways: 

1. System performance at minimum throttling (wide 
open throttle cone) 

2. Static pressure parameter 

3. Static pressure variation along the mixing tube at 

a selected total flowrate 

4. Velocity profiles 

5. Cluster pressure loss characteristics 
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Entrainment Ratio Results 


The system performance at minimum throttling is presented in figure 51 
where secondary flow rate is plotted as a function of primary flow rate. Except for 
the loss characteristics of the nozzle clusters themselves, the system loss charac- 
teristics are identical for all three nozzle configurations when the cone is in the wide 
open position. Both of the modified nozzle clusters show equal improvement with 
respect to the performance of the original cluster. Both have an increased second- 
ary flow rate for the same primary flow rate. 

The LD #2-4 cluster nozzles were designed for a slightly lower pri- 
mary flow rate at a given pressure and temperature than the Case 4 Nozzles. This 
difference in design accounts for the shift in data points along the curve for the two 
reduced drag clusters. 

Jet Pump Static Pressure Parameter Results 

The analytical and experimental static pressure parameters are com- 
pared for the three nozzle clusters in the following figures: 

Original Nozzle Cluster: figure 48 

LD#1 - Low- Drag Nozzle Cluster: figure 52 

LD#2 - Low-Drag Nozzle Cluster: figure 56 

The best agreement between analytical and experimental results was obtained for the 
LD#1 cluster (figure 52). The higher nozzle cluster losses in the original cluster 
(figure 48) and the less-complete mixing obtained with the 7-nozzle arrangement of 
the LD#2 cluster (figure 56) are the causes of the poorer agreement between analy- 
sis and test results for these clusters. 

Velocity Profiles at the Diffuser Discharge 

The velocity profiles for the original cluster without and with nozzle ex- 
tensions are presented in figures 49 and 53. The profiles had two velocity peaks with 
a dip in the center amounting to 30 to 60 fps. The velocity profiles for the LD#2 
cluster are presented in figure 57. Three nearly-equal velocity peaks were obtained 
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with the LD#2 cluster when traversing along a line passing through three nozzles, 
and one velocity peak was obtained when traversing along a line passing between noz- 
zles. The difference between the profiles for the original cluster and for the LD#2 
cluster appears to be related to the number of nozzles used. The LD#2 nozzle clus- 
ter has seven equally-spaced nozzles, each having an equal share of the mixing tube 
flow cross section to energize. The original cluster has two additional nozzles 
around the outside, leaving the center nozzle with a larger percentage of the mixing 
tube flow area to energize. The velocity dip probably was not caused by nozzle clus- 
ter losses (i.e. , a wake effect) because the LD#1 and LD#2 clusters apparently had 
similar loss characteristics (figure 51), 

Several tests were completed with and without the throttle cone in place. 
No change in the velocity profile was detected. Therefore, the presence of the throt- 
tle cone at the diffuser exit does not seem to be the cause of the dip at the center of 
the velocity profile when the original nozzle cluster was used. No satisfactory ex- 
planation for the dip was developed during this program. 

Mixing Tube Static Pressure Distribution 

The static pressure variation along the mixing tube is shown for all three 
clusters in figure 55. The static pressure levels are dictated primarily by the total 
flow rate, but they are also slightly affected by changes in the mixing process which 
accompany modifications of the nozzle cluster geometry. The original nozzle clus- 
ter produced similar pressure distributions with and without the nozzle extensions. 
The pressure distribution for the LD#2 cluster shows higher static pressures near 
the mixing tube inlet because of reduced blockage, and lower static pressures at the 
mixing tube exit because of less complete mixing (the result of the reduction in num- 
ber of nozzles from 9 to 7). 

C luster Drag 

The cluster loss tests were discussed in section 4. 2. 3 . Both of the low- 
drag clusters showed equal reduction in loss characteristics relative to the original 
cluster (figure 41). 



Selection of the Optimum Nozzle Cluster 




The measured nozzle cluster losses (figure 41) and the jet pump per- 
formance characteristics with minimum throttling (figure 51) show that the two low- 
drag nozzle clusters yield similar results. The comparison of experimental and 
analytical static pressure parameters and the comparison of static pressure distri- 
butions along the mixing tube for both low-drag clusters show that the original clus- 
ter with 2 inch nozzle extensions (LD#1) is slightly preferable to the LD#2 cluster. 

The original cluster with 2 inch nozzle extensions (LD#1) was selected 
as the best cluster for the remainder of the test program for the following reasons: 

1. The nozzles for the remainder of the test program 
were already available for this cluster. 

2. This cluster gives better agreement between experi- 
mental and analytical static pressure values because 
of the more complete mixing upstream of the measur- 
ing station. 

3. The LD #2 cluster offers no advantages in comparison 
to the original cluster with extensions; the measured 
losses are equal for both clusters and the flow rate 
curves at the wide-open throttle position are the same. 

4.4.3 Comparison of Short and Extended Mixing Tubes 

Performance data was obtained for both the short and the extended 
constant-area mixing tubes using each of the four nozzle sets with the LD#1 nozzle 
cluster. The test results are indexed in table 6 . 

The extended mixing tube was found to change the jet pump performance 
in the following ways: 
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1. The static and stagnation pressures both reached a 
maximum value within the extended mixing tube. 



This indicates that complete mixing was achieved. The 
extended mixing tube was superior to the short mixing 
tube, particularly with respect to static pressure recovery 
at low entrainment ratios. 

2. The wall friction losses in the extended mixing tube were 
significant in comparison to the stagnation pressure rise 
developed by the jet pump. 

3. The velocity profile at the diffuser inlet was improved by 
extending the mixing tube. 

Identical duct loss characteristics existed for the maximum-entrainment 
ratio runs for each test number. The throttle cone was withdrawn to a fixed location 
for these runs. Thus, the performance characteristics ( of the jet pump as influenced 
by nozzle design, mixing tube length, and primary flow pressure and temperature 
can be determined by comparing these maximum-entrainment runs. 

Entrainment Ratio Results 

A Comparison of the maximum entrainment ratio achieved with the long 
and short mixing tubes operated with the same nozzles and the same primary flow 
conditions shows that the short mixing tube configuration achieves a slightly larger 
entrainment ratio in all cases. The mixing tube extension section allows more com- 
plete mixing and improves the diffuser inlet velocity profile. However, the extended 
tube introduces an additional frictional loss which becomes significant at large flow 
rates (figure 43). The reduction in entrainment for the extended mixing tube shows 
that the extra friction losses in the longer tube cancel the effects of improved mix- 
ing. The use of a mixing tube length longer than the short tube and shorter than the 
extended tube would probably lead to a higher entrainment ratio than was developed 
by either of the tested lengths. 
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Jet Pump Static Pressure Parameter Results 


Table 6 lists the figure numbers which show the static pressure para- 
meters for the eight nozzle and duct configurations tested. The experimental values 
shown are based upon measured static pressures (P m or P max ) and upon secondary 
flow stagnation pressures (P ) calculated from the measured inlet losses, cluster 
losses, and mixing tube losses (figures 41 and 43). The analytical values shown in 
the figures were calculated by the computer using the ideal jet pump analytical model 
which neglected inlet, cluster, and mixing tube losses. 

The results show that the extended mixing tube test data corresponds 
more closely to the analtycal predictions. The biggest difference between the short 
and long mixing tubes occurs at low entrainment ratios because the extended mixing 
tube significantly increases the static and stagnation pressure recovery at low flow 
rates. At high entrainment ratios, no improvement is produced by the extended mix- 
ing tube because the increased frictional losses cancel the potential gains from more 
thorough mixing. 

Inclusion of the inlet and mixing tube losses in the analytical model would 
give lower (more negative) values of the static pressure parameter which would agree 
more closely with the experimental values. The biggest changes would occur at high 
entrainment ratios where the curves are presently furthest apart. 

Jet Pump Stagnation Pressure Parameter Results 

Table 7 lists all of the experimental stagnation pressure rise parameters 
which were calculated from the traverse data taken. Traverses were made either 
at the mixing tube discharge or in the conical diffuser at a station where the area is 
1. 61 x the inlet area. The stagnation pressure at the traverse station was calculated 
by the mass-momentum method presented in section 4. 2. 2. Each result represents 
the integration of one or two traverses as indicated in the table. The value of P mo 
was determined from the diffuser exit traverse values (P^ Q ) by correcting for the 
diffuser losses using equation (92) of section 4. 2. 2. 
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Figure 59 compares the analytical predictions and experimental values of 
the stagnation pressure rise parameter for the extended mixing tube tested with the 
Case 4 nozzles at 300 psia and 1150° F primary flow conditions. The experimental 
values were obtained as follows: 

'0 The stagnation pressure value P mQ was calculated 
using the diffuser exit stagnation pressure traverse 
data and correcting for diffuser losses by using 
equation (92) of section 4. 2. 2. 

( if 0 The stagnation pressure P mo was calculated by 
using stagnation pressure traverse data obtained 
at the mixing tube exit. 

^3 An additional correction was made to the^'Q data 
to account for mixing tube wall friction effects. The 
wall friction pressure losses are taken from figure 43 
to make this correction. 

$ ^ The wall friction pressure losses in the mixing tube 

(figure 43) were added to the Q^data. 

The. correction of the stagnation pressure rise parameter to account for duct losses 
places three of the four experimental points within 6 psf of the analytical value. The 
differences which remain may be due to the factors discussed in section 4.4.1 where 
similar results were presented for the short mixing tube test. 

Additional comparisons of the analytical and experimental stagnation pres- 
sure values are made in section 5.1 of this report. 

Velocity Profiles at the Diffuser Discharge 

Table 6 lists the figure numbers showing the velocity profiles obtained 
with both the short and long mixing tubes. The effect of extending the mixing tube 
can be seen clearly on figure 73 where results for both mixing tubes are plotted for 
the Case 2 nozzles. The centerline depression at the diffuser discharge is much 
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smaller when the extended mixing tube is used since the flow has a longer flow path 
for mixing. The tests with the Case 3 and Case 4 nozzles gave the same result. The 
longer mixing tube improves the symmetry of the flow entering the diffuser and re- 
duces the likelihood of flow separation in the diffuser. 

For the Case 3 nozzles with the extended mixing tube, two traverses 
each were taken at the mixing tube exit and at the diffuser exit. One traverse at 
each location passed along a diameter through the wake of three nozzles. The sec- 
ond traverse in each location passed through the wake of only the centerline nozzle 
(figure 79). The four traverses show the flow to be reasonably symmetrical. A re- 
duced velocity exists along the centerline while two peaks of velocity appear on ei- 
ther side. The velocity profiles in the two planes at one station were nearly iden- 
tical. Further discussion of the velocity profiles is included in section 4.4.4. 

Mixing Tube Static Pressure Distribution 

Table 6 lists the figures showing the static pressure variation along the 
mixing tube. These variations are plotted only for the extended mixing tube tests. 

The static pressure at each location in the duct is a result of the interaction of the 
following factors: 

1. The local Mach number of the flow, 

2. The stagnation pressure rise achieved by mixing of 
the two streams, 

3. The frictional losses on the walls. 

The data in figures 71, 75, and 81 show that the static pressure reaches a maximum 
at the middle or near the end of the mixing tube extension. In the duct upstream of 
the location of maximum static pressure, the increase in stagnation pressure due to 
mixing is larger than the reduction in stagnation pressure due to wall friction. Be- 
yond the maximum point, the added wall friction loss becomes dominant. The loca- 
tion of the maximum static pressure point is closer to the mixing tube inlet for the 
high flow rate test points because of the increased magnitude of frictional pressure 
loss which accompanies the increased velocity levels. 
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These results show that an optimum mixing tube length exists which will 
produce the maximum static pressure. This optimum length is a function of flow 
rate. 

4.4.4 Velocity Profile Investigations 

As mentioned in other portions of section 4.4, the velocity profiles 
measured at the mixing tube exit and the discharge of the initial diffuser section have 
shown various degrees of distortion. This distortion is due to the following causes: 

1. Angular and/or radial misalignment of the nozzle cluster 
and mixing tube centerlines. 

2. Non-uniform pumping. 

The alignment of the nozzle cluster elbow with the mixing tube center- 
line was accomplished with a special centering plug which was inserted snugly into 
both the elbow and the mixing tube. The original centering plug aligned the elbow 
and tube concentrically with the elbow at room temperature. Traverses taken with 
low primary air temperatures (from ambient to 200° F) gave velocity profiles which 
were quite symmetrical about the centerline, thus indicating good alignment. How- 
ever, when the primary flow was heated to 1150° F, the elbow flange-to-centerline 
dimension increased by about 1/16" due to thermal expansion. The high velocity 
region in the velocity profile shifted noticeably in the same direction (an example 
is given by figure49 ). 

To compensate for the thermal expansion, an offset plug was made to 
position the two centerlines 1/ 16 inch apart when the elbow was at room tempera- 
ture. The offset plug was used to position the cluster for all of the high tempera- 
ture tests performed subsequent to test No. 16. Even with the offset plug, some of 
the short constant area mixing tube tests showed distorted velocity profiles, indicat- 
ing that the mixing process in the short tube is highly sensitive to slight misalign- 
ments. All of the velocity profiles for the extended mixing tube tests were reason- 
ably symmetrical showing that the mixing process in the longer tube is relatively 
insensitive to misalignment. 
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Most of the velocity profiles exhibited a slight depression of the center- 
line velocity below the velocity of the surrounding flow. This depression was less 
pronounced when the extended mixing tube was used (an example is given by figure 
73). To determine the reason br the slight depression of velocity along the center- 
line, the following series of tests were run: 

1. Velocity profiles were obtained with the Kiel probe with 
and without the throttle cone in place. The same primary 
flow conditions were maintained for both tests. 

2. A wedge probe was used to measure local static pressure 
and flow direction for comparison with a Kiel probe measure- 
ment. 

3. A nozzle arrangement using 8 Case 2 nozzles around the cir- 
cumference of the cluster and 1 Case 4 nozzle located in the 
center was tested to explore the effect of increasing the mo- 
mentum in the center of the mixing tube. 

The first test described above showed conclusively that the throttle cone 
does not influence the flow at the traverse station. There was no detectable dif- 
ference in the traverse results with and without the cone. 

The traverse results for tests 2 and 3 above are shown on figure 85. The 
wedge probe results, like the Kiel probe results, show the depression in velocity 
along the centerline. Associated with this depression in velocity was a reduction in 
static pressure of about 1 inch of water and a departure of the velocity from the 
axial direction by roughly + 2 ° all along the diameter. These results indicate that 
a slight swirl exists in the flow. 

The third test above employed a center nozzle with a throat area about 
50% larger than that of the surrounding eight nozzles. The velocity profile still 
shows a centerline depression. Thus, the depression cannot be accounted for as 
only a primary flow momentum deficiency in the center of the mixing tube. 
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The small amount of swirl that appears to be present in the flow is not 
likely to have a significant effect on the test results. This small amount of swirl 
could be caused by one or more of the primary nozzles being bent at a small angle 
to the axis, by the presence of the primary nozzle elbow, or by wakes shed off of 
objects in the laboratory outside the jet pump. 

No satisfactory explanation for the centerline depression in the velocity 
profile remained after these tests were completed. The cause of the depression is 
unknown. 
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Section 5 


COMPARISON OF ANALYTICAL AND EXPERIMENTAL RESULTS 


5. 1 Jet Pump Stagnation Pressure Rise 

The experimental measurements of the jet pump stagnation pressure 
rise parameter, AP^.*, generally fell below the analytical predictions. Examples 
are given in section 4. 4. 1 (figure 48) and section 4. 4. 3 (figure 59). Additional 
values of A.P^* were measured in other tests and the results are given in table 7. 

The analytical values of AP^.* were predicted by neglecting wall friction 
losses in the mixing tube, by assuming complete mixing, and by neglecting suction 
duct losses and conical diffuser losses. Thus, the analytical values represent 
"ideal" jet pump performance. The measured performance fell below the "ideal" 
values for the following reasons: 

• the experimental value is very sensitive to traversing 
thoroughness and accuracy 

• mixing tube wall friction losses were not included in 
the analysis; these losses are not negligible 

• in many tests, the stagnation pressure traverses were 
made in the conical diffuser. The test results were 
corrected to account for diffuser losses between the 
measuring section and the mixing tube exit. These 
corrections may be inaccurate. 

• to determine the measured value of AP, = P -P , the value 

t 111 o so 

of P gQ had to be estimated by subtracting suction duct and 
nozzle cluster losses from the atmospheric pressure. These 
losses cannot be measured under the conditions which exist 
during jet pump operation so the corrections may be inaccurate. 
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the adjustment of the supersonic primary flow to match 
the secondary flow static pressure was assumedtobe isen- 
tropic in the analysis. When the primary flow at the noz- 
zle exit is appreciably over- or under-expanded, the flow 
is non-isentropic and stagnation pressure losses occur in 
the primary flow adjustment process. 


The effects of mixing tube wall friction can be estimated by using the 
loss measurements of figure 43. The analytical values of AP^.* were reduced by 
the appropriate loss taken from figure 43 to derive the column in table 7 entitled 
"Corrected for Duct Friction". The corrected analytical values were used to pre- 
pare the right-hand column which shows the percentage error between the corrected 
analytical value and the test measurement of AP^*, i.e. : 




analytical 

corrected 


analytical 

corrected 



experimental 


P 


so 


The table shows that the difference between the corrected analytical value and the 
test value of AP^.* is 10% or less when traverses in two directions were made during 
the test. The difference can be 20% or more if only one traverse was made, par- 
ticularly when the short mixing tube was used or distorted velocity profiles were 
observed. 


The remaining differences between the corrected analytical value and 
the test value are due to incomplete traversing and possible inaccuracies in the 
conical diffuser and suction duct loss corrections. Furthermore, the pressure loss 
measurements shown in figure43were made by drawing air through the mixing tube 
with a blower. The wall friction losses may be different when the jet pump is opera- 
ting because the mixing action changes the velocity profiles considerably. Thus, 
even the mixing tube wall friction corrections may not be exact. 
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5. 2 Jet Pump Static Pressure Rise 




The experimental measurements of the jet pump static pressure rise 

parameter, AP *, fell consistently below the analytical predictions. The measured 
s 

static pressure values at the mixing tube exit were lower than predicted. The rea- 
sons for this discrepancy were as follows: 

• mixing tube wall friction losses were not included in 
the analysis; these losses are not negligible 

• in the case of the short mixing tube, the tests showed 
that mixing was not completed within the tube. Thus, 
the static pressure jilid not reach its mixed-out value, 
the value which the computer program seeks to predict. 

• the distorted velocity profiles existing in the actual mix- 
ing tube lead to lower static pressures than the uniform 
velocity profiles assumed in the computer analysis. 

• to determine the measured value of AP = P - P , the 

s m so 

value of P go was estimated by subtracting suction duct and 
nozzle cluster losses from the atmospheric pressure. 

These loss corrections may not be exact. 

• the primary flow is not isentropic in the accommodation 
region when the primary flow is over- or under-expanded 
at the primary nozzle exit. Isentropic accommodation is 
assumed in the analysis. 

The differences between the experimental measurements and the analy- 
tical predictions are most pronounced at low entrainment ratios when poor mixing 
occurred, when the short mixing tube was used, and particularly with the 7-nozzle 
LD#2 nozzle cluster. The extended mixing tube improved the recovery of static 
and stagnation pressure at low flow rates. At high flow rates (high entrainment 
ratios), the increased frictional losses in the extended mixing tube cancel the po- 
tential gains from more thorough mixing. 


76 



The computer program was modified as described in sections 3. 1. 5 and 
appendix B. 1 to allow inclusion of mixing tube wall friction and suction duct losses 
in the analytical performance predictions. Preliminary values of and K sc j can be 

obtained by loss tests as in figures 41 and 43, or by estimates using duct loss cor- 
relations available in the literature. More accurate values of K^^and K g ^ must 
be obtained empirically from jet pump testing. 

An example of the empirical approach is shown in figure 86, Test 

values of AP * are replotted for the 260 psia operating condition from figure 72. 
s 

A number of alternative values of K MT were assumed and used as input data for 
the computer performance analysis. The value of K g( ^ was assumed to be 0. 006 
for all calculations. The results show that selection of = 0. 055 makes the 
analytical prediction agree closely with the test results. The loss coefficient 
as used here includes not only the effect of mixing tube wall friction, but also the 
effects of the other sources of discrepancies mentioned above. The mixing tube 
loss coefficient measured by drawing air through the mixing tube with a blower 
(figure 43 ) was = 0. 053. 

This good agreement between the value of K MT measured in the blower 
test and the value of deduced from the computer calculations suggests that the 

blower test method may afford a simple and accurate way to determine K MT for 
a new jet pump design. However, the mixing action which occurs when the jet pump 
is operating may cause the mixing tube wall friction characteristics to vary when 
the primary flow pressure and temperature are changed. To determine whether 
such variations of K MT are significant, more of the test results of section 4 could 
be analyzed by the computer to determine the appropriate values of K^ T . These 
"jet- pump- derived" values of K MT could be compared to as measured by the 

blower test method to determine whether serious discrepancies can occur. This 
report includes sufficient data to make such comparisons. 
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Section 6 


CONCLUSIONS 


The conclusions which may be reached as a result of this investigation 
listed below: 

1. Validity of the Analysis: 

The analytical model developed in section 3. 1 for high-entrainment 
compressible-flow jet pumps with constant area mixing tubes is based 
upon the simplifying assumption that the supersonic primary nozzle 
flow adjusts isentropically to match the secondary flow static pres- 
sure. The analytical predictions of jet pump pressure rise and thrust 
augmentation based upon this assumption agree closely with test re- 
sults when mixing tube wall friction effects are taken into account. 

The agreement is good over a very broad range of operating condi- 
tions. 

2. Mixing Tube Wall Friction: 

The stagnation pressure losses which occur in the mixing tube 
due to wall friction are significant in comparison to the stagnation 
pressure rise developed by a high- entrainment jet pump. Measure- 
ment of the wall friction losses by drawing air through the mixing 
tube with a blower may yield values of the friction loss coefficient 
which are accurate enough for design purposes. Further analysis 
of the data in this report is required to check this point. 

3. Jet Pump- Duct System Matching: 

A technique for selecting the optimum design for a jet pump to 
match given operating conditions was presented in section 3.4. For 
the same primary flow conditions, the jet pump geometry to achieve 
maximum entrainment was shown to be different from the geometry 



required to achieve maximum thrust augmentation. In order to match 
a jet pump to its associated duct system to obtain peak performance, 
it is essential that the loss coefficients of the duct components be es- 
timated as accurately as possible. 

4. Optimum Mixing Tube Length: 

The mixing tube length must be selected to balance the increased 
pressure recovery resulting from more complete mixing in a longer 
tube against the increased wall friction losses in the longer tube. The 
results of this investigation suggest that no simple length- to- diameter 
rule is applicable to mixing tube design. The optimum mixing tube 
length is a function of the primary flow conditions and the entrainment 
ratio at the operating point. This entrainment ratio is set by the loss 
characteristics of the jet pump duct system. 

5. Nozzle Cluster Design: 

The original position of the nozzle cluster close to the mixing tube 
inlet led to increased losses and inferior jet pump performance. When 
the nozzle cluster was moved upstream away from the mixing tube inlet, 
its pressure loss and blockage effects were minimized and the perform- 
ance of the jet pump was measurably improved. A special "low-drag" 
nozzle cluster design was no better than the original cluster when both 
were-positioned away from the mixing tube inlet. 

A reduced centerline velocity appeared in most of the velocity pro- 
files measured at the mixing tube exit and in the conical diffuser. The 
cause of the reduced centerline velocity is unknown. 

6. Performance of the Wagner Jet Induced Lift System: 

The loss characteristics of the suction duct and nozzle cluster 
arrangement used in the NAS 2-2518 test program were measured 
during this program. The high suction duct and nozzle cluster losses, 
incomplete mixing in the too-short mixing tube, and non-optimum match- 
ing of the jet pump to the duct system all acted to reduce the entrainment 
ratios achieved with the system substantially below the peak values attain- 
able. 
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7 . 


Mixing Tube Design: 


The original mixing tube geometry, which had a small reduction 
in cross-sectional area along its length, produced a slightly higher 
entrainment ratio than the constant-area mixing tube under similar 
test conditions. Further analytical and experimental work should 
be carried out to determine the performance characteristics of a 
variety of mixing tube shapes. Significant performance improvements 
may be possible if mixing tubes other than the conveniently-analyzed 
constant area and constant pressure designs are used. 



APPENDIX A 

Listing of the Computer Program 


An before a line indicates that it can be removed if there is to be 
no plotting with the EZPLOT subroutine. 
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0001 



DIMENSION PRISE! 11* 15) , THRUS! 11,15) ,ENT( 11 ,15 ) , AR AY ( 1 1 , 15 ) 

0002 



DIMENSION COUNT! 15) ,X ( 15 ) , Y( 1 5) , PT ITLE < 15 ) ,TT ITLE ( 15 ) ,NC <9 ) 

0003 



DIMENSION CHARI! 11,6 ) ,CHAR2< 11,5) 

0004 



DATA NC/38, 63, 16,55,44, 19,24, 52/ 


C 


CLEAR ARRAYS 

0005 



DO 26 K= 1 , 15 

0006 



COUNT ( K)=0 

0007 



X! K)=0 

0008 



Y ( K ) = 0 

0009 



PT ITLE ( K ) =0 

0010 



TT ITLE ( K ) =0 

0011 



DO 25 J=l, 11 

0012 



PRI SE ( J , K ) = 0 

0013 



THRUS! J,K)=0 

0014 



ENT! J,K)=0 

0015 



ARAY! J, K ) =0 

0016 


25 

CONTINUE 

0017 


26 

CONT INUE 


C 

CALL IDFRMV REMOVED FROM PROGRAM HERE 

0018 

Line 1 


READ(5,350)NTP, (PTITLE! I ) , 1=1, 15) 

0019 

Line 2 


READ! 5,450)NTT, (TTITLE! I) ,1=1,15) 

0020 

Line 3 


RE AD ( 5 , 5 2 ) N 

0021 



G=1.4 

0022 



G0=3 2.2 

0023 



R = 53 . 3 5 

0024 



CONV =778 . 16 

0025 



K = 0 


C 


BEGIN OUTER LOOP, EACH LOOP USES A NEW SET OF INITIAL CONDITIONS 

0026 


1 

WRITE(6, 100) 

0027 



DO 761 J=l, 11 

0028 



DO 760 M= 1 , 5 

0029 



CHARI! J,M)=0 

0030 



CHAR 2 ( J , M ) = 0 

0031 


760 

CONT INUE 

0032 


761 

CONTINUE 

0033 



CHARI! 11,6 )=0 

0034 



K = K+1 


C 


READ DIMENSIONLESS INITIAL CONDITIONS 

0035 

Line 4 


READ! 5,50,END=300)PBAR,TBAR,ABAR1, ABAR2 


C 


READ DIMENSIONAL INITIAL CONDITIONS 

0036 

Line 5 


READ(5,51)PSOI ,TSO,AM 

0037 

Line 6 


READ(5,455)FDUCT,FTUBE 

0038 



WRITE(6, 101)PBAR,TBAR 

0039 



WRITE (6, 102) ABAR1, ABAR2 

0040 



WRITE (6, 103) PSOI ,TSO, AM 

0041 



WRITE! 6,456 IFDUCT, FT UBE 


C 


CALCULATE OTHER DIMENSIONAL VALUES FROM DIMENSIONLESS VALUES 

0042 



PPO =P BAR*P SO I 
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0043 

0044 

0045 


0046 

0047 

0048 

0049 

0050 

0051 

0052 


0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 
0061 

0062 

0063 

0064 

0065 

0066 

0067 

0068 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 

0079 

0080 
0081 

0082 

0083 

0084 

0085 

0086 


C 

C 

Line 7 


C 


Line 8 
Line 9 C 


Line 10 


C 


TP0=TBAR*TS0 

ATH=ABAR1*AM 

AP=ABAR2*AM 

BEGIN INNER LOOP, EACH TIME THROUGH GIVES A SET OF SOLUTIONS 
FOR A DIFFERENT VALUE OF ENTRAINMENT RATIO. 

DO 13 J=l,ll 
ENTR=7.0+J*3.0 
WR I TE ( 6 , 105 ) J» ENTR 

CALCULATE PRIMARY AND SECONDARY MASS FLOW RATES 
WP=G*GO*( 1/R)* <2.0/( G+l) )** ( (G+l) /( G-l) ) 

WP=SQRT(WP )* 144. 0*60. 0* PPO* ATH/SQRT < TPO ) 

WS=ENTR*WP 

PSO=PSOI-FDUCT*WS*WS*R*TSO 
1 / ( 3 600. 0* 2. 0*G 0*P SOI * AM* AM* 144. 0- 144 . 0 ) 

ITERATE TO FIND SECONDARY MACH NUMBER 
820 1=0 

GUE SS=0. 2 

Z =W S* SO R T ( T SO ) /(60. 0*144. 0*PS0* ( AM-AP) *SORT (G-GO/R ) ) 

2 PARAM=1+(G — 1) /2. 0— GUE SS^GUE SS 
ZCALC=GUESS/PARAM**< ( G+ 1 ) / < 2. 0* ( G- 1 ) > ) 

IF (GUESS-1. 0) 14, 14, 15 

14 I F ( 1-100) 3,3,5 

3 IF(ABS(ZCALC-Z)-.0005)6,6,4 

4 DERIV=PARAM**< -(G+l) / ( 2 • 0*{ G- 1 ) ) > -( G+ 1 ) /2 .0*GUE SS*GUESS 
1-PAR AM’S* ( ( 1-3. 0*G) /( 2.0* (G-l) ) ) 

GUE SS = GUESS-( ZCALC-Z ) /DERI V 
1 = 1+1 
GO TO 2 

5 WR IT E ( 6 , 66 ) 

GO TO 405 

15 WRITE (6, 115) 

GO TO IB 

6 CONTINUE 
SMO K=GUE SS 

PS = PSO/< 1+ ( G- 1 ) /2.0*SM0K*SM0K)*=MG/(G-1) ) 

TS = T SO / ( l+( G-l ) /2. 0- SMOK# SMOK ) 

VS = 2.0*G*G0*R*< TS0-TS)/(G-1> 

VS= SORT (VS) 

PSOR=PSOI— F DUCT* VS* V S- P S/ ( 2.0*G0*R*TS ) 

IF ( ABS< ( PSOR-PSO ) /PSOR ) -. 0005) 800,800,810 
810 PSO = P SOR 
GO TO 820 
800 PSO=PSOR 

WRITE(6»457)PS0 
W R I T E ( 6 ,600 ) I 

MAKE CORRECTION FOR UNDER OR OVER EXPANSION 
1=0 

16 PMOK=( 1+ ( G — 1 ) /2. 0* SMOK-SMOK ) * ( PPO/PSO)**( (G-l)/G)-1.0 
PMO K= SQRT ( 2. 0*PM0K/ ( G-l ) ) 

AM2=W S* SOR T ( TSO )*< l + ( G-l ) /2. 0*SM0K*SM0K) ** ( (G + l )/(2.0=MG-l) ) ) 

1 / ( P SO* SMO K ) 

AM2=AM2+WP*SQRT< TPO )* ( 1+ ( G- 1 ) /2. 0* PMOK*PMOK ) ** ( < G + l ) / ( 2 . 0* (G- 1 ) > > 
l/( PPO*PMOK) 
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Ik 


008 7 
0088 
00 8 9 

0090 

0091 

0092 

0093 

0094 

0095 

0096 

0097 

C 

0098 

0099 

0100 
Old 
0102 

0103 

0104 

0105 

0106 

0107 

Line 11 C 

0108 

0109 

0110 
0111 
0112 

0113 

0114 

0115 

0116 

0117 

0118 

0119 

0120 
0121 
0122 

0123 

0124 

0125 

0126 

0127 

0128 

0129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 


AM 2= AM 2* SORT ( R / ( G*G0) ) / ( 1 44. 0* 6 0 . 0 ) 

I F ( ABS ( AM2-AM) -.00005) 19, 19, 17 

17 SMOK=SMDK+ ( ( AM2- AM) /ABS( AM2-AM) ) * F X P ( 5 .2-SMOK-SMOK ) *5 . E - 5 
I F ( 1-200) 18,5,5 

18 CONTINUE 
1 = 1 + 1 

GO TO 16 

19 CONTINUE 
WRITE(6,601) I 

AP= ( 1+ ( G- 1 ) / 2. 0*PM0K*PM0K)** ( ( G+ 1 ) / ( 2 .0* (G-l > ) ) 
AP=AP*WP*SQRT( TPO ) /( 144.0*60. 0*PPO*PM0K*SQRT (G*G0/R ) ) 
CALCULATE CONDITIONS AT END OF ACCOMODATION REGION 
PS2=PS0 / ( 1 + ( G-l ) /2.0*SM0K*SM0K)**(G/(G-1) ) 

T S2 =T SO / ( l+IG-l) /2.0*SM0K*SM0K) 

VS2= 2. 0*G*GO*R*< TS0-TS2) /( G-l ) 

VS2=S0RT (VS2) 

PP2=PS2 

TP2=TPO*(PP2/PPO)**< (G-l) /G) 

VP2=2.0*G*G0*R*< TP0-TP2) /(G-l) 

VP2=S0RT (VP2) 

C2 = S0RT< G*G0*R*TP2) 

PMO K= VP 2/C2 

BEGIN ITERATION TO FIND OUTLET CONDITIONS 
W M = W P * ( 1 + ENTR ) 

CP2=. 24914 
CS2=. 24914 
1 = 0 

HP2 = CP 2* TP 2 + VP 2* VP 2/ ( 2 . 0*G0*C0N V ) 
HS2=CS2*TS2+VS2*VS2/(2.0*G0*C0NV) 

I NT = 0 

PMG = P S2+ . 4 

7 VM=VP2+ENTR*VS2+60. 0*144. 0*G0*AM*( PS2-PMG)/WP 
VM=VM*WP/WM 

TM=60. 0* 144. 0*PMG*VM* AM/ ( R*WM) 

CM=. 24914 

PM=(HP2+ENTR*HS2)*R*WP/( 60. 0* 1 44 .0* VM*AM*CM > 
PM=PM-R*WM*VM/< 2. 0*144. 0*60. 0*GO*CONV*CM*AM ) 

I F ( I -800 ) 8,8,15 

8 I F ( ABS( PMG-PM) -.01) 11, 11,9 

9 I F ( PM .GT. PMG) GO TO 20 

I F ( INT .EO. 1 ) GO TO 11 

PMG=PMG+ ( PMG-PM) /PMG 

1 = 1 + 1 
GO TO 7 

20 INT= 1 

PMG= PMG-. 001 
1 = 1 + 1 
GO TO 7 

10 WRITE ( 6,66 ) 

GO TO 405 

11 CONTINUE 
WRITE(6,602) I 
PM=PMG 
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0138 

Line 12 

0139 

0140 

0141 

0142 

0143 

0144 

0145 

0146 

0147 

0148 

0149 

0150 

0151 

0152 

0153 

0154 

0155 

0156 

0157 

Line 13 

0158 

0159 

0160 
0161 
0162 

0163 

0164 

0165 

0166 

0167 

0168 

0169 

0170 

0171 

0172 

0173 

0174 

0175 

0176 

0177 

0178 

0179 

0180 
0181 
0182 

0183 

0184 

0185 

0186 


EM0K=VM/SQRT(G*G0*R*TM) 

C CALCULATE THE EFFECT OF WALL FRICTION 

GAM 1=1.0+ (G-l) /2.0* EMOK* EMOK 

ALPH= 144.0* PM* (GAM1** ( G/IG-l) ) -FTUBE* EM0K*EM0K*G/2 . 0 ) 
BETA=ALPH/SQRT( 144. 0* PM* 144. 0*PM*EM0K*EM0K*GAM1 ) 

1 = 1 

GMOK=EMOK 

700 G AM 2= 1 • 0+ ( G-l ) /2. 0*GM0K*GM0K 

A=< l/GMOK)*GAM2**( ( G+ 1 ) / ( 2. 0* ( G- 1 ) ) ) 
DA=(G+1)/2.0*GAM2**( ( 3-G ) / ( 2. 0* ( G-l ) ) ) 

DA=D A-< 1 / ( GMOK*GMOK ) ) *GAM2** ( ( G+ 1 ) /( 2.0* (G-l ) )) 

IF ( I .GT. 200)60 TO 10 
IF (GMOK-1. 0)701,701,15 

701 I F ( ABS( A-BETA) -.0001 >703,703,702 

702 1=1+1 

GMOK=GMOK-( A-BETA) /DA 
GO TO 700 

703 CONTINUE 
WRITE(6»604) I 

P2 = ( 1/144. 0)*ALPH/GAM2** (G/( G-l) ) 

TM=TM* (P2/PM)**2. 0*( GMOK/EMOK ) **2 .0 
C CALCULATE OUTLET PARAMETERS 

VM=GMOK*SORT(G*GO*R*TM) 

PM = P2 
EMO K=GMOK 

PMTOT = 144. 0* P M* ( 1 . 0+ ( G— 1 ) / 2 • 0* EMOK* EMOK ) ** (G/ (G-l ) ) 

D EL P = PM TOT— 144. 0*PS0 
DDEL P = DE LP / ( 144.0*PSO) 

DST AT = 144 . 0* ( P M— P SO) 

DDSTA=(PM— P SO ) /P SO 
AUG=WM*VM / ( W P* VP2 ) 

ENER= 144. 0*PM*VM*VM/( 2.0*G0*R*TM) 

DENER=ENER/( 144.0*PSO) 

NOEL P = DE LP 

IF ( NDELP .LE. 0) GO TO 1 
C STORE SOLUTIONS FOR TABLE PRESENTATION 

SUM = J 

CHAR1( J, 1 ) =ENTR 
CHAR1( J,2)=DDELP 
CHAR1( J, 3) =DDSTA 
CHARlf J , 4 ) = E MO K 
CHARK J, 5 ) =DENER 
CHAR1( J,6)=AUG 
C H AR 2 ( J , 1) =ENTR 
CHAR2 ( J , 2 ) = PM 
CHAR2( J, 3) =DELP 
CHAR2 ( J ,4 ) = DST AT 
CHAR2I J,5 )=ENER 

C STORE SOLUTIONS IN ARRAYS TO BE PLOTTED LATER 

PRISEt J,K)=DDELP 
THRUSt J , K ) = AUG 
ENT ( J, K) =ENTR 
COUNT ( K ) = J 


0187 

0188 

0189 

0190 

0191 

0192 

0193 

0194 

0195 

0196 

0197 

0198 

0199 Line 14 

0200 
0201 
0202 

0203 

0204 

0205 

0206 

0207 

0208 

0209 

0210 
0211 
0212 

0213 

0214 

0215 

0216 

0217 

0218 

0219 

0220 Line 15 

0221 
0222 

0223 

0224 

0225 

0226 

0227 

0228 

0229 

0230 

0231 

0232 

0233 


WR IT E ( 6, 106) 

WRI TE ( 6 , 107 ) 

WRITE(6»108)PP2,PS2,PM 
WRITE(6»109)TP2»TS2* TM 
WR ITE ( 6, 110)VP2»VS2*VM 
WRITE (6, 11DPM0K, SMOKtEMOK 
WRITE (6* 112)WP,WS,WM 
WRI TE ( 6 * 1 16 ) AP 
WRITE(6*113)DELP 
WRITE(6»120)DDELP 
WRITE(6*603)DDSTA 
WR ITE ( 6 * 114) AUG 
C END OF INNER LOOP 

13 CONTINUE 

C PRINT OUT TABLE OF RESULTS 

765 WRITE(6»707) 

WRITE<6,710) 

WRITE (6,75 1)PBAR,TBAR, ABAR1 
WRITE( 6,707 ) 

WR ITE ( 6 ,720) 

WRITE(6,730) 

M= SUM 

DO 725 J=1,M 

WRITEI6,704)CHAR1( J , 1 ) , CHAR 1 ( J , 2 ) , CHAR 1 ( J , 3 ) , CHARI ( J , 4 ) , CH AR 1 ( J , 5 ) 
1,CHAR1( J,6) 

725 CONTINUE 
WRITE(6,707) 

WRITE!6,750) 

WRITE(6,751)PBAR,TBAR,ABAR1 
WRITE(6,752)PS0, TSO,AM 
WRITE(6,707) 

WRITE (6,753) 

WRITE(6,754) 

DO 726 J= 1 , M 

WRITE (6, 755 >CHAR2( J, 1 ) , CH AR2 ( J, 2 ) , CHAR2 ( J , 3 ) ,CHAR2 ( J ,4 ) ,CH AR2 { J , 5 ) 

726 CONTINUE 

C END OF OUTER LOOP 

I F ( K .GE. N) GO TO 300 
GO TO 1 

300 CONTINUE 

C FIND LARGEST VALUES IN ARRAYS PRISE AND THRUS 

PMAX=TOP ( PRI SE) 

TMAX=TOP ( THRUS ) 

C ENTER PLOTTING SECTION OF THE PROGRAM 

K= 1 
NA=17 
N0= 27 
NF = 1 

GO TO 303 

301 NTP=0 
N A= 0 
N0= 0 
NF = 2 
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0234 

0235 

0236 

0237 

0238 

0239 

0240 

0241 

0242 

0243 

0244 

0245 

0246 

0247 

0248 

0249 

0250 

0251 

0252 

0253 

0254 

0255 

0256 

0257 

0258 

0259 

0260 
0261 
0262 

0263 

0264 

0265 

0266 

0267 

0268 

0269 

0270 

0271 

0272 

0273 

0274 

0275 

0276 

0277 


302 K=K+ 1 

303 M=C0UNT(K> 

DO 305 J=l,ll 
X< J)=0 

Y ( J ) = 0 
305 CONTINUE 

DO 304 J= 1 » M 
X< J)=ENT< J,K) 

Y( J)=PRISE< J,K) 

304 CONTINUE 
NP=COUNT ( K) 

C CALL EZPLOT REMOVED FROM DECK HERE 
IF ( K .GE. N ) GO TO 400 
IF(K-2)301»302»302 

400 K= 1 

N A= 17 
N0 = 28 
NF = 1 

GO TO 403 

401 NTT=0 
NA=0 
N0 = 0 
NF = 2 

402 K=K+1 

403 M=COUNT(K) 

DO 406 J-1,11 
X( J)=0 

Y ( J ) =0 

406 CONTINUE 

DO 404 J= 1 * M 
X( J)=ENT( J,K) 

Y( J)=THRUS( J,K) 

404 CONTINUE 
LAST=1 

I F ( K .GE. N ) LAST=2 

C CALL EZPLOT REMOVED FROM DECK HERE 
I F ( K .GE. N) GO TO 405 
I F t K-2 ) 401 » 402* 402 

C CALL PLTND REMOVED FROM DECK HERE 

405 CONTINUE 

50 FORMAT < 4F 10.5) 

51 FORMATt 3F10.5 ) 

52 FORM AT ( I 10 ) 

100 FORMAT ( / /25Xf65H*=*’S”S <, S ! THE FOLLOWING CASES WILL USE THESE INITIAL C 
10ND I T IONS *****) 

101 FORMAT ( /2X, 15HPRESSURE R AT I 0= , E 10 .4 , 1 OX , 18HTEM PER ATURE RATIO=, 
1E10.4) 

102 FORMAT ( /2X , 12H ( ATH*CW/AM) =, E10. 4, 1 3X , 6HAP / AM= , E 10.4 ) 

103 FORMAT (/2X»19HSEC0NDARY PR ESSURE= , E 10 .4 , 3X , 22HSECOND ARY TEMPERATUR 
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0278 


0279 

0280 
0281 
0282 

0283 

0284 

0285 

0286 

0287 

0288 

0289 

0290 

0291 

0292 

0293 

0294 

0295 

0296 

0297 

0298 

0299 

0300 


0301 


0302 

0303 

0304 

0305 


0306 


0307 

0308 

0309 

0310 

0311 


&- 


1E=,E10.4,3X,17HMIXING TUBE AREA=,E10.4) 

105 FORMAT ( /40X, 10H*#** CASE , I 2 , 2X, 1 8HENTRA INMENT R AT 10= ,F 10.5, 

15 H ***=*> 

66 FORMAT ( /2X» 19HT00 MANY ITERATIONS) 

106 FORMAT ( /39X, 19H ACCOMODATION REG ION , 10X , 1H* , 10X ,6H0UTLET ) 

107 FORMAT! 34 X,7HPRIM ARY, 1 IX , 9HS ECONDARY ,7X , 1H* ) 

108 FORMAT! /2X , 13HPRE SSURE , PS I A, 15X , E 10 .4 , 10X , E 10 .4 , 8X , 1H* ,6X ,E10.4> 

109 FORMAT! /2X, 11HTEMP ERA TURE, 1 7X , E 10 .4 , 10X , E 10 .4 , 8X , 1H* , 6X , E 10.4 ) 

110 FORMAT ( /2X * 15HVEL0CITY, FT/SEC, 1 3X , E 10 .4 , 1 OX , E 10 .4 ,8X , 1H* ,6X ,E10.4) 

111 FORMAT! /2X, 11HMACH NUMBER, 17X,E10.4,10X,E10.4,8X,1H*,6X,E10.4) 

112 FORMAT! /2X,24HMASS FLOW R ATE , POUND/M I N ,4X , E 10 .4 , 1 OX , E 10 .4 , 8X , 1H* , 
16X,E 10.4) 

113 FORMAT ( /5X, 20HT0TAL PRESSURE R I SE= , E 10 ,4 , 2X , 14HP0UND /SQ . FOOT) 

114 FORMAT! /5X, 15HM0MENTUM R AT I 0= , E 10 .4 ) 

115 FORMAT! /5X,47HTHIS VALUE OF ENTRAINMENT RATIO IS INACCESSIBLE) 

116 FORMAT ( /5X,43HTHE PRIMARY STREAM AREA AFTER ACCOMOD AT ION = , E 10. 4, 

1 IX , 8H SQ . FEET) 

120 FORMAT! /5X , 28HDI MENSI ONL E SS PRESSURE RISE=,E10.4) 

350 FORMAT! I 15, 15A4) 

450 FORMAT ( I 15, 15A4) 

600 FORMAT ( /5X, 23HMACH NUMBER I TERAT I ONS= ,13) 

601 FORMAT! /5X , 24HP R I MARY AREA I TE R AT I ONS= , I 3 ) 

602 FORMAT! /5X, 20HPRESSURE I T ER AT I ONS= , I 3 ) 

603 FORM AT ( /5X, 37HDIMENSI ONL ESS STATIC PRESSURE CHANGE = ,E 10.4 ) 

707 FORMAT < /IX, 116H* si:************************ 

710 FORMAT! 15X,54HDIMENS IONLESS SOLUTION USING THESE INITIAL CONDITION 
IS, ) 

720 FORMAT! /3X, 11HENTRAI NMENT, 4X, 1H* , 2X , 1 9HD I MENS I ONLESS TOT AL , 3X , 1H* , 
12X.20HD I MENS IONLESS STATIC, 2X, 1H*,4X,4HMACH,6X,1H*,2X , 13HD IMENS I ON 
2LESS»4X,1H*,3X,8HM0M ENTUM ) 

73 0 FORMAT! 3X,5HR AT 10, 1 OX , 1H* , 2X , 1 3HPRE SSURE R I S E ,9X , 1H* , 2X , 15HPRE SSUR 
IE CHANGE, 7X, 1H*, 4X , 6HNUMBER , 4X , 1H* ,2X , 14HK I NET I C ENERGY ,3X , 1H*, 3X , 
25HRATIO) 

704 FORMAT! /8X,F4. 1 , 6X , 1H* , 4X , F 10. 5 , 10X , 1H* , 5X , F 10 . 5 ,9X , 1H* ,4X , F5 .3 , 

15 X, 1H*,3X,F7.4,9X, 1H*,3X,F7.3) 

750 FORMAT! /5X , 5 2HDI MENSI ONAL SOLUTION USING THESE INITIAL CONDITIONS, 
1) 

751 FORMAT (5X, 15HPRESSURE RATIO= , F 5 . 2 , 3X , 1 8HTEMPERATURE R AT 1 0 = , F5 . 2 , 3X 
1,12H< ATH*CW/AM)=,F8.6> 

752 FORMAT! /5X, 19HSEC0NDARY PR ESSURE= , F5 .2 , 1 X , 3HPS I ,3X ,22HSECONDAR Y TE 
1MPERATURE = ,F7. 2, IX , 14HDEGREE R ANK I NE , 3X , 17HM I X I NG TUBE AREA=,F8.6, 
2 IX ,7HSQ. FT.) 

753 FORMAT! /3X, 11HENTRAI NMENT, 2X, 1H*,2X ,11HEX IT STATIC, 5X,1H#,2X,14H TO 
IT AL PRESSURE, 2X, 1H*,2X, 15HSTATIC PRE SSURE , 2X , 1H* , 2X , 14HK I NET I C ENE 
2RGY ) 

754 FORMAT!3X,5HRATIO,8X, 1H~ » 2X , 14HPRESSURE 1 PS I ) , 2X , 1H* , 2X , 10HR I SE !P 
1SF) ,6X, 1H#,2X, 12HCHANGE (PSF),5X,1H*,5X,5H(PSF)) 

755 FORMAT! /8X,F4. 1,4X, 1H* , 4X , F5 . 2 , 9X , 1H* ,3X ,F7. 2 ,8X , 1H*, 3X , F9 .4 ,7X , 

1 1H* ,5X,F7.2) 

604 FORMAT! /5X, 20HFRICTI0N I T ER AT I ONS= , 1 3 ) 

455 F0RMAT(2F10.5) 

456 FORMAT! /2X,23HIMLET LOSS COE FF IC I ENT= , E 1 0.4 , 3X , 24HM IX ING LOSS COEF 
1FICIENT=,E10.4) 
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0312 

0313 
0001 
0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 


457 F0 RM AT (/5X»39H SECONDARY TOTAL PRESSURE BEFORE M IX ING= , E 10.4 ) 
END 

FUNCTION TOP ( AR AY > 

DIMENSION AR AY ( 1 1 t 15 ) 

TO P= ARAY ( It 1) 

DO 501 L= 1 t 11 
DO 500 M= 1» 15 

IF ( TOP .LT. ARAY(LtM)) TOP= ARAY < L t M ) 

500 CONTINUE 

501 CONTINUE 
RETURN 
END 
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APPENDIX B 


B. 1 Discussion of Computer Program by Blocks 

Block 1 - Declare and Clear Arrays, Set Values of Constants: 

To make sure that the storage arrays used in the program are all empty, 
zeroes are placed in every location by this block of the program. 

Numerical values of G, GO, R, and CONV are also defined in this sec- 
tion of the program. The numerical values are not read in as input; they are de- 
fined within the program. To change them, the appropriate cards in the deck must 
be changed. 

Block 2 - Read Initial Conditions and Initialize Parameters: 

In the segment of the program between line 1 and 6, the initial conditions 
are read in and prepared for further calculations. These initial conditions include 
data required by the plotting routine. 

Titles for plots of dimensionless pressure rise and momentum ratio as 
functions of entrainment ratio (PTITLE, TTITLE) and the number of cases to be 
solved (N) are introduced in lines 1, 2 and 3 respectively. The number preceding 
each title is the number of characters and spaces in the title. Dimensionless initial 
conditions (PBAR, TBAR, ABAR1 and ABAR2) are read in on line 4. Although 
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ABAR2 has not been used in the formulation, it is used as an initial guess for an 
iteration loop as described in more detail in the discussion of block 5. Since a value 
of k has been specified in block 1, the initial conditions for a dimensionless solution 
are complete at this point. 

Values of PSOI, TSO and AM are read in on line 5, and values of fric- 
tion coefficients (FDUCT, FTUBE) are introduced at line 6. This completes the in- 
formation required to begin calculations. Specific values of PPO, TOI, ATH and 
AM are obtained by multiplying the dimensional initial conditions by the non- 
dimensional initial conditions. 

Block 3 - Calculation of Suction Duct Pressure Loss: 


The performance of a jet pump is dependent upon the stagnation pres- 
sure of the secondary flow at the exit of the primary nozzles, PSO. It is difficult 
to measure this pressure in an experimental jet pump. Instead, the suction duct 
inlet stagnation pressure can be measured and the loss between the inlet section 
and the primary nozzle exit section can be accounted for by the following equation: 


A P. 


= K 


P V 
s s 


suction duct 


sd 2 g 


(97) 


The loss coefficient K , can be evaluated by drawing air through the 

so 2 

suction duct with a blower and plotting (AP.) vs. p o V /2g on a 

t suction duct s b ° 

Cartesian graph. The slope of the resulting curve is equal to K gc j- 


The correction for suction duct pressure loss begins on line 8 of the 


program. The equation for this correction is given below. 

P V 2 

P so = P soi - K sd 


(98) 


where P so ^ is the secondary flow stagnation pressure at the suction duct inlet. 
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Blocks 4 and 10 - Specification of Entrainment Ratio: 


Line 7 initiates a DO loop which encloses all statements down to line 14 
(statement number 13). The entrainment ratio (ENTR) is initially set equal to 10. 0 
and the loop is repeated, increasing ENTR by 3. 0 each time through, until a value 
of 40. 0 has been reached. The integer J is used to count the number of loops com- 
pleted. 

Block 5 - Solution of Equations for the Accommodation Region: 

Block 4 extends between line 9 and line 11. Primary and secondary mass 
flow rates (WP, WS) are calculated using equations ( 1) and (2 ). The value of the 
secondary flow Mach number (SMOK) is obtained by solving equation (57). Equation 
(57) is solved using Newton's method. A series of values of are tried until one 
is found which makes the left side of equation (57) sufficiently close in value to the 
right side. The numerical criterion for acceptable convergence is agreement with- 
in 0. 0005, which gives a secondary flow rate accurate to within 0. 5 lbm/min. 

The primary mass flow rate and the entrainment ratio are specified in 
this program before the secondary Mach number is evaluated. There exists a maxi- 
mum possible secondary mass flow rate corresponding to a secondary Mach number 
of 1. 0. Thus, for given initial conditions, the jet pump will have a maximum per- 
missible value of entrainment ratio, which can be calculated from equation (64). 

The statement on line 10 was inserted to recognize trial values of Mg which are 
greater than 1. 0. If the program tries to obtain a solution for an entrainment ratio 
greater than the maximum possible for the system, a message "This value of en- 
trainment ratio is inaccessible" is printed and the program proceeds to block 9 (line 
number 15) to look for another set of initial conditions. 

When the secondary Mach number has been determined, the static pres- 
sure at the end of the accommodation region (P 9 = P 0 ) is calculated using P and 

pz sz so 

M^ in equation ( 3). Then equation ( 3 ) is used with P^ o and P^ to calculate M^- 
This procedure neglects the effect of primary stream expansion or contraction upon 
the secondary stream area. Equation ( 6 ) is used to correct for primary flow area 
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changes as follows. The calculated values of M^ 2 and M g2 are substituted into the 
left side of equation (6 ) and the result is compared with the actual mixing tube area. 
If the values are not sufficiently close, the value of M g2 is changed slightly, a new 
value of Mpg is determined from equation (3 ), and the test to see if equation (6 ) is 
satisfied is repeated. When values of M g2 and M^ 2 which satisfy equation (6 ) are 
determined, equations (7), (8 ), (9 ), (10), and (11) are employed to calculate the 
pressure, temperature and velocity for each stream at the end of the accommodation 
region. 


The solution of equation (57) for M g2 requires a value of A^, the area of 
the primary nozzle flow, which is introduced as an initial condition through the di- 
mensionless area ratio ABAR2. However, this value of M g2 is used only temporarily; 
it is eventually modified as the effect of primary area changes is considered. There- 
fore, even though the analytical formulation did not make use of ABAR2, the program 
is made more efficient by using it to obtain a preliminary value of M . 


Block 6 - Solution of Equations Representing the Mixing Region: 


The segment of the program which solves the mixing region equation ex- 
tends from line 11 to line 12. The solution is obtained by an iteration technique in 
which the value of one variable is assumed and the equations are solved to obtain a 
calculated value of that variable. When the trial value is sufficiently close to the 
calculated value, an acceptable solution has been determined. 


Values of stagnation enthalpy (HP2, HS2) are calculated for the primary 
and secondary streams using their properties at the end of the accommodation re- 
gion. A trial value of PM (PMG = PSO+ 0.4) is chosen and equations (19), (18), (17) 
and (20) are solved sequentially to determine a calculated value of PM. If the trial 
and calculated values satisfy the criterion below, the trial value is accepted as the 
solution. 

j PMG - PM) < 0.01 


If this convergence criterion is not satisfied, a new trial value is chosen and the pro- 
cedure is repeated. The convergence criterion corresponds to an error of less than 0.1%. 
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The integer I keeps track of the number of iterations performed in order 
to limit their number to 400. Convergence is generally obtained in less than 100 it- 
erations. In some cases, the trial value will oscillate about the calculated value 
without quite converging. When this occurs, the trial value is within a region very 
close to convergence. Therefore, the logical variable INT is used to detect this con- 
dition, stop the iteration, and accept the current trial value as a satisfactory solu- 
tion. 


In the course of the development of a satisfactory value of PM, corres- 
ponding values of VM and TM are determined. These three values completely speci- 
fy the conditions at the end of the mixing region. If no frictional effects are included, 
all the desired jet pump performance parameters such as stagnation pressure rise, 
static pressure change and momentum ratio can be determined from PM, VM, and 
TM. 


For the temperature ranges of interest, the variation of specific heat 
with temperature has been found to have a negligible effect on the results (figure5). This 
computerprogram, therefore, has beenwritten to treat specific heatand k as constants. 

As a result, equation (17) has not been used in the program. If it is desired to in- 
clude a variable specific heat, an equation such as (17) can be added to the existing 
program without altering its basic structure. 

Block 7 - Solution of Equations Representing Mixing Region Friction Effects: 

Block 7, from lines 12 to 13, corrects the values of M , P m , T^ and 
V m for frictional effects. Equation (38) is used to accomplish this. For a specified 
value of mixing duct friction coefficient, K.^,, both a and ft can be determined 
from the results of block 6. Equation (38) is then solved using Newton's method to 
give a new value of Mach number at location m. A new value of pressure is obtained 
from equation (34). Corrected values of temperature and velocity are then determined 
using equations (39) and (33). 

Block 8 - Print Solutions for One Value of Entrainment Ratio: 

When a solution for a particular value of entrainment ratio is obtained, 
all variables are printed together with appropriate titles. Data internal to the program 
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is also printed, e .g. , the number of iterations required for convergence in each sec- 
tion in which an iterative method is used. When a new set of initial conditions are 
used, this is indicated by an appropriate statement and a list of the new initial condi- 
tions. 

Block 9 - Store Solutions for Plots and Tabular Presentation: 

The results which are printed as each solution is obtained are also sum- 
marized in tabular form. To allow the tabular form to be printed, arrays are filled 
with the numbers from each solution. 

Arrays of the same type are used to store solutions for the plotting routine. 
These arrays are larger because the plotter is used for several sets of initial con- 
ditions whereas the tabular results are printed for each individual set of initial con- 
ditions. 

Block 11 - Print Solutions in Tabular Form: 

After a set of solutions for one set of initial conditions has been obtained 
and stored in the arrays of block 9, the results are printed in a table. Values of the 
desired performance parameters are printed vs. entrainment ratio. Two tables are 
printed; one presents dimensionless variables and the other presents dimensional 
variables. The dimensionless solutions are independent of the dimensional initial 
conditions unless frictional effects have been included in the solution. 

Block 13 - Plot Dimensionless Pressure Rise and Momentum Ratio: 

Values which have been stored in the arrays PRISE and THRUS in block 9 
are plotted as functions of entrainment ratio. Two plots are obtained; dimensionless 
pressure rise and momentum ratio vs. entrainment ratio. If it is desired to plot 
variables other than these, this may be done by storing the desired variables in the 
PRISE and THRUS arrays in place of the dimensionless pressure rise and momen- 
tum ratio values. 

The plotting system is too complex to discuss in detail here. Unless the 
user has access to the EZPLOT routine used in this program, it is unlikely that he 
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will be able to use the plotting section directly. However, the storage arrays devel- 
oped within the program should be useful for providing data to any other plotting de- 
vice which may be available. 

B. 2 Using the Computer Program 

This section discusses the various options which may be exercised when 
using the program and the tasks which must be performed to set them up. 

The Complete Program: 

The complete program includes suction duct and mixing tube frictional 
effects and plotting of the results. Values of PBAR, TBAR, ABAR1, ABAR2, PSOI, 
TSO, and AM must be provided as input data together with empirically-determined 
values of FDUCT and FTUBE. A value of N equal to the number of cases to be 
solved must be included to maintain control within the plotting block. It is best to 
vary only one parameter such as A in a single set of solutions. The plot titles 
(which are provided as input data) can then denote this parameter on the plots as il- 
lustrated in the sample solutions provided in section 3. 3. 

If variables other than DDELP and AUG are to be plotted, this may be 
done by storing them in the PRISE and THRUS arrays. 


Omission of Frictional Effects: 

Frictional effects may be omitted by inputting zero values for FDUCT 
and/or FTUBE, depending on which frictional effect is to be eliminated. 

Omission of the Plotting Section: 

The plotting function may be temporarily omitted by punching a charac- 
ter in the first column of the "CALL EZPLOT" cards. The necessary input data for 
plotting must still be read-in or else the input format will not function correctly. 


96 



If the plotting section is to be removed permanently, it is best to re- 
move all cards associated with the plotting function. This includes READ state- 
ments, logic for data manipulation, storage arrays, and the function TOP (ARRAY). 
The cards which may be removed are noted by an asterisk on the program listing 
of Appendix A. 
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APPENDIX C 


Typical Computer Solutions 



***** THF Foil OWING C ASF S WILL 'IS F TH*SF INITIAL CONDITIONS 


PRESSURE «ATins0.3000F 0? tcmPFPATIJRF RATI 0*0.' *QOE <11 

< ATH*CW/AM)*0.’ OOOF-CP AP/ AM=0.R000E-0? 

S cr ON^ARY PRFSSURF-O.UfOF 0? SECONDARY TE Mpf R A TI|R E *0. R100 p 01 HIDING TUBE APF A-0. R77AF-0 1 

**♦* CAST » 6NTR A I NMFNT PATIO* 10.00000 **** 

MACH NUMBER ITERATIONS- 2 
PRIMARY AREA ITFR AT !ONS= 10 
PRFSSURF ITERATIONS. 7? ' 



ACCOMODATION RFGlON 


* outlet 


PRIMARY 

SECONDARY 


PR F S SURE ♦ PS l A 

0.1439E 02 

0.1A39E 

02 

• 0.1493E 0? 

TEMPERATURE 

0.2996E 03 

0.M2 78F 

01 

• 0.549RE 03 

VEL0C1TY.FT/SEC 

0.2441E 04 

0.1624F 

03 

4 O.1T06F 03 

MACH NUMBER 

0.287SE 01 

0.1AA1F 

00 

• 0.15B3F 00 

MASS FLOW RATE.POUND/MIN 

0.623QE 01 

0.62 30E 

07 

• 0.6B53F 0? 


THF PRIMARY STREAM AREA AFTFR ACCOMOOAT ION-0. 3262E-0* SQ. FFFT 
TOTAL PRESSURE RISE-0. B384E 0? POUNO/SO. FOOT 
DIMENSIONLESS PRESSURE R I SF-O. 3988E-01 
DIMFNSIONLESS static pressure CHANGE-0. ?25lF-0l 
MOMFNTU** PATIG-O.P049E 00 

**** C ASF 2 ENTRAINMFNT PATIO- 13.00000 *••• 


MACH N'J“RER ITERATIONS- 1 
PRf M A°Y AREA I TFR AT IONS* 0 
PRFSSURF ITFRATIONS* A 



ACCOMOOAT I ON PFGION 
PRIMARY SECONDARY 

PRESSURE, PS1A 

0.1A24F 

0? 

0.1474F 02 

TPMOf raTURF 

0.?oR7 F 

03 

0.5762F 01 

VFLOCITY. FT/SEC 

0.?4*l c 

0* 

0.2J74P 0* 

MATH »|1JMB C P 

0.2RR2F 

01 

0.1RRRF 00 

MASS FlOW c AT F . POUND/M IN 

0.6230E 

01 

0. RODS'' (p 


TH c PPTMAPY S TP F AM A B F A AFTFR ACT OMOP AT 1 ^N = 0. 3 301 F-07 SQ. 
TOTAL PRF ^ SURF D TSE=0. B 251F 0’ POUN^/SO. FT OT 


r F F T 


tlUTLFT 

0.1476F 07 
0.*433F o? 
0.2272F O'* 
0. 1 9P7F 00 
n.pnrr or 


0 | *iS 1 UNI c S* PR c S S' JR F ° I SP = 0. n?cp_oi 


to 

CD 


'MM* NSlflNl rs^ STATjr PQCSSURF r.n ANC.e = O . 11 0? F - 01 
MOMENTUM PAT rO.Q.' 107F 01 
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**** CASE 3 FNTRAINMFNT RATIO' 16.00000 **** 

MACH NUMBFr ITERATIONS' 1 
PRIMARY AREA I T FR A T I ONS' 6 
PRESSURE ITERATIONS' A 



ACCOMOOAT] 

PRIMARY 

! ON REGION 
secqno ary 

• outlet 

PRESSURE. PSIA 

0.1A05E 

02 

0.1A0SE 

o? 

• 0.1A55E 02 

TEMPERATURE 

0.2975F 

03 

O.S2A2F 

03 

• 0.5383E 03 

VELOCITY, FT/SEC 

0.2AA6E 

04 

0.26A1E 

03 

• 0.2773E 03 

MACH NUMBER 

0.2891F 

01 

0.23S2E 

00 

• 0.2A37E 00 

MASS FLQH RATE,PQUNP/MlN 

0.6230E 

01 

0.9968F 

02 

• 0.1059E 03 


THE PRIMARY STREAM AREA AFTER ACCOMODAT ION-0. 3331 E-03 SO. FEET 
TOTAL PRESSURE RISE'0.B131E 02 POUNO/SQ. FOOT 
DIMENSIONLESS PRESSURE R I SE«0. 3867E-01 
DIMENSIONLESS STATIC PRESSURE CHANGE =-. 3367F-02 
MOMENTUM RATI 0-0.1 92 7 E 01 

**** CASE A ENTRAINMENT RATIO- 19.00000 **•* 

MACH NUMBER ITERATIONS' 2 
PRIMARY AREA ITERATIONS' 0 
PRESSURE ITERATIONS' 6 




ACCOMOOATION region 


• OUTLET 


PRIMARY 

SECONDARY 


PRF SSURE » PS ! A 

0.1381F 

02 

0.1381E 

0’ 

• 0.1429E 02 

TEMPERATURE 

0.2961F 

03 

0.5216E 

03 

• 0.5336E 03 

VELOCITY, ET/SEC 

0.2449E 

04 

0.31 77F 

03 

• 0.3291E 03 

MACH NUMBER 

0.2903E 

01 

0.2836E 

00 

• 0.2905E 00 

MASS FLOW RATEf p OUNn/MTN 

0.6230E 

01 

O.U 84F 

03 

• 0.12A6F 03 


THF PRIMARY STREAM AREA AFTFR ACC OMOO A T ION'0 . 336H E-03 SQ. FFFT 
TOTAL PRFSSUOF RISE'O.ROTOE 0? POUNO/SO. FOOT 
0 I MF NS I ONL ESS PRESSURE R I S E = 0 . 3 R1 9 F- 0] 

OIMFNSIONLF'S STATIC PRESSURE C HANGF =- . 20 ROF -0! 


MOMENTUM PAT ln=0.?6RRF 01 
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**** tasf s fntrainmfnt KATjn= ??.oooon *♦** 

MACH NUMRF R I T Fp AT I0NS = ? 

PRIMARY ARFA ITFPAT!nNS= 0 


PRESSURE ITERATIONS 3 * 



accomodation region 


♦ outlet 


PRIMARY 

SECONnARY 

* 

PRESSURE, PSI» 

0.1 1 El E 0? 

0.1 3S1C 

07 

* 0.J399F 07 

TEMppuiTURF 

0 . 2 Q 4 3E 03 

0.S1R4F 

03 

« 0. c ?BQF 07 

VELOCITY, FT/SEC 

0.24S4E 04 

0.373SE 

03 

% O.’P’AF 03 

MACH NUMBER 

0.201 7E 01 

0.3345F 

00 

% 0.3399E 00 

MUSS EL DM PATE, POUND 7 M I N 

0.6230F 01 

0.1 3 71 E 

03 

• 0-1433F 03 


THF PRIMARY STREAM AREA AFTER ACCnM0DATinN=0.3*14E-CP SO- FEET 

TOTAL PRESSURF RIS6=0.7Q46E 0? POUND/SQ- FOOT 

DIMENSIONLESS PRESSURF R I SE=0. 3780E-01 

DIMENSIONLESS STATIC PRESSURE CHANGE=“. 41 9TF-01 

MOMENTUM R AT 1 0=0 . 7 5°3E 01 

**** CASE * FNTPAINMPNT RATIO= 29.00000 **** 

MACH NUMRFR itepations= ? 

PRIMARY AREA ITERATIQNS= 0 
PRESSURE 1 TFR AT tON$= 10 




ACCOMODATION REGION 

* OUTLFT 


PRIMARY 

SECONDARY 


PRFSSUREtPSIA 

0.131 5E 

02 

o.niEF 02 

• 0.1362E 02 

TCMPFRATUPF 

0.20’OE 

03 

0.S144F 03 

• 0.92*7? 03 

VELOCITY, FT/SEC 

0. 245R F 

04 

0.4326F 03 

% 0.A40RF O’ 

MACH NUMB FP 

0.2O3SE 

01 

0.3889E 00 

♦ 0.3927F 00 

MASS FLOW RATE,70uN0/MIn 

0.4230R 

01 

0. 1. 5 *'8E 03 

% 0.1620F 0’ 


TH C PRIMARY ^TRFAM ARFA AFTER ACC OMO^A T I HN=0 . 347? F-O’ SO. FFFT 
TOTAL PPFSSURF PISE=0.7RB4E 0? POljNP/SQ. F P n T 
01 MPNS IONt F$S p»FSSIJPE R I SE = 0. 37S0F-0 1 
0 1 MF NS I n NI F S S STATIC PRFSSUP^ C HA NGF (S 71 7F - 01 

MOMCNTIJM P AT TO = O.^A60F 0! 
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***« r.ASF 7 ENTRAINMENT RATTO= 70.00000 **** 


MATH N'l^npR ITERATIONS* 3 
PRIMARY AREA T T F R AT I DNS = 0 


PRESSURE I TFRAT l ON S = 




ACCOMODATION REGION 


* 

OUT L r T 


PRIMARY 

SECONDARY 

% 



PRPSSURF T PS 1 A 

0.1772E 

0? 

0.1272F 

O'* 

% 

0.1 ?1»E 

0? 

temperature 

0.7RR2F 

03 

0. 509SF 

07 

* 

0. k ’70F 

03 

veldcity»ft/sec 

0.2466F 

04 

0.496SC 

0^ 

* 

0.5024F 

0? 

MACH NUM0FP 

0.PR57F 

01 

0.44fl6E 

00 

* 

0.4S02E 

00 

MASS FLOW RATE, POUNO/MIN 

0.6230F 

01 

0.1T44F 

07 

% 

0.1 R07F 

o-> 

THF PRIMARY STREAM AREA 

AFTER ACCOMOOAT ION=0 

. 3F47E-0’ 

SO. 

FEET 




TOTAL PRFSSURF R I SE = 0 .7048E 0? POUNO/SQ. FOOT 
01 HFNSinNLESS PRESSURE R I SE = 0 . 77? 3E-01 
OIHPNS IONIESS STATIC PRESSURE CHANGF^-. 07 34F- 01 
MOMENTUM P AT I 0=0. S90RF OT 

C A SF 0 FNTR AT NMF NT RATIO* 31.00000 **** 

MACH NUMRFR ITERATIONS* 3 
PRIMARY AREA I TFR AT TDNS= 0 
PRESSURE ITERATIONS* 16 




accomodation region 

* OUTLET 


PRIMARY 

SECONDARY 


PRESSURE, PSI A 

0.1 219E 

0? 

0.1219E 

0? 

* 0.1264F 02 

TEMPERATURE 

0.2R17P 

O'* 

0.3033F 

07 

♦ 0.3111c 03 

VELOCITY, FT/SEC 

0.2475F 

04 

0.56F9F 

O'* 

» 0.F707E 03 

MACH NIJMREr 

0. 29949 

or 

O.ll 33F 

00 

* 0.F143F 00 

*4SS *L n W RAT E , POUNO/M IN 

0 . 6 7 30F 

01 

0.193IF 

07 

* 0.1°04 c 03 


TH C PRIMARY STRFAM AO F A A F TFP ACC OMDOA T IDN = C . 344* F-0* SO. FFF T 
rnTAl PPFSSUPF PISF=0.7fl^0F 0? PmiNP/SQ. FpnT 
D l ME N* 1 n^( F SS PPPSSU«F R I S F*0 . 0 7 ■> QT -Q 1 
01 MF NS I UNI T S S STATIC PMFSSURE CHANGE -- . 13 A OF CO 


»Jin V|T Nj Ti |m o aT IM*0 # 7 r 01 



SOI 


**** T4CF C FAITPAINM^NT r at t n= . 00000 ♦ *** 


MACH N» F R TTFOATJOM^- 

■a 






P P I m A R y APEA [TPRATInMC 

= 0 






p P c SSUP r I T F p A T 1 0^1 5 = 7i 








AC CflMOOATI ?N R F G I ON 


* 

OUT L FT 


PRIMARY 

SFCflN n ARY 

* 



PRF S SI#R «= » °S 1 A 

0.J’90 C 02 

0.11 30 c 

07 

♦ 

O.i 1 one 

0? 

tcmppratuRF 

0.2R10P 03 

0.4951= 

OT 

* 

0.50?6F 

07 

VELOCITY, FT/SEf 

0.RAR6P 04 

0.54 73F 

OT 

* 

0.F474F 

Q7 

MACH NUMMFR 

0.3024F 01 

0.5937F 

00 

• 

0. c 88°F 

no 

MASS FLOW RATF*PnUWO/Mpt 

0.6230F 01 

0.2113F 

07 

♦ 

0.?lfll c 

0^ 

THF PRIMARY STRFAm AREA 

aftfr AcrnMOPATioN=o 

. P7R0F-0’ 

SO. FFFT 




TOTAL PPFRSURF P I ">£ = 0 . 7373 F 0? POUNP/SQ. FPOT 





n I MP NS ! ONL E S3 PPPSSURC 

R I SF=0. T^fSF-o! 






Q I M c NS I ONI F S ^ STATIC ppFSS'JPF CHAnGFs- . \ 79fcF 

00 





MOMENTUM pATin=0.011 C F 

o 1 







**** TASF 

10 FNTRAINMFNT PATIO- 

77.00000 

*♦** 

MACH N'Jmp^R ITFPATI0NS = 

■> 






primary APFA I TF D A T I On<s 

= l 






PRFSSU^E I TFPAT ! f, NS= 








ACCflMnOATinN RFGIf)^ 


* 

OUT 1 FT 


PR I M A R Y 

SFfONn ARY 

* 



PR c SS'JRF t p S I A 

0 . 10 S R F 07 

0.10SRF 

O 7 

♦ 

0.1110F 

07 

TCMDFRATUR r 

0.7764P O’ 1 

0. «>fl ^4F 

07 

* 

0.^°l ? F 

0-5 

V c LOCI TY.FT/^Ff 

0.?^07f 04 

0. 74 R7F 

O’ 

* 

0. 74 1 ^ F 

07 

MArn n'imrcp 

0 . 1 0 R0 C 0- 

0. 4 O 3 Qp 

00 

* 

0- 4 °71 F 

00 

MASS Fl -1W PATF, POUNf /M IN 

0 . A. *> i O r 01 

0. 57 0 5F 

07 

* 

0.7^(-7F 

07 


TH C PRI^APY STRF A w APF* \cttp Af r [ )Mnn a T 1'1 N*0 . 7oa 7 f-0 * SQ. F r r T 
T 0 T A | Ppr<su p F P 1<F = 0 .7ao->c 07 i»niNn/S0. r f M )T 
o i nr»js I i c s^ pf c ssupf o I s r : o . ’-vnr-ni 
nf |HM| c S T A T I r P” r S S l JP r r u a mt, r r - . 7 » Q 7r 00 

MfiMPfjTijM P 4 l[Msy,n^F 0 ? 
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**** CASE 11 FNTRAINMENT RATIO 40.00000 **** 

MACH NUMRFR ITERATI0NS= 4 
PRIMARY AREA I T ER AT I ONS = l 
PRESSURE I T ER AT I ONS = 80 

ACCOMODATION REGION * OUTLET 

PRIMARY SECONDARY * 


PRESSURE, PSI A 

0.8932E 

01 

0.B932F 

01 

* 

0.97?7F 

01 

TEMPERATURF 

0. 2*1 4F 

03 

0.46 06F 

03 

♦ 

0.4723E 

03 

VELOCITY, FT/SEC 

0.2F33F 

04 

0.91.37E 

03 

♦ 

0.8776E 

03 

MACH NUMBER 

0.3194F 

01 

0.8681F 

00 

* 

0.8234E 

00 

MASS FLOW RATE, POUND/M IN 

0.6230E 

01 

0.2492E 

03 

* 

0.2554F 

03 


THE PRIMARY STREAM AREA AFTER ACCOMOD ATI 0N=0. 4*45 F-03 SQ. FEET 
TOTAL PRESSURE R I S E=0 . 8 36S E 02 POUNO/SQ. FOOT 
OIMFNSIONLESS PRESSURE R I SE=0. 3979E-01 
DIMENSIONLESS STATIC PRESSURE CHANGF=-. 3337E 00 
MOMENTUM P AT 1 0 = 0 . 1 420E 0? 

******************************************** 
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'«S0. „ I 


D l MF NS I ONI ESS SOLUTION USING THE SF INITIAL CONDITIONS, 
PRESSURE RATI 0= TO • 00 TEMPFRATURE RATIO* 1.50 ( ATH*C W/ AM ) =0 . 001 000 


******** 

********* 

***** 

*********** 

* * 

****** 

********** 

******* 

ENTRAINMENT 

RATIO 

* DIMENSIONLESS 

• PRESSURE RISE 

• 

TOTAL • 

• 

DIMENSIONLESS static 
PRESSURE CHANGE 

* 

• 

MACH 

NUMRFR 

• DIMENSIONLESS 

* KINFTIC fnergy 

momentum 

ratio 

10.0 

• 0.03988 

• 

0.02251 

• 

0.1 F5 

• 0.0173 

0.805 

1 °.o 

• 0.03925 

• 

0. 011.02 

* 

0.199 

• ,0.0? 80 

1 .302 

16.0 

• 0.03867 

• 

-0.003’7 

• 

0.244 

• 0.0414 

1 .927 

1°.0 

• 0.03819 

• 

-0.02089 

• 

0.29? 

* 0.0579 

2.688 

2?.0 

* 0.03780 

* 

—0. 04! 97 

• 

0.340 

* 0.0775 

3.593 

25.0 

* 0.03750 

• 

-0. 0671 7 

• 

0.393 

• 0.1007 

4.660 

78.0 

* 0.03733 

• 

-0.097’ A 

• 

0.450 

* 0.1281 

5.908 

31 .0 

♦ 0.03729 

• 

-0.13397 

• 

0.5?4 

* 0.1604 

7.373 

34.0 

* 0.037*6 

• 

-0.17957 

• 

0.589 

• 0.1992 

9.115 

3T.0 

• 0.03801 

• 

-0. 2397? 

• 

0.682 

* 0.2476 

11.259 

‘0.0 

* 0.03979 

• 

-0.33376 

• 

0.823 

• 0.316? 

14.205 

******** 

********* 

***** 

*********** 

* * 

****** 

********** 

****** 


DIMENSIONAL SOLUTION USING THESF INITIAL CONDITIONS, 

PRESSURE P AT 10=30.00 TEMPERATURE RATIO= 1.50 ( ATH*CW/ AM ) =0. 00? 000 

SECONDARY PP c SSURE= 1 A. 60 PSI SECONDARY TEMPERATURE^ 530.00 DEGREE RANKlNE MIXING TURE AREA=0. 0"7?60 SO. FT. 

********************************************************** 
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ENTRAINMENT 

RATIO 

♦ 

♦ 

EXIT STATIC 
PRFSSURF (PSD 

* 

* 

TOTAL PRESSURE 
R I S F IPSF) 

* 

* 

STATIC PRESSURE 
CHANCE (PSP) 

* 

* 

KINETIC FNFRGY 
(PSF) 

10.0 

* 

14.93 

* 

8^.84 

* 

47.3240 

* 

36.70 

13.0 

* 

14.76 

* 

8?. 51 

* 

23.168* 

* 

58.76 

i6.0 

* 

14.55 

* 

»1 .31 

* 

-7.0784 

* 

37.09 

1R.0 

* 

1 4. ?o 

* 

80.30 

* 

-43 . QR4 4 

* 

1 21 .64 

22.0 

* 

17.09 

* 

70.46 

* 

-8R.2290 

* 

1 52.93 

? 6 . 0 

* 

13.6? 

* 

78.84 

* 

-1 41.2129 

* 

2H.76 

'XJ 

TO 

• 

o 

♦ 

13.18 

* 

78.48 

* 

-204.651 8 

* 

269.21 

31.0 

* 

12.64 

* 

78.39 

* 

-931 .6621 

* 

337.16 

?4.0 

* 

11.93 

* 

78.75 

* 

-377. 5278 

* 

41 8.71 

37.0 

* 

11.10 

* 

79 . °? 

* 

-504.013? 

* 

5 ?0 . 54 

o 

• 

o 

* 

9.73 

* 

83.65 

* 

-701 . 6 53 ■a 

* 

*64.84 
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Jet Pump Total 

Mixing Tube . 

Mixing Tube 
Dynamic Head 
psf 

/n V 2 \ 

Mixing Tube 

Pressure Rise 

Outlet Static 

Exit Mach No. 

psf 

Pressure 

r m V m l 

M 

(AP t ) 

psia 


m 

(P ) 
' nr 




10. 0 

101.5 

15. 32 

40.83 

. 163 

13. 0 

99. 59 

15. 13 

65. 30 

. 206 

16. 0 

97. 86 

14. 90 

96. 04 

. 255 

19. 0 

96. 35 

14. 62 

133.47 

. 300 

22. 0 

95. 05 

14. 29 

178. 30 

. 352 

25. 0 

93. 96 

13.88 

231.09 

.407 

28. 0 

93. 12 

13.40 

293. 75 

. 469 

31. 0 

92. 61 

12. 81 

368. 10 

. 535 

34. 0 

92. 60 

12. 06 

458. 39 

. 615 

37. 0 

93. 52 

11. 03 

573. 92 

. 722 


Table 2 - Jet Pump Performance Characteristics from Computer Output 
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Nozzle 

Nozzle 

Throat 

Diameter 

inches 

Nozzle 

Throat 

Area 

ft 2 xl0 4 

Nozzle 

Type 

Design 

Pressure 

psia 

Design 

Temp. 

°F 

Flow 

Coefficient 

C 

w 

Flow Rate at 
Design 

W 

P 

Case 2 

0.047 

1.089 

converging- 

diverging 

350 

200°F 

.938 

6. 39 lbm/min 

Case 3 

0. 1109 

6.05 

converging- 

diverging 

100 

1200°F 

.965 

6. 52 lbm/min 

Case 4 

0.0596 

1.75 

converging- 

diverging 

350 

1200°F 

.965. 

6. 60 lbm/min 

Case 4A 

0.063 

1.952 

converging 

350 

1200°F 

.935 

7. 11 lbm/min 

Case 

LD-4 

0.0658 

1. 650 

converging- 

diverging 

350 

1200°F 

.941 

6. 08 lbm/min 


Table 3. Primary Nozzle Characteristics 
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■m.. 



Flow 

Parameter 

Instrumentation 
Used to Measure 
Parameter 

How 

Recorded 

Required for Determining 

Data Reduction Procedure 

Primary Flow 

P 

po 

Bourdon Tube Gage 

Manually 

and 

Photographically 

Jet Pump Input Conditions 

None needed 


T 

po 

Thermocouple 
and Bridge 

Manually 

Jet Pump Input Conditions 

None needed 


W 

P 

Orifice Flow Meter 
and 

Panel Gage 

Manually 

Jet Pump Input Conditions 

Standard calibration curves 
provided by flowmeter manu- 
facturer 

Secondary Flow 

T 

so 

Dial Gage 
In Suction Duct 

Manually 

Secondary Flow Temperature 

None needed 


P atm 

Mercury Barometer 

Manually 

Atmospheric Pressure 

None needed 1 

' 


P b 

Manometers 

Manually 

and 

Photographically 

Secondary Flow Rate 

See below 


w 

s 

Calibra ted 
Bellmouth 

Manually 

Secondary Flow Rate 
in lb/ min 

Equation (91) 

Mixing Tube 
and Diffuser 

P vs. length 

Manometer Board 

Photogra phically 

Mixing Tube and Diffuser 
Static Pressures 

None needed 


P mo | 
P do ) 

Kiel Probe 
Traverse 

Manually 

Discharge Stagnation Pres- 
sure for Jet Pump 

"Mass-momentum'' method 
for Compressible Flow 
(See Text) 


T 

eo 

Dial Gage Near 
End of Diffuser 

Manually 

Jet Pump Discharge 
Temperature 

None needed 


Table 4. Measured Parameters and Instrumentation 
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Test Series #1: Reference Test 


Tes 

t No. 

^ Un Nozzle Set 
No. 

Primary Flow 
Pressure 

Primary Flow 
Temperature 

Secondary 
Flow Rate 

Discharge 

Configuration 

Cluster 

Configuration 




(psia) 

(° F) 






1 

1-5 4 

1 

350 

1085° 

4 values 

NAS2-251! 

3 mixing tube 

Original 


2 

6-9 


300 

1160° 

4 values 



1 


3 

10-14 


260 

1185° 

5 values 





4 

15-18 


225 

1145° 

4 values 



♦ 


Test Series #2: Nozzle Cluster Performance Comparison 


5 

19-22 

4 


350 

1130° 

6 

23-26 



300 

1160° 

7 

27-31 



260 

1165° 

8 

32-35 



225 

1158° 

9 

36-39 

4 

350 

1160° 

10 

40-43 



300 

1180° 

11 

44-47 



260 

1200° 

12 

48-51 



225 

1200° 

13 

52-55 

LEW2-4 

350 

1130° 

14 

56-59 



300 

1195° 

15 

60-63 



260 

1190° 

16 

64-67 



225 

1180° 


4 values Constant-area 

I short mixing tube 


Original 


4 values 


4 values 


LD#1 


I 

LD#2 


Test Series #3: Final Performance Tests 


17 

96-99 

18 

100-103 

19 

104-107 

20 

108-111 

21 

80-83 

22 

84-87 

23 

88-91 

24 

92-95 

25 

68-71 

26 

72-75 

27 

76-79 


4 


350 

300 

260 

225 

300 

260 

225 

175 

260 

225 

175 


1110 ° 

1150° 

1200 ° 

1200 ° 

750“ 


455° 

460° 

460° 


4 values Constant-area 

j long mixing tube 


LD#1 


Table 5. Jet Pump Test Program 
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Test No 

. Nozzle Set 

No. 

Primary Flow 
Pressure 

Primary Flow 
Temperature 

Secondary 
Flow Rate 




(psia) 

( 

° F) 



28 

242-245 4A 

350 

1169° 

4 values 


29 

246-249 (converging) 

300 

1150” ■ 

-1100° 




30 

250-253 


260 

1150” 




31 

254-257 


225 

1150° 




’ 32 

140-143 


350 

1160” 




33 

144-147 


300 

1150° 




34 

148-151 


260 

1095”- 

-1031° 




35 

152-155 


225 

1020° 



_ 

36 

124-127 


300 

750° 



1 

37 

128-131 


260 


1 




38 

132-135 


225 

1 





39 

136-139 


175 

745° 




40 

112-115 


260 

455° 




41 

116-119 


225 

455° 




^42 

120-123 

i 

175 

460° 

i 



r 43 

180-183 

2 

350 

200° 

4 values 


44 

184-187 


300 






45 

188-191 


260 






46 

192-195 


225 


f 




47 

196-199 


400 

455° 



i 

48 

200-203 


350 

455° 




49 

204-207 


260 

450° 




50 

2 08-211 


225 

450° 




51 

212-215 


400 

460° 




52 

216-219 


350 

455° 




53 

220-223 


260 

450° 




.54 

224-227 


225 

450° 


i 


55 

168-171 

3 

100 

1160°- 

-1115° 

4 values 


50 

172-175 


85 

1140° 




57 

176-179 


70 

1145° 




58 

156-159 


85 

755° 




59 

160-163 


70 

750° 




60 

164-167 


55 

745° 



(61 

228-231 


100 

1150°- 

-1110° 




C 62 

232-235 


85 

1155° 




(63 

236-239 


70 

11 

r r 0 

55 




Table 5. Jet Pump Test Program (Continued) 


Discharge 

Configuration 


Constant-area 
short mixing tube 


Constant-area 
long mixing tube 


Constant-area 


long mi 


xing tube 


Constant-area 
short mixing tube 


Constant area 
long mixing tube 


Constant-area 
short mixing tube 


•«£!. 

Cluster 

Configuration 

LD#1 


LD#1 


LD#1 



Test 

Date 

Test 

Nos. 

Run 

Nos. 

Constant 

Area 

Mixing 

Tube 

Nozzle 
Case No. 

Tabulated 

Data 

Static 

Pressure 

Parameter 

P -P 
m so 

Stagnation 

Pressure 

Rise 

Parameter 

P -P 
mo so 

Velocity 

Profiles 

Static 

Pressures 

vs. 

P and m 
po 

Static 

Pressures 

vs. 

Distance 

Thrust 

Augmentation 

Parameter 

T 

P 

SO 

P 

SO 

25 July 
1968 

1-4 

1-18 

short 

4 

no extension 

Table 8 


Table 7 
Figure 48 

Figure 44 
Figure 45 

Figure 46 

Figure 47 


31 July 
1968 

5-8 

19-35 

short 

4 

no extension 

Table 9 

Figure 48 

Table 7 
Figure 48 

Figure 49 

Figure 50 

Figure 47 


13 August 
1968 

* 

9-12 

36-51 

short 

4 

Table 10 

Figure 52 

Table 7 

Figure 53 

Figure 54 

Figure 55 


19 August 
1968 

13-16 

52-67 

short 

LD#2-4 

Table 11 

Figure 56 

Table 7 

Figure 57 

Figure 58 

Figure 55 


28 August 
1968 

17-27 

68-111 

extended 

4 

Table 12 

Figure 59 
Figure 60 

Table 7 
Figure 59 

Figure 61 
Figure 62 

Figure 63 

Figure 64 

Table 12 

13 Sept. 
1968 

28-31 

242-259 

short 

4A 

Table 13 

Figure 65 

Table 7 

Figure 66 

Figure 67 



3 Sept. 
1968 

32-42 

112-155 

extended 

4A 

Table 14 

Figure 68 
Figure 69 



Figure 70 

Figure 71 

Table 14 

11 Sept. 
1968 

43-50 

180-211 

extended 

2 

Table 15 

Figure 72 

Table 7 

Figure 73 

Figure 74 

Figure 75 

Table 15 

11 Sept. 
1968 

51-54 

212-227 

short 

2 

Table 16 

Figure 76 

Table 7 

Figure 73 

Figure 77 



10 Sept. 
1968 

55-60 

156-179 

extended 

3 

Table 17 

Figure 78 

Table 7 

Figure 79 

Figure 80 

Figure 81 

Table 17 

12 Sept. 
1968 

61-63 

228-241 

short 

3 

Table 18 

Figure 82 

Table 7 

Figure 83 

Figure 84 




Table 6 

Index to Test Results 
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Test 

No. 

Run 

No. 

Traverse 

Location 

and 

Number of 
TraveraM 

Nozzles 

Cluster 

Con- 

figuration 

Duct 

Con- 

figuration 

P 

po 

psia 

T 

po 
* F 

m 

P ,lo 

psf 

Diffuser 

Experimental 

Anahtienl 

l> . |> WP i -iP j 

niU p — — ’*nio "''nm' ana * j Comments on Results 

1 ’ nnal ex |P ) - P 

^oTmTSI m0 anal 80 1 

ld “‘ . Friction t»t ■ \ - 

P 

mo 

psf 

| P -1> 

p mo so 

so P 

80 

pflf 

ixing Tube — *■ 

1 

1 

Diffuser 

l 

Case 4 

Origt 

Clus 

nal 

er 

Sho 

Conve 

Jtxinar 

rt 

rging 

JZlltUL 

350 

10W 

2G. 2 

2169. 5 

2187.9 

2093.3 .04519 

.0538 1 0503 : 111.70 *>1 ». 2 

short mixing tube 
one Inverse 
incomplete mixing 

1 

2 

Diffuser 

1 

Case 4 





350 

1086 

22.5 

2171.2 

2184.0 

2100.3 

.03985 

.0540 ' .0514 24.26 *22. 5 



1 

4 

Diffuser 

1 

Caae 4 





350 

1084 

18.4 

2180.6 

2189. • 

2107.0 

. 03891 

.0555 .0536 30.95 *27.4 


1 

5 

Diffuser 

1 

Case 4 





350 

1085 

15.8 

2146.9 

2203.0 

2110.7 

.04372 

.0561 .0547 | 23.18 *20. 07 


> 

7 

Diffuser 

1 

Case 4 





300 

1158 

154 

2169.4 

2180.8 

2102.4 

. 03729 

.M59 1 .0435 13.06 .14.2! 


3 

11 

Diffuser 

1 

Case 4 





260 

1186 

29.7 

2157.8 

2169.1 

2102.2 

.03182 

1 

.0391 .0387 10.28 413.30 


5 

10 

Diffuser 

1 

Caae 4 

Original 

Cluster 

Short Con- 
stant Area 
Mlxii*T\ib+ 

350 

1131 

22.4 

2214.6 

2225.8 

2128.0 

.04595 

• 0M1 ’ .0515 11. !1 .10.7! 


• 

14 

Diffuser 

2 

Case 4 





300 

1174 

24.9 

2216.4 

2225.8 

2129.4 

.04527 

M60 .; .0437 , -3.34 - 3.59 1 


7 

n 

Diffuser 

1 

Caae 4 





140 

1184 

24.13 

2192.5 

2201.8 

2128.1 

.03463 

.1 

.0393 ! .0371 

5.26 

+ 8.88 


' 10 

41 

Dtffuser 

1 

Case 4 

LDfl 

Short Con- 
stant Area 

300 

1187 

25.8 

2192.1 

2202.4 

2117.2 

.0402 

• 0458 . 0434 

! 6.96 

+ 7.10 


i u 

57 

Diffuser 

2 

Csse 

LD2-4 

LD#2 

afcrt Con- 
stant Area 
MtalMTub. 

300 

1198 

27,4 

2216.3 

2226.5 

2105.8 

.0573 

. 0427 , . 0403 

; -35. 80 

-42. 18 

Not as well mtimi as Mhar etaeasra. 
Traversa tkrsegfc aaaslea ilpiffliiMlj 
lirfir tki traverse sat through — niaa 

18 

101 

Mixing 

Titos 

2 

Case 4 

LDfl 

Cwituni 

An* 

MlilocTA* 

300 

1154 

28.5 


2204.1 

2126.7 

. 03839 

. 0455 ; . 0376 

2.57 

+ 3.22 

extended mixing tuba 

two traverses 

18 

101 

Diffuser 

2 

Case 4 





300 

1154 

26.5 

2190.4 

2202.1 

2126.7 

.0355 

.0455 .0376 

4.47 

+ 5.38 



18 

142 

Mixing 

Tt*e 

2 

Case 4 





300 

1158 

22. 1 


2206.9 

2131.25 

.0355 

.0465 

.0409 

11.51 

+13.20 



IS 

101 

Diffuser 

2 

Case 4 





300 

1158 

22. 1 

2204. 9 

2212.9 

2131.3 

. 03828 

. 0465 

.0409 

5.58 

+ 6.40 



54 

173 

Mixing 

Tikx 

2 

Csse 3 





85 

1138 

23.6 


2188.7 

2127. 8 

. 02862 

.0360 

.0298 

2.51 

4- .1. 96 



58 

173 

Diffuser 

2 

Case 3 





85 

1138 

26.6 

2178. 1 

2187.6 

2127.8 

. 02782 

. OJCO 

. 0298 

4.21 

+ 6.64 



48 

2N 

Diffuser 

1 

Case 2 





350 

444 

28.5 

2124.0 

2136.1 

2088.6 

. 02274 

. 0350 

. 0269 

8. 69 

+13.418 

extended mixing tube 
one traverse 

52 

217 

Diffuser 

1 

Case 2 



Short Con- 
stant Area 
fixing 1 ube 

350 

452 

27.3 

2156.1 

2166.6 

21189.2 

. 03704 

. 11349 

0325 

•». 4S 

13. 97 

short mixing tube 
one traverse 

63 

240 

Diffuser 

1 

Case 3 





85 

1159 

26.3 

2142.5 

2153.7 

2088. 7 

.03111 

. 0354 

032- 

3.53 

t 5.15 


29 

258 

Diffuser 

1 

Case 4A 






300 

1 1 o:j 

26.6 

2113. 

2157.4 

2093. 9 

. 031132 

.0490 

. O l Vi 

32 i.J 

+33 94 

short mixing mhe 
one traverse 
distorted irlwih profile 


Table 7. Tabulation of Stagnation Pressure Results 





Table 8 

Jet Pump Test Results 
NAS 2-2518 Mixing Tube 
Original Nozzle Cluster 
Case 4 Nozzles 
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P 

T 

Test No. 

Run 

po 

psia 

po 
° F 


° F lbm/min 


27 

26 

29 


30 


31 



Table 9 


Jet Pump Test Results 
Short Constant Area Mixing Tube 
Original Nozzle Cluster 
Case 4 Nozzles 











Test No. 

Run 

P 

po 

T 

po 



psia 

°F 


lbm/min 


6. 82 
6. 82 
6.82 
6. 82 


lbm/ min 


182. 0 
149. 5 
125. 0 
107.4 


170.0 
146. 9 

125.0 
98. 5 


160.0 
139.7 
112.5 
93. 6 


151.0 

129.0 
113. 1 

95.5 
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Table 11 


Jet Pump Test Results 
Short Constant Area Mixing Tube 
LD#2 Nozzle Cluster 
LD#2-4 Nozzles (7 Nozzles) 
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Test No. 

Run 

P 

po 

psia 

T 

po 
° F 

W P 

lbm/min 

T 

so 

°F 

W 

s 

lbm/min 

m 

T 

25 

68 

260 

453 

6.48 

78 

163.4 

25.2 

5. 06 



69 



452 

6.48 


140. 4 

21.7 

3. 64 



70 



453 

6.46 

1 

122.9 

19.0 

2.77 

1 


71 

1 

’ 

455 

6.48 

T 

110.5 

17.1 

2. 24 

2 

6 

72 

2‘ 

;5 

454 

5. 67 

78 

153.2 

27.0 

5. 05 



73 



456 

5. 67 

79 

139.3 

24.6 

4. 11 



74 



460 

5. 65 

79 

124.7 

22.1 

3.28 

i 


75 

i 


462 

5.64 

77 

104.2 

18.5 

2.28 

2 

7 

76 

V 

5 

459 

4.43 

78 

133.8 

30.2 

4. 92 



77 



459 

4.43 

79 

121.6 

27.5 

4.04 



78 



463 

4.43 

79 

106.8 

24.1 

3. 09 

1 

’ 

79 


' 

463 

4.42 

79 

94.3 

21.3 

2.42 

2 

L 

80 

300 

753 

6.63 

82 

172.3 

26.0 

4.96 



81 



753 

6. 60 

1 

153.2 

23.2 

3. 81 



82 



752 

6.59 


130.0 

19.7 

2.71 

' 

r 

83 

’ 


752 

6.59 

♦ 

111.4 

16. 9 

1.99 

22 

84 

260 

749 

5.66 

81 

161.8 

28.6 

5. 07 



85 



745 

5. 65 

80 

142.2 

25.2 

3. 81 



86 



749 

5.64 

80 

123.5 

21.9 

2. 86 



87 

! 


752 

5.62 

80 

107.6 

19.2 

2. 19 

2; 

3 

88 

225 

753 

5.03 

82 

150.4 

29.9 

4.89 



89 



750 

5.02 

80 

134.3 

26. 7 

3. 83 



90 



749 

5.01 

81 

117.8 

23.5 

2.94 



91 



748 

5.02 

81 

103.9 

20.7 

2. 28 

2 


92 

1’ 

’5 

744 

4.0 

82 

132.2 

33.1 

4. 77 



93 



745 

3.96 

80 

119.9 

30.3 

3.92 



94 


1 

745 

3.94 

83 

107.3 

27.2 

3. 15 


f 

95 


' 

748 

3.94 

80 

95. 7 

24.3 

2.50 

17 

96 

350 

1116 

6.82 

88 

179.8 

26.4 

4. 85 



97 



1115 

6. 81 

88 

153.0 

22.5 

3.35 


I 

98 


1 

1112 

6.76 

88 

129. 7 

19.2 

2. 39 


r 

99 


f 

1105 

6.77 

89 

117.3 

17.3 

1.96 

li 

3 

100 

3C 

>0 

1147 

5.78 

88 

167.3 

28.9 

4. 76 



101 



1149 

5.78 

87 

153.1 

26.5 

3.91 



102 


r 

1154 

5.77 

86 

127.7 

22.1 

2.68 

T 


103 

1 


1159 

5.73 

86 

112. 7 

19.7 

2.08 

19 

104 

2e 

;o 

1195 

5.01 

83 

156.7 

31.3 

4. 67 



105 



1198 

4.98 

83 

137.4 

27.6 

3.54 



106 


r 

1199 

4. 98 

83 

115.1 

23.1 

2.47 



107 



1203 

4.96 

80 

106.7 

21.5 

2. 13 

20 

108 

225 

1198 

4.35 

83 

145.5 

33.4 

4. 61 

I 


109 



1203 

4.35 

82 

134.7 

31.0 

3. 92 

1 


110 



1205 

4.26 

81 

119.0 

27. 9 

3.09 

t 


111 



1204 

4.30 

81 

101. 0 

23.5 

2. 22 


Table 12 


Jet Pump Test Results 
Extended Constant Area Mixing Tube 
LD#1 Nozzle Cluster 
Case 4 Nozzles 
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w 

p 

Ibm/ min 


7.25 

86 

175.3 

7.24 

86 

156.4 

7.25 

86 

134.3 

7.24 

85 

105.5 

6.31 

85 

166.0 

6. 24 

85 

148.0 

6. 23 

85 

127.3 

6.23 

84 

101.0 

6. 26 

87 

166.8 

6.30 

88 

167.7 

5.44 

85 

157.7 

5.41 

85 

142.5 

5.38 

85 

122.0 

5.38 

84 

100.5 

4.70 

84 

150.0 

4. 69 

84 

136.6 

4. 66 

84 

119.0 

4.64 

83 

95.3 


Table 13 

Jet Pump Test Results 
Short Constant Area Mixing Tube 
LD#1 Nozzle Cluster 
Case 4A Nozzles 




Test 

No. 

Run 

P 

po 

psia 

T 

po 

°F 

W 

P 

lbm/min 

1 

SO 

S F 

w 

s 

lbm/min 

m 

T 

40 

112 

260 

453 

6.96 

78 

165.6 

23.8 

4.93 



113 



457 

6. 92 

78 

142.8 

20. 6 

3. 57 



114 



459 

6.91 

79 

125.2 

18.1 

2.72 



115 



462 

6.89 

79 

108.3 

15. 7 

2.05 

41 

116 

225 

456 

6.05 

79 

155.4 

25.7 

4.96 



117 



455 

6.02 

80 

138.5 

23. 0 

3.89 



118 



456 

5. 98 

80 

124.3 

20.8 

3. 12 



119 



457 

5.98 

80 

103.2 

17.3 

2. 16 

42 

120 

Vi 

rs 

458 

4.81 

80 

140.2 

29.1 

5. 06 



121 



462 

4.81 



123. 5 

25.7 

3. 89 



122 



463 

4.78 



107.3 

22.5 

2. 94 



123 



464 

4.78 



96.5 

20.2 

2.38 

36 

124 

3( 

)0 

753 

7.10 

84 

172. 1 

24.2 

4. 73 



125 



750 

7.10 



150.4 

21.2 

3. 50 



126 



749 

7.08 



128. 7 

18.2 

2. 54 

1 


127 



750 

7.07 



111. 7 

15. 8 

1. 92 

37 

128 

2( 

;o 

751 

6.20 

85 

161.9 

26.1 

4. 75 



129 



752 

6.19 



140.5 

22. 7 

3. 49 



130 



754 

6.17 



123. 7 

20.1 

2. 70 



131 



756 

6.17 



106.8 

17.3 

2. 01 

38 

132 

22 

i5 

745 

5.35 

86 

152. 6 

28. 5 

4. 86 



133 



743 

5.35 

86 

134.2 

25. 1 

3. 69 



134 



744 

5.36 

85 

117. 8 

22. 0 

2. 82 



135 



746 

5.35 

85 

103.2 

19.3 

2. 17 

39 

136 

175 

739 

4.26 

86 

137. 0 

32.2 

4.94 



137 



743 

4.24 



119.7 

28.2 

3.73 



138 



746 

4.23 



105.4 

24. 9 

2. 89 


f 

139 


r 

748 

4.22 



95. 9 

22. 7 

2. 40 

32 

140 

3i 

.0 

1175 

7.26 

92 

174. 6 

24.1 

4. 34 



141 



1162 

7.18 



151.1 

21.1 

3. 17 



142 



1158 

7.17 



124. 5 

17.4 

2. 13 



143 


r 

1140 

7. 18 



115.5 

16.1 

1. 83 

33 

144 

300 

1146 

6.26 

94 

165.1 

26.4 

4. 43 



145 


i 

1149 

6. 23 



145.2 

23.3 

3.35 



146 



1148 

6.28 



124. 6 

19. 8 

2.43 



147 


r 

1124 

6.27 



110.7 

17. 7 

1. 93 

3 

4 

148 

260 

1094 

5.52 

94 

156. 8 

28.4 

4. 60 



149 



1057 

5.55 



129.7 

23.4 

3. 05 



150 



1048 

5.56 



112.8 

20.3 

2.30 



151 



1031 

5.58 



105.1 

18. 8 

1. 99 

3 

5 

152 

225 

1021 

4. 95 

94 

148. 7 

30. 0 

4. 62 



153 



1021 

4. 94 



129.7 

26.3 

3. 47 



154 



1021 

4. 94 



114.5 

23.2 

2. 70 


f 

155 



1021 

4. 93 



101. 9 

20.7 

2. 14 


Table 14 

Jet Pump Test Results 
Extended Constant Area Mixing Tube 
LD#1 Nozzle Cluster 
Case 4A Nozzles 
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Test No. 

Run 

P 

po 

psia 

T 

po 

°F 

W 

P 

Ibm/mtn 

T 

so 

°F 

W s 

lbm/min 

m 

T 

43 

180 

350 

201 

6.35 

78 

158.7 

25.0 

5. 55 



181 



199 

6.38 



140.2 

22.0 

4.21 



182 



199 

6.39 



123.3 

19.3 

3. 21 



183 

1 


200 

6.39 



106.4 

16.7 

2.38 

44 

184 

300 

200 

5.51 

7! 

9 

148.8 

27.0 

5. 53 



185 



205 

5.47 



137.0 

25.0 

4. 68 



186 



204 

5.45 



123.6 

22.7 

3. 79 



187 


\ 

205 

5.44 



101.3 

18.6 

2. 53 

45 

188 

260 

201 

4. 80 

79 

138.4 

28.8 

5. 52 



189 



203 

4. 79 



126.3 

26.4 

4. 56 



190 



201 

4.78 



114.1 

23.9 

3. 70 


r 

191 


’ 

201 

4.78 



96.4 

20.2 

2. 63 

46 

192 

225 

197 

4. 14 

79 

128.4 

31.0 

5. 52 



193 



197 

4. 14 



118.9 

28.7 

4. 70 



194 



196 

4. 14 



108. 1 

26. 1 

3. 86 



195 



199 

4. 12 



93.9 

22.8 

2.92 

47 

196 

400 

456 

6.30 

8 

1 

167.7 

26.6 

5.43 



197 



453 

6. 30 



147.3 

23.4 

4. 06 



198 



455 

6.29 



129.3 

20.6 

3.09 

\ 


199 



454 

6.29 



110.5 

17. 6 

2. 24 

48 

200 

350 

444 

5. 56 

81 

158.3 

28.5 

5. 43 



201 



454 

5. 54 



140.8 

25.4 

4. 21 



202 



454 

5. 54 



126.4 

22.8 

3. 36 



203 

1 


454 

5. 52 



104.8 

19.0 

2. 31 

4 

9 

204 

260 

449 

4. 07 

82 

137.0 

33.7 

5. 50 



205 



450 

4.08 



126.7 

31.0 

4. 63 



206 



451 

4.08 



115.4 

28.3 

3. 82 

1 


207 



452 

4.08 



97.2 

23.8 

2.70 

50 

208 

225 

447 

3.58 

82 

127.8 

35.7 

5. 45 



209 



447 

3. 58 



116.2 

32.5 

4. 47 



210 



449 

3. 59 



106.5 

29.7 

3. 73 



211 



450 

3.58 



93.2 

26.0 

2. 85 


Table 15 


Jet Pump Test Results 
Extended Constant Area Mixing Tube 
LD#1 Nozzle Cluster 
Case 2 Nozzles 


Test No. 

Run 

P 

po 

psla 

T 

« po 

F 

W 

P 

lbm/mln 

T 

so 

°F 

W 

s 

Ibm/min 

m 

S 

1 

212 

4 


456 

6.25 

85 

171.4 




213 



459 

6. 26 

85 

149.6 




214 



458 

6.25 

85 

127.4 




215 



458 

6. 24 

86 

,.3.7 

15.8 

5 

•2 

216 

3 

50 

451 

5.52 

86 


29.2 



217 



452 

5.52 

85 

150.8 

27.3 



218 



456 

5. 50 

84 

128.3 

23.3 



219 



457 

5.49 

83 


17.6 

5 

3 

220 

2< 

30 

451 

4. 12 

81 

142.5 

34.6 



221 



451 

4. 12 

81 

128.6 

31.2 



222 



451 

4.09 

82 

109.5 

26.8 



223 



451 

4. 09 

82 

92.1 

22.5 

54 

224 

225 

449 

3.57 

83 

133.4 

37.4 



225 



448 

3.58 

83 

120.9 

33.8 



226 



448 

3.58 

83 

105.9 

29.6 



227 



448 

3.58 

83 

88.9 

24.8 


Table 16 


Jet Pump Test Results 
Short Constant Area Mixing Tube 
LD#1 Nozzle Cluster 
Case 2 Nozzles 
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Test No. 

Run 

P 

po 

psia 

T 

po 

°F 

W 

P 

lbm/min 

T 

so 

°F 

W 

s 

lbm/min 

m 

T 

5 

8 

156 

8 

5 

755 

6. 64 

540 

158.3 

23.8 

4. 95 



157 



758 

6.59 

541 

138.6 

21.0 

3. 75 



158 



758 

6. 59 

542 

124.3 

18.9 

3.02 



159 



758 

6.58 

543 

105.5 

16.0 

2. 18 

5 

9 

160 

70 

751 

5.44 

542 

141.9 

26. 1 

4.97 



161 



748 

5.42 

542 

128.9 

23.8 

4. 07 



162 



748 

5.43 

542 

113.4 

20.9 

3. 13 



163 



749 

5.41 

542 

98.0 

18. 1 

2.35 

6 

0 

164 

5 

i5 

743 

4.40 

542 

121.6 

27.6 

4. 68 



165 



743 

4.40 

543 

114.2 

26.0 

4. 16 



166 



744 

4.40 

543 

105.0 

23.9 

3. 52 



167 

1 


746 

4. 40 

543 

90.2 

20.5 

2. 59 

5 

5 

168 

100 

1158 

6.61 

550 

161.2 

24.4 

4. 62 



169 



1140 

6.66 

550 

146. 1 

21.9 

3. 68 



170 



1128 

6. 68 

550 

128.3 

19.2 

2. 81 



171 


f 

1116 

6.67 

550 

107.9 

16.2 

2. 00 

5 

6 

172 

85 

1141 

5.74 

549 

149.4 

26.0 

4. 58 



173 



1138 

5. 75 

550 

134.5 

23.4 

3. 67 



174 



1138 

5.76 

550 

120.7 

21.0 

2.94 



175 


t 

1138 

5.76 

550 

104.3 

18.1 

2. 20 

57 

176 

7 

'0 

1141 

4.81 

549 

133.4 

27.7 

4. 41 



177 



1146 

4. 80 

548 

120.5 

25. 1 

3. 58 



178 



1148 

4.76 

548 

110.7 

23.3 

3. 06 



179 



1154 

4.77 

549 

96.4 

20.2 

2.33 


Table 17 

Jet Pump Test Results 
Extended Constant Area Mixing Tube 
LD#1 Nozzle Cluster 
Case 3 Nozzles 
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Primary Flow 

P , T , W i — Section 1 

po po p i 



Accommodation Region 


Mixing Region 


Blowing Duct Diffuser ► 


Figure 2 

Jet Pump Definitions 


LLLliJi 



Primary Flow 


Section 1 



Section b 


V 


Figure 3 

Mixing Tube Wall Friction Approximation 



( OPTIONAL) 


( OAT/OHJAL) 
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Figure 4 

Block Diagram of the Computer Program 


















e 5 

it Specific Heat Assumption 










DIMENSIONLESS PRESSURE RISE 


PUMP CHARACTERISTICS. PARAMETER IS PRESSURE RATIO 



entrainment ratio 

Figure 6 

Jet Pump Performance Characteristics 
Dimensionless Pressure Rise (t^P^*) vs. Entrainment Ratio 

A = 0.001 
T = 1. 5 

P = 25, 30, 35, 40, 45 
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MOMENTUM RATIO- EXIT/PRIMARY 


MOMENTUM RATIO. PARAMETER IS PRESSURE RATIO 



Figure 7 

Jet Pump Performance Characteristics 
Momentum Augmentation (t) vs. Entrainment Eatio 

A_= 0.001 
T = 1.5 

P = 25, 30, 35, 40, 45 
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DIMENSIONLESS PRESSURE RISE 


PUMP CHARACTERISTICS. PARAMETER IS PRESSURE RATIO 



Figure 8 

Jet Pump Performance Characteristics 
Dimensionless Pressure Rise (AP^.*) vs. Entrainment Ratio 

A = 0. 001 
T = 3.5 

P = 40, 45, 50, 55, 60, 65 
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MOMENTUM RATIO- EXIT/PRIMARY 


MOMENTUM RATIO. PARAMETER IS PRESSURE RATIO 



10 20 30 40 


ENTRAINMENT RATIO 

Figure 9 

Jet Pump Performance Characteristics 
Momentum Augmentation (t) vs. Entrainment Batio 

A = 0.001 
T = 3.5 

P = 40, 45, 50, 55, 60, 65 
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Dimensionless Pressure Rise 



Entrainment Ratio 
Figure 10 

Jet Pump Performance Characteristics 

Dimensionless Pressure Rise (A P t *) vs. Entrainment Ratio 
A = 0.001 
T = 8.0 

P = 60, 70, 80, 90 
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Momentum Ratio - Exit/Primary 



Entrainment Ratio 
Figure 11 

Jet Pump Performance Characteristics 

Momentum Augmentation (t) vs. Entrainment Ratio 
A = 0.001 
T = 8.0 

P = 60, 70, 80, 90 
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oihensiqnless pressure rise 


pump characteristics, parameter is pressure ratio 



ENTRAINMENT RATIO 

Figure 12 

Jet Pump Performance Characteristics 
Dimensionless Pressure Rise (AP^.*) vs. Entrainment Ratio 

A = 0.003 
T = 1.5 

P= 9, 10, 11, 12, 13, 14 
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MOMENTUM RATIO- EXIT/PRIMARY 


MOMENTUM RATIO. PARAMETER IS PRESSURE RATIO 
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ENTRAINMENT RATIO 

Figure 13 

Jet Pump Performance Characteristics 
Momentum Augmentation (t) vs. Entrainment Ratio 

A = 0.003 
£ = 1.5 

P= 9, 10, 11, 12, 13, 14 
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OIMENS IONLESS PRESSURE RISE 


PUMP CHARACTERISTICS. PARAMETER IS PRESSURE RATIO 
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entrainment ratio 

Figure 14 

Jet Pump Performance Characteristics 
Dimensionless Pressure Rise (AP^*) vs. Entrainment Ratio 

A.= 0.003 
T = 3 # 5 

P= 13, 15, 17, 19, 21 




MOMENTUM RATIO- EXIT/PRIMARY 










Dimensionless Pressure Rise 


Si 



Entrainment Ratio 
Figure 16 

Jet Pump Performance Characteristics 
Dimensionless Pressure Rise (AP t *) vs. Entrainment Ratio 

A = 0.003 
T = 8.0 

P = 20, 22, 24, 26, 28, 30 
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DIMENSIONLESS PRESSURE RISE 





MOMENTUM RATIO- EXlT/PRlMARY 


MOMENTUM RATIO. PARAMETER IS PRESSURE RATIO 



Figure 19 

Jet Pump Performance Characteristics 
Momentum Augmentation (t) vs. Entrainment Ratio 

4= 0.007 
T= 1.5 

P = 4, 4.5, 5, 5.5, 6 
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dimensionless pressure rise 


PUMP CHARACTERISTICS. PARAMETER IS PRESSURE RATIO 



Figure 20 

Jet Pump Performance Characteristics 
Dimensionless Pressure Rise (APj*) vs. Entrainment Ratio 

A = 0.007 
T = 3.5 

P = 5, 6, 7, 8, 9 
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MOMENTUM RATIO- EXIT /PRIMARY 



MOMENTUM RATIO. PARAMETER IS PRESSURE RATIO 



Figure 21 

Jet Pump Performance Characteristics 
Momentum Augmentation (t) vs. Entrainment Ratio 

A = 0.007 
T = 3.5 

P = 5, 6, 7, 8, 9 
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dimensionless pressure rise 




MOMENTUM RATIO- EXlT/PRIhARY 


MOMENTUM RATIO. PARAMETER IS PRESSURE RATIO 




ENTRAINMENT RATIO 


Figure 23 

Jet Pump Performance Characteristics 
Momentum Augmentation (r) vs. Entrainment Ratio 

A = 0. 007 
T = 8. 0 

P = 9, 10, 11, 12, 13, 14 
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DIMENSIONLESS PRESSURE RISE 


PUMP CHARACTERISTICS. PARAMETER IS AREA RATIO 



Figure 24 

Jet Pump Pressure Rise in Relation to 
Geometry (3^) and Mixing Tube Exit Mach Number 
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MOMENTUM RATIO. PARAMETER IS AREA RATIO 



Figure 25 

Jet Pump Momentum Ratio in Relation to Geometry (A^) 
and Mixing Tube Exit Mach Number 
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A R) D t A P 

Total Pressure Change — — or 




Figure 26 

Jet Pump System Operating Point 
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Figure 29 

Chart for Determining (Equation 66) 
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Entrainment Ratio 
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1 2 3 4 5 

Blowing Duct Diffuser Area Ratio A^/A^ 

Figure 32 

Influence of Blowing Duct Diffusion Upon Entrainment Ratio 
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Jet Pump Thrust Augmentation, 





Figure 33 

Jet Pump Thrust Augmentation Characteristic 
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System Thrust Augmentation, 


i ii i iiiiiiiii 



0 1 2 3 4 5 


Blowing Duct Diffuser Ratio, A^/A m 


Figure 34 

Influence of Blowing Duct Diffusion Upon System Thrust Augmentation 
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Entrainment Ratio, m 


Figure 35 

Relationship of Thrust Augmentation to Entrainment Ratio 



Primary Air Flow 
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Extension 


Original 
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Section 


Throttling 

Cone 


Figure 36 
Jet Pump Test Rig 
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Figure 37 

Original Mixing Tube Geometry 







Static Pressure Taps 




Figure 40 

Dynatech Jet Pump Test Instrumentation 





Stagnation Pressure Loss - psf 



Bellmouth Inlet; 


O No Cluster or Elbow 

Q Original Cluster and 
Nozzles 

A Original Cluster with 

2 Inch Nozzle Extensions 

□ Low Drag Cluster #2 


NAS 2-2518 Suction Duct; 

V Original Cluster and Nozzle 
Q No Cluster or Elbow 


40 50 60 80 100 


150 200 


300 


W - Secondary Flow Rate - lbm/ min 

s 


Figure 41 

Pressure Loss Characteristics 
Four Inlet Duct and Nozzle Cluster Configurations 

(Measurements Made at End of Short 
Mixing Tube^ L/d^ = 1, 35) 
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Pressure Loss - psf 


IjMllliffithtrii ti 


TOTuy .*.( u||iu 1 1 4|ii| 
: I : I ! ! i Mtiilit T 1 1i!it 


t(||:|)l|ji 'ijir {tr “ Stagnation Pressure 

jMSsBinlt !ii! Loss 


! I IM Mm !i 
f ! Mi li 


i-i- 
. 1 

ill 

-L-l 

! !!!l| 

i 1114)2 

■ 1 ■- n 

-it: 

iillii 

1 ' 1 I 

:3i 


111 Sii It I ; 

; : :ii tc: L:i 

-^"Ii jr-i iu; : 
rul:rt L.:i uM . . 


^ " : _ 


.:':t an; 

It; a; i i I ! 








O Static Pressure Change 
Corresponding to the 
Stagnation Pressure 
Loss 


60 80 100 150 200 300 400 

W g - Secondary Flow Rate - lbm/ min 


Figure 43 

Pressure Loss Characteristics of the 

Extended Constant Diameter Mixing Tube 

(Pressure Measurements Taken 
2. 72 Duct Diameters Apart ; 
d = 4.0 Inches) 
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Entrainment Ratio 
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22.5 
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18.4 
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15.6 


traverse probe 
direction 


Figure 44 

Velocity Profiles at Diffuser Discharge 

A 

d = 5.084" — t“ — - 1.71 

A 

m 

NAS 2-2518 Mixing Tube 
Original Nozzle Cluster 
Case 4 Nozzles (350 psia, 1090°F) 


W = 6. 69 lbm/min 
P 
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Axial Velocity fps 
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Velocity Profile with 
Heated Primary Flow 


Symbol 

© 

a 

Run No. 

7 

ii 

P po (psia) 

300 

260 

T (°F) 
po 

1158° 

1186° 

Wp(lbm/ min) 

5. 66 

4. 87 

W g (Ibm/min) 

143. 5 

144. 8 

m 

25.35 

29. 7 


Velocity Profile with 

Cold Primary Flow 

P = 191 psia 

po 

T = 83 "F 

po 

W = 6. 06 Ibm/min 
P 

W = 139. 00 Ibm/min 
s 

m = 23. 2 


in 



Local Radius Ratio Squared 


out _ traverse probe direction 


Figure 45 


Velocity Profiles at Diffuser Discharge 


d = 5. 084" 



1. 71 


NAS 2-2518 Mixing Tube 
Original Nozzle Cluster 
Case 4 Nozzles 
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Static Pressure At - Inches Water p ~ Static Pressure 0.81 Inches Down- 

Mixing Tube Discharge Gage stream of Primary Nozzle Discharge 

Inches of Water Gage 
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m - Entrainment Ratio 
Figure 46 

Static Pressure in NAS 2-2518 Mixing Tube 
Original Nozzle Cluster 
Case 4 Nozzles 


40 
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Measured Static Pressure: Inches Water Gage 
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mo 

from two Per- 
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NAS 2-2518 Mixing Tube 


® P = 350 psia 
po 


m - Entrainment Ratio 


Figure 48 

Jet Pump Pressure Rise Parameters Calculated at 
Mixing Tube Discharge 
Original Nozzle Cluster 
Case 4 Nozzles 
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Figure 49 

Velocity Profiles at Diffuser Discharge 

Aj 

d = 5.084" -r 2 - = 1.61 

m 

Short Constant Area Mixing Tube 
Original Nozzle Cluster - Case 4 Nozzles 
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Static Pressure At - Inches Water P - Static Pressure 0. 81 Inches Downstream 

Mixing Tube Discharge Gage of Primary Nozzle Discharge - Inches of 

Water Gage (tap #18) 


i 

ft 



Figure 50 

Static Pressures in Mixing Tube 
Short Constant Area Mixing Tube 
Original Nozzle Cluster 
Case 4 Nozzles 


p 

T 

po 

po 

psia 

6 F 

© 350 

1130® 

□ 300 

1160° 

A 260 

1165° 

^ 225 

1158° 
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3.0 4.0 5.0 6.0 7.0 


□ Original Cluster 
with Case 4 Nozzles 


A Original Cluster 
with Case 4 Nozzles 
and 2" Extensions 
(LD#1) 


* Low Drag Nozzle 
Cluster with Case 
LD#2-4 Nozzles 
(LD#2) 


Wp - Primary Flow Rate - 


lbm 

min 


Figure 51 

Secondary Flow Vs. Primary Flow 
Wide Open Throttle Condition 
Short Constant Area Mixing Tube 
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Figure 52 


Jet Pump Static Pressure Parameter Calculated at 
Mixing Tube Discharge 

Short Constant Area Mixing Tube 
Low-Drag Nozzle Cluster #1 
Case 4 Nozzles 
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Velocity Profile with 

Heated Primary Flow 

Run 41 

P = 300 psia 

po 

T = 1187 °F 

po 

Wp = 5.751bm/min 

W = 146. 9 lbm/ min 
s 

m = 25. 6 


Velocity Profile with 
Unheated Primary Flow 


po 
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po 

W = 6. 26 lbm/ min 

P 
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s 


m = 22. 85 



Local Radius Ratio Squared 
Figure 53 


out - traverse probe direction 


Velocity Profiles at Diffuser Discharge 


d = 5.084" 



1. 61 


Short Constant Area Mixing Tube 
LD#1 Nozzle Cluster 
Case 4 Nozzles 
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- Static Pressure a.t Mixing Tube P - Static Pressure 0.81" downstream of 

Discharge - Inches Water Gage primary nozzle discharge (tap #18) 


ft 



Primary Conditions 
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Tp 0 

(psia) 

°F 

O 350 

1160° 

□ 300 

1180° 

A 260 

1200° 

V 225 

1200° 


m - Entrainment Ratio 
Figure 54 

Static Pressures in Mixing Tube 

Short Constant Area Mixing Tube 
LD #1 Nozzle Cluster 
Case 4 Nozzles 
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Figure 55 

Variation of Static Pressure in Mixing Tube and Diffuser 
for Three Nozzle Clusters 

P = 350 psia 

po 

T = 1100° F to 1200° F 

po 

W m = 124. 9 lbm/min 
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Analytical Results 

Experimental Results 
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Figure 56 

Jet Pump Static Pressure Parameters Calculated at 
Mixing Tube Discharge 
Short Constant Area Mixing Tube 

Low-Drag Nozzle Cluster #2 
LD #2-4 Nozzles (7 Nozzles) 
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Figure 57 

Velocity Profiles at Diffuser Discharge 


a h 
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Short Constant Area Mixing Tube 
LD #2 Nozzle Cluster 
LD#2-4 Nozzles 
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Static Pressure at Mixing Tube p _ static Pressure 0.81” downstream of primary nozzle 

Discharge - Inches Water Gage discharge - Inches Water Gage (tap #18) 
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Primary Conditions 

P T 

po po 

(psia) °F 

O 350 1130° 

□ 300 1195° 

A 260 1190° 
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m - Entrainment Ratio 

Figure 58 

Static Pressures in Mixing Tube 
Short Constant Area Mixing Tube 
LD #2 Nozzle Cluster 
LD #2-4 Nozzles 
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Figure 61 
Velocity Profiles 

Extended Constant Area Mixing Tube 
LD#1 Nozzle Cluster — Case 4 Nozzles 


189 


- Axial Velocity - fps 


t 



500 


400 


300 


200 


100 


Traverse 1 


— □ — Traverse 3 



400 


300 ^^ 


200 fc 


100 


(r/r,J - Local Radius Ratio Squared 

w 


.8 1.0 

out - traverse probe 
direction 


Figure 62 
Velocity Profiles 

Extended Constant Area Mixing Tube 
LD #1 Nozzle Cluster - Case 4 Nozzles 
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Figure 63 


Measured Static Pressure - Inches Water Gage 



m 

Run No. 

17.3 

99 

19.2 

98 

22.5 

97 

26.4 

96 

P 

po 

= 350 psia 

T 

po 

= 1112°F 


Suction 

Tube 


V 


Mixing Tube Extended 


Diffuser 


r 


Figure 64 

Variation of Static Pressure in Mixing Tube and Diffuser 
Extended Constant Area Mixing Tube 
LD #1 Nozzle Cluster - Case 4 Nozzles 
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Figure 65 


Jet Pump Static Pressure Parameter Calculated 
At the Mixing Tube Discharge 
Short Constant Area Mixing Tube 
Ljj #1 Nozzle Cluster - Case 4A Nozzles 
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Figure 66 

Velocity Profiles at Diffus.er Discharge 

d = 5. 087” 

A/A = 1.61 
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Short Constant Area Mixing Tube 
LD #1 Nozzle Cluster - Case 4A Nozzles 
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Figure 67 


Static Pressures in Mixing Tube 
Short Constant Area Mixing Tube 
LD #1 Nozzle Cluster 
Case 4A Nozzles 


195 


Static Pressure Parameter 


ft 

i 


o 

cn 


ft 



m - Entrainment Ratio 
Figure 68 

Jet Pump Static Pressure Parameter 


Extended Constant Area Mixing Tube 
LD #1 Nozzle Cluster 

Case 4A Nozzles 


196 



Static Pressure Parameter 


Analytical Results 
Experimental Results 



Figure 69 

Jet Pump Static Pressure Parameter 
Extended Constant Area Mixing Tube 
LD#1 Nozzle Cluster 
Case 4A Nozzles 
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Figure 70 

Static Pressures in Extended Constant Area Mixing Tube 
LD #1 Nozzle Cluster - Case 4A Nozzles 
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Figure 71 

Variation of Static Pressure in Mixing Tube and Diffuser 
Extended Constant Area Mixing Tube 
LD #1 Nozzle Cluster 
Case 4 A Nozzles 


T 

Run ^ 
143 1140 

142 1158 

141 1162 

140 1175 

= 350 psia 


199 



Static Pressure Parameter 



m - Entrainment Ratio 
Figure 72 

Jet Pump Static Pressure Parameter 
Extended Constant Area Mixing Tube 
LE>#1 Nozzle Cluster - Case 2 Nozzles 
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Figure 73 
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Figure 74 
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Static Pressures in Extended Constant Area Mixing Tube 

LD #1 Nozzle Cluster 
Case 2 Nozzles 


Measured Static Pressure - Inches of Water Gage 
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Figure 75 


Variation of Static Pressure in Mixing Tube and Diffuser 
Extended Constant Area Mixing Tube 
LD #1 Nozzle Cluster - Case 2 Nozzles 
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Figure 76 

Jet Pump Static Pressure Parameter Calculated 
at the Mixing Tube Discharge 
Short Constant Area Mixing Tube 

LD #1 Nozzle Cluster 
Case 2 Nozzles 
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Figure 77 

Static Pressures in Mixing Tube 
Short Constant Area Mixing Tube 
LD #1 Nozzle Cluster 
Case 2 Nozzles 
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Figure 78 

Jet Pump Static Pressure Parameter 
Extended Constant Area Mixing Tube 
LD #1 Nozzle Cluster 
Case 3 Nozzles 
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Figure 81 

Variation of Static Pressure in Mixing Tube and Diffuser 
Extended Constant Area Mixing Tube 

LD #1 Nozzle Cluster 
Case 3 Nozzles 
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Figure 82 

Jet Pump Static Pressure Parameter Calculated at 
Mixing Tube Discharge 


Short Constant Area Mixing Tube 
LD #1 Nozzle Cluster 
Case 3 Nozzles 
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Figure 83 


Velocity Profiles at Diffuser Discharge 
Short Constant Area Mixing Tube 

d = 5.084" A, /A = 1.61 
d m 

LD #1 Nozzle Cluster 
Case 3 Nozzles 
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Figure 84 

Static Pressures in Mixing- Tube 
Short Constant Area Mixing Tube 
LD #1 Nozzle Cluster 
Case 3 Nozzles 
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Figure 85 

Velocity Profiles with Unheated Primary Flow 
Short Constant-Area Mixing Tube 

LD#1 .Nozzle Cluster 
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Figure 86 

Influence of Duct Loss Coefficient Upon 
Predicted Jet Pump Performance 


Case 2 Nozzle, Extended Mixing Tube 
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